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§ 1. Introduction. 

The present paper continues my researches in the theory of gamma functions. 
Previously to a certain extent I obtained known results by new methods : none of 
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the sricceediiig investigations however have, I believe, been undertaken or suggested 
by other mathematicians. 

In the first paper'"^ published in the connection I, attempted to give a homogeneous 
theory of the ordinary gamma function, considered from the point of view of 
WErERSTRASs' function theory. I introduced a parameter o), and showed that the 
theory was subordinate to that of a function satisfying the difference equation 

s being any complex quantity. 

That theory led naturally to the considerationf of the G iunction, satisfying the 
difference equation 

G{z + 1) ^ r(^)G(^), 

and substantially a function all of whose properties coidd be obtained by differentiating 
the simple gamma function with respect to the parameter, 

I next considered^ an extended function G{z/t) satisfying the two functional 
relations 

f{z + 1) - r^''^)/(^) ; f{z + r) = r(.) (27rp' r-^VX^), 

and reducible to the G function when r = 1. Several points in that paper suggested 
the formation of a symmetrical douhle gamma function, in which r should be replaced 
hy the quotient of two parameters w^ and w^^ In the present investigation such a 
function is defined, and its theory developed in, I hope, complete detail. The 
function is the natural extension to two parameters of the simple gamma function 
V,{z I 0)). 

It is necessary for a complete exposition, of the theory to consider the properties of 
wliat I propose to call double BernouUian numbers and functions : functions which 
are the natural extension to two parameters of the simple BernouUian functions, 
considered in Part II, of the earliest paper of the series^ Such a theory is develo|)ed 
in Part I. of the present paper. 

In Part 11. I consider the elementary theory of the double gamma function. It 
is shown that certain symmetrical modular constants arise as finite terms of 
asymptotic expansions in a manner exactly analogous to the origin of Eulek's 
constant y. 

Such considerations lead naturally to Part III., in which are deduced from a 
contour integral, which is a double generalisation of Riemann's £ function, certain 
noteworthy asymptotic approximations, of which the most important is an extension 

'"" Bi^RNES, ''The Theory of the Gamma Function," 'Messenger of Mathematics,' voL 29, pp. 64-128. 
t Barnes, " The Theory of the Q Function," ' Quarterly Journal of Mathematics,' vol. 31, pp. 264-314. 
X Barnes, " Genesis of the Double Gamma Function," ' Proceedings of the London Mathematical 
Society,' vol. 31, pp. 358-381. 
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of Stirling^'s theorem. By the aid of this theory it is possible to express the 
logarithm of the double gamma function and the double gamma modular constants as 
contour integrals similar to those given in Part III. of the ' Theory of the Gamma 
r unction. 

In Part IV. I consider the multiplication and transformation theories of double 
gamma functions as well as certain curious integral formula), which correspond to 
Raabe's theorem for the simple gamma function, and are elementary cases of a 
general theorem connecting successive similar transcendents of higher orders. 

In Part V. the asymptotic expansion of the double gamma function is obtained, 
and it is shown that the function cannot arise as the solution of a differentia.1 equation 
whose coefficients are more simple transcendents. 

There exist similar functions of any nu.mber of parameters, and these transcendents 
I propose to call multiple gamma functions. I reserve the formal expression of their 
properties for publication elsewhere. I have worked out the theory for double 
gamma functions independently inasmuch as, the complex variable being tw^o 
dimensional, there are many points in which a higher analogy breaks down : and also, 
since many proofs in the higher theory are, in their simplest form, inductive and, 
to be rigorous, require a knowledge of the theorem, for the two simplest cases. Not 
only so, but in the case of the double gamma functions it is possible to give easily an 
alo'ebraical theory (such as that worked out in Part 11. ), which is more simple than 
if one derived all the formulae from the fundamental consideration of certain contour 
integrals. 

I append a statement of the notation adopted in this paper, mentioning the place 
in the present series of investigations where such notation is used for the first time. 



'Derivation. 


Name. 


Symbol. 

Bn(cf-' to) 


First occurrence. 


Algebraic solution of 

f{a + «) - /(«) = a» 


Simple Bernoullian func- 
tion 


'' Gamma Function," § 11. 


7h 


Simple Bernoullian num- 
ber 


AM 


"Gamma Function,'' § 15. 


Algebraic solution of 

f{a + oji) ^ f{a) 

== S. {a (02) + ^'^^±^ ^^ "^ 
^ n+l 


Double Bernoullian 
function 


2Sn{a (01, OJ2) 


For the case of equal para- 
meters : 
" G Function," § 15. 

In general : 

" Double Gamma Function," 

'^' 


2B'n(0 0)1, 0)9) 

n 


Double Bernoullian num- 
ber 


2B,j, + l((0i, CO2) 


" Double Gamma Function," 


^-1 J'-\ 

(0 W l\(jy> 


Simple gamma function 


T,{z (o) 


" Gamma Function," § 2. 




- 
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Derivation. 






]ogri(,^|w) 



Solution oi j\z + 1) :== V{z)f{z) 
Solution of /(;^ + 1) - Tl^)f{p) 



Vide §§ 18-24 



-^-^ log V^{z I (01, wg) 



Name. 



Symbol 



Logarithm derivative of 
simple gamma function 



^i^'^'^C^ I ^0 



G function 



(•'^) 



Unsymmetrical double 
gamma function 



Double gamma function 



Logarithm derivatives 
of double gamma func- 
tion 



tr(l -^ s) 



m 



e 



-az 



1 



/9^(o:, 



-- ( " zy-'^dr. 



27r JCI ~-e^^y'')(l -e~^ 



Simple Riemann (' (zeta) 
function 



L 



(1 -^"^'^i'^)(l -e~^^) 



Double Riemann {" func- 
tion 



1 + !m_^ 



CO 



0> 



Finite terms of certain asymp- 
totic limits 



Do, 



do. 



Do. 



do. 




27r 



to 



Limit of a certain definite in- 
tegral 



Constants which take the values 
Oi ± 1, according to the dis- 
tribution of 0)1 and 0)2 



Simple gamma modular 
form 



Unsymmetrical double 
gamma modular forms 



Symmetrical double 
gamma modular forms 



Glaisher-Kinkelin con- 
stant 



Simple Stirling modular 
form 



Double Stirling modular 
form 



G{PJ I r) 



^2(2^ I wi, (.02) 



^2*'^('^^|t'^l, W2) 



First occurrence. 



" Gamma Function," § 2. 



a n. ly 



G Function," § 3. 



u 



Genesis," &c., § 1, 



" Double Gamma Function," 



C(s,a,o)) 



^(s, a I (Di, C02) 



rn(to) 



C(r) 
I){r) 



721(^1, C02) 

722(^1, 0)2} 



pi(co) 



P2(wi, OJ2) 



* « 



m 



\i 



} 



u 



Double Gamma Function 



3? 



'^ Gamma Function," § 23 



For ecpial parameters : 

'' G Function," § 23. 

Lr general : 

'' Double Gamma Function," 

§39. 



" Gamma Function," § 2. 



a 



Genesis," &c., §§ 3 and 4. 



" Double Gamma Function," 
§§ 21 and 23. 



a 



G Function," § 3. 



^' Gamma Function j" § 31, 



" Double Gamma Function," 

§ 4a 



" Double Gamma Function," 

§2L 
" Double Gamma Function," 

§39. 



The symbolic notation by which F,2[f{z + w)] is written for 

f{z + «i + 0J2) ^f{z + co^) -f{z + m^) 



IS introduced in §49, 
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The Theory of Dotihle Bernoullicm Functiofis and Niimhers. 

§ 2. In the " Theory of the Gamma Function," Part IL, we have defined the simple 
BernouUian function S,,(a | o>) as that solution of the difference equation 

f{a + (^) -- /(^) = <^^'\ 

where n is a positive integer^ which is such that it is an algebraical poljniomial and 
S,^(o I ct}) ==: 0. And it was proved that such a solution does exist. 

In exactly the same manner it may be proved that the difference equation 

has an algebraic solution, which is a rational integral polynomial of degree n + 2. 

The difference between any two solutions will be a simply periodic function of 
period w, and will therefore be a constant if the solutions are both algebraic 
polynomials. 

There thus exists a unique algebraical polynomial of degree n -\- 2, which is a 
solution of the difference equation. 

f{a + CO , ) »^ /(a) rz. S,{a \ m,) + ^ ""^ p , 

vnth the condition /(o) = 0. 

This solution we call the double BernouUian function of a with parameters o)^ and 
0J2, and we denote it by 2^^,^{a\(i)j_,o)^). By symmetry with this notation the simple 
BernouUian function would be denoted by xS,,(a | o)). 

We shall often omit the parameters m^ and m^, when there is no doubt as to their 
existence, and write the function simply 2S;^(a). 

§ 3, We now proceed to show that the double BernouUian function of a of order 7i 
is also the unique algebraical polynomial which is the solution of the difference 
equation 

J- [a + oig) -J {a) rr h,{a 0)^)+ --^;"_^-Y' , 

with the condition y(o) =: 0. 
For since 

we have at once 



w 



Bn{ct' + coo j ft)o) — S„{a I «o) "= ct" = S.;,(a +6) J cd J — S^,(a ] ojj). 
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And therefore if we pvit 

/,(«) ~ 3S„(a -f 0J3) -■ gS,(a) — S,(a 



'1) 



we shall have 



/.{a + 0)^) —fn{a) ^ 0. 



Now 2S«(a) and S,4a/o)|) are algebraical polynomials of a, and therefore /;,(a) is 
also such a function. And^ therefor e, since it is simply periodic of period o),, it must 
be a constant. Thus we have 



,S«(a + (jt)^) — 2S,,(a-) == S,,(ajcfjj) + constant 



o 8 a o 



I ^. ) 



Again, integrating the relation (l) with respect to a between and w^, w^e have 



f COi +W.) 







2^ 



•^0 -'0 " 



smce 



'Wq 



D 



Sn{€l I (D^da =: — COg 



72 H- 1 



Integrating the relation (2) in the same manner between and ojj, w^e obtain for 
the value of the constant 






a\(j){)da 



^'n.yi{0\0)i) 



And thus oS;,(a) is the unique algebraic solution of the equation 



f{a + dig) ~/(a) :r= S;,(alw|) + 



S^ 



71 + t 



with the condition /"(o) ==: 0, 

From the symmetrical nature of the equations which give c^n{a } ci)^, Wo), wo see that 
this fimction itself must be symmetrical in m^ and cy^. 

§ 4. If now we assume 



2S,,(a|oji, C03) = oe,,^,2a^^+2 ^ ^^^^^^^«+i ^ 



-4° Ct-^Cf,-, 



the calculation of the highest coefFicients may be readily effected. 
For we have 

n -f- 1 



S«(a + c/)j) — 9^n{G() ■= ^,i{a I a)^,^ d 



a 



2 

M + l 



(72+ l)<i)o 



a''^ . /''/2 \ I3i , 






see 



Hence if we substitute the assumed expansion for oS^,(a) and equate corresponding 
powers of a, we find 

1 

{n + 2) £o^a,.^3 = 



(72 -h 1) &)o 



(72 + 2) (72 + 1) 



1 9 



<^/' ^^M~3 + (^'^ + 1) ^1 ^n-Vi 



*> 5 
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(n + 2)(n+l)n „ , (n+l)n „ fn\Bi 
— o)y^ ci„+ci + oj/ «„+] + nwj a„ = j—o). 

and so on. 
On solving these equations successively we readily obtain 



3? 



^;/ + 2 



^//-M 



and, since B, =: i 



-L e^ 

2(r^ -1- 1) W^Wg ' 



2 1 O 



^)i — 12 



y COj + tOg'" + OCOiCO^ 



Thus 






^ ^ 1 ^' ^^ Oi-fl) (r^ + 2)coiWo 2(?i. -f l)*^^^, 12wiWo ^ 



Further terms can be calculated if necessary. It will be seen, however, that they 
form what we propose to call double BernouUian numbers, v/hose properties may be 
investigated without the necessity of their formal evaluation. 

Corollary, We note that 



And hencc; 



^Pi{a\o}i, Mc^) ™ ,r -— —-- - - -^ a - - - 

0p ^[a o)^, 0J3) :— ~ -^ ■^~- H -^ —^ , 

S/c,\/ \ \ ^^ Ct)i 4" (J^o 
^ ^^^\a 0)1, Wg) = — " " , 

WiCt)o 

It will be found that these expressions are of constant occurrence in the course of 
the present investigation. 
Note also that 

2(o^o)o 2ft)i«o 

§ 5. We will now prove that, if n — k^O, Jc > 0, 

We have, M'hen n — /c > and h < 0, 

,S„W (a + coO - gS/) (a) = S,/« {« I 0,2) 

irr: S,/^^^ (0 I co^) + 7 yr: S,,_i.(a | w^) 

('' Theory of the Gamma Function,'^ § 14). 

VOL. CXCAa, — A. 2 N 
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and therefore 

aS^W (a + CO,) - ,S/) (a) 

= S/)(o|c.2)+ T~-^b«-.(a + coi) -™ ,S,„,(a) ^ ^^i:=^±M^) 

^ ' ^ (71 — k)l ^ ' ^ ^ n— k + 1 

Thus if we write 

f(a) = ,S,;^0 (a) ^ -^- ,S,^, (a) 
we shall have 

= (^^ Gamma Function/^ § 15). 

Similarly 

f{a -p coo) — f{a) =:= 0. 

and therefore since /(a) is an algebraical polynomial in a, it is a constant. 
On making a = and remembering that ^n{o) =^ 0^ we obtain 

.S«^'K«) - (7^ .S«_. («) =. ,S/) (o), 

which is the required result. 

§ 6. We are now in a position to prove that 







,S„ (a) c?« = - ~f^ ,S'„,, (o 1 «i, «,) + ^^^^ _^_ ^^ _ ^^^^ ^ 2^ 



^ Ci / \ 7 ^^ C\f /I \ 1 ^ ?J-f V^ ^ij 

and at the same time the important relation 

Sn\/\ \ yiv *T' 'i') rN / / \ 

. n + A: (0 1 Oil, Wg) == /^^:^~3y ' 2^ nvA^ ^1. ^2)- 

Since 28^^(0) ™ we see from the fundamental difference equation that 
Hence, if we put a =:^ w, in § 5, we see that when n--h > 0, and h > 0, 



Take now the relation 

(7? ^" ^o 



,S/>(a) z^,;^--^-^^,,,S,.,(a) + 2S/)(o). 
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and integrate with respect to a between and coi, we obtain 



,S/-i^ M - ,S, «-!> (o) 



n 



0}\ 



(71 -h)] Jo 



B,,^k («) da + 0)1 ^S/^ (o), 



so that 



n I 



rtoi 



(?^ - /^ 1 . 



(J 



*pfi-^k \(y del 



n I 



^1 sS/^ (^) + (,, :|:2i^^) I ^'^^^■3-^(^1 ^s)- 



Write now rt. for {n—h), as is evidently allowable, since both n and I^ are positive 
integers, and we have 



•wi 



n \ 







..r, ^Slfl, (o) + 



b u\-2 \P I ^'3/ 



3S, (a) ck --- -»- coj ^^^ _^^ ^^^ ^ ,^,+, ^.; , ^^^ _^_ ^^ ^ ^^^^ _^_ 2)- 



We thus see that 



71 I 



^nlk (^1^1? ^s) 



is independent of ^, since this is the only time in the relation just obtained which 
depends on h 

Putting then ^ = 1, we have when ^ > and n > h. 



(m + 1) 1 



which is one of the relations required. 
And also 



Wi 







gS;, (a) da =r: 



ft) 



^i + 1 •' 



?S^+i (o) + 



8'« + 3 (^ I ^2) 



o\> 



(t?. + 1) . (n + 2) 



another of the given relations. The second integral formula of course may be 

written down by symmetry. 

We notice that in the notation formerly introduced ('^ Gamma Function/' §15) 

we have 

S'm +q (q I C03) __ iB»+2 {(^2) 

(^ + 1) . (^^ 4- 2) """ ■ 7h + 2™ ' 



and therefore that each of these expressions 







n 



""" (?^ -f 1) . (^^ + 2) 

Thus we see that when n is odd 



when n is odd. 



when n is even. 



w, 



ft)l 







3S. («) ^^^ = - :f^ S'n^i (0). 
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§ 7. We now inti-oduce double BernouUian immbers analogous to the simple 
Bernoiillian numbers introduced in the theory of the gamma function. 

In that theory the simple nth BernouUian number was defined by the relation 

__ S^n{0 I ft)) 

and now the rith double BernouUian number is given b}^ 



B,, (w) 



^B,, {a)i, 0)^) 



n 



We note that by the theorem of § 6 we may put 



n ' 



jB«+i (^n ^2) — Z: r'iA\ 3^^^-^ (^ 1 ^1' ^2). 



and therefore 

71 ! 






§ 8. At this point we may conveniently note the reduction which takes place in 
the double BernouUian functions and numbers when the parameters are equal to one 
another. 

If we put mi = 03.^=^(0 we have as the single difference equation of the ^ith double 
BernouUian function the relation 



f{a + m) —f{a) =:: 8u{a\a)) + -Ji,,^^{a)), 

and. the function is now defined as the algebraical solution, of this equation with the 
condition 2^,^(0 j co^ co) == 0. 
Put now 

f(^a) ::= .— S;, ^.^ (a 1 0)) + (« ^ «^) S,,,(a | m) + a ' "' '-^^- 5 



and we have 



t[€C + 0)) ■— r{a) = (Dbj\a (i)) -\- T , 7 "- 



Hence, the other definition conditions being satisfied, we see that 

Snialoj, oj) :=:^- — b;, (a U)) — b„.,^ (a ft)) + -^- r"n^? 

that is, the double BernouUian function when the parameters are equal reduces to 
simple BernouUian functions and numbers. 

It will be seen later that it, is for this reason that it was possible to obtain aU the 
expansions in the theory of the G function in terms of simple BernouUian functions, 
ote that the above relation may be written 

"1 

cS;;(a ft), ft)) = — — - S;, (a I ft)) •— - S;m . (a ft)) + "- iB,^+i (o)). 

■^ ^ ' CO ^ ^ ^ CO 0) 



DOUBLE GAMMA FUNCTION. 



277 



On differentiation we have 



-I -j 



so that 



and therefore 



gS',, (O I COj O)) ==: — S',, (O I Oj) — - ^,,.^1 (w), 



2B,, (o), Oi) 



CO 



§ 9. We now see at once that 



Sn {Cl 1 <i>i5 CO^) 



« 



?H-3 



0)1 -f- ft)o , 1 

. __ ^ _ ,__e ^ ^fH-t-] 

(n +1) (n + 2) 6>]^ coo 2 (n + 1) «! w.^ 



<^'^'' + A(^o^2)^^' 



J 2B, (o),, 0;,) a-i + y 2B3 (o),, o),) a-^ + 



and so complete the expansion of § 4. 
For by Maclaurin's theorem we have 



^^u 



(ajcoi, W3) = %S';, (0) + 



2^w, "^ V^j 



2! 






shice the hiofher differentials vanish. 



From the few terms found in S 4 we see that 



^ {n + 2) ! '' ' 



Now when ^2 > k and />; > 0, we have 



n I 



m^o)^ 



n ! (ft), + ft)p) 



2Sf(o) 



?l 



(7^ - ZO 1 



^JBn^k^^l{^\, (^2)^ 



and thus we have the expansion in question. 
§ 10. It may now be shown that 



(") 



M-l 



,S„(a) = {-)" 2S„(a>i + CO,- a) + ^-^Y [S;,, (ojci) + S; + i((>|co,)], 



or, as we may write it, 



Remembering the value of, iB,,+i(a)) we thus prove that 

^S,, (a) =3: ^S,, (wi + 0)3 •- a), when n is even, 

K- 1 

^ ^ ^ B^^+j^{o){' + 0)/), when n is odd. 



2S,, (a) ^ — ^Bn (a)i + 0)3 — a) + 



?^ -{- 1 



2 
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Substitute w^— a for a in the fundamental difference equation and we find 



Sn{0)l + % — a) — 28,4^3 "~ ^0 = S,,(0).3 



a 



%) + iB/^ + i(%)= 



and therefore, since S;,(o)2 — ajoi^) :==z (»-)«+i ^,,[a\M.^), we have 



Sn{0)i + % -- ^0 



glO;^ I U^3 



(oi^ ^ a) := (^)- ' S,(a|a)o) + iB, + i{co^). 



If therefore we put 



n 



oS« (ft);^ + 0)^ — a) + -^ i-fi + j (^^9.) 



+ T> / \ 



we see 



that f{ci) is an algebraic solution of the difference equation 

and therefore can only differ by a constant from ^nifl)^ 

Determine this constant by making a = and we have the relation 



sS^a) ™ (-»,)« [oS,,(cui + % »™. a) — .S,,(co_^ + 0)3) 









When n is odd the last two terms cancel each other, and when % is even ^^,^^{0^^ 
vanisnes. 



J Lt/il-v_/t5 



%^ 



S;,(ala)i, 0).) ™ ( — y^ [38,(0), + 



^3 



a) -- Sn{(^i + %). 



From the fundamental difference equations we see at once that 






n -\- 1 



n + 1 



l^u-hli^l) + l^n + l{0)%): 



and therefore we have the relation stated. 

§ 11. We may now show that, when n is even, 






2'?i 



a simple expression for the even double BernouUian numbers which corresponds in 
some degree to the fact that the even simple BernouUian numbers vanish. 

On differentiating with regard to a the result of the previous paragraph we find 



2 



Wn{^\o)l, 0)o) + { — y.S'ni^l + % ™ ^0 ^ ^* 



Fi-om the fundamental difference equations we have 

^Bn (a + wi + 0)^) — 2,8;, (a + wi) »»• gS,, {a + w.^) + ^8,, (a) = a", 
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and hence ^S^^ (coj^ + cog) == S\i (<^i) + 2^^n{o>^) •— S^n (o) 

Thus, since ('' Gamma Function/' § 15) 

S',, {o\o)) = when n is odd 



n 



( — ) ''^ B,; ct)T^* ^ when ^^ is even. 



2" 



we see that 



gS'^^ (ojj^ + ^^s) = 2^'u (o), when n is odd ; 



''-1 



and c^8\,{a)i + Wg) = ^S',, (o) + ( — )' ^^{^i'' ^ + ^i' ^)^ ^^en 7^ is even. 



But our former relation gives us, when n is even, 



2 



S'„ {a>, + cu.) 



2 



S'. (o). 



Hence, when n is even, 



'rt I \n 



B 



^8',,{o) = {~f .-,/ . «-l + <-!), 



/J 



which is equivalent to the relation required. 
§ 12. We now proceed to show that 



aS,(a) da = V"^T ^ ^' 2^^/+] (^1? ^2)- 



We have 







n -j- 1 



oS,, (a + ftjj ) — gS,, (a) =r: S^^ (a \ ojg) + 



n -]- 1 



Hence, integrating with respect to a 

'a + 0}} 







2 



S^, (a) da 



a 







2 



S;, (a) cfo. 



fWi 







2^ 



a 



S,, (a) 6^a + Su (a \ ojg) da + <^ 



71+1 



But (''Gamma Function," § 19) 



so that, if /(cf) 
equation 



" S, (a 1 0),) da+a ^^^M^ -.^ ^±i(^ K) . 
2S,, (a) da^ this function is an algebraic solution of the difference 



•' 



f(a + CO,) —/(a) 



gb^; (^aj aO' + ■ 







n ~\- 1 



But this diflPerence equation is evidently satisfied by 



a 



^ „ (a) da — -; — ^-; V "N 
«i LJo ^ ^ (7i + 1) (71 + 2) _ 



1 



+ ,r+ 1 2^"+i ^") 



re 



CDi 



.. .!_ 1 2P^^ + 1 (^) 



'?^ f 1 



+ "T^ by § 6, 



r^ h 1 



n H- 1 



a gB^^^j^i (ojj, fti.^). 
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Hence as these two solutions both vanish when a := 0, we have 



•a 







^S,, (a) da = ^J;J^ ^" — a ^B,.^^^ {m^, co^). 



As a corollary we have on differentiation 

gS',,,.^ (a) - (n + 1) ^S, (a) + ^8',_,j (o). 
§ 13. Tiie multiplication theory of double Bernoullian functions may be con- 



veniently expressed by the formula 



m"' gS,, ( a 



C0| (Oc, 



) 



From the fundamental difference relation we have 



2 



S,, -<^ m I a + 



COi 



Qll 



6) 2, CO 2 



gS,, {ma 



oj^, r/jg_j 



'^ Wh 



11 



s, 



O;^ I a 



u 



COo 

A/ 

7]i 



Hence ^7/2^^/^ {ma\ w^. oj^ satisfies the difference equation 

Jlv 



fia-i-fj-fia) 



H 



a 



\ 



0)0 



rih / ^ r i 



/ 



and is the onhy algebraical solution such that/(o) =^ 0. 



Hence 



1 



^; gS, (ma I o)^, ojg) ^.BJa 



m m 



OOc 






'lb + 1 -J 



which is the relation required. 

As a corollary we see that the ^ith double Bernoullian number is homogeneous 
and of degree (^—1) in the co's. 

For in \ 9 we have seen that in the expansion of 

the part of the coefficient of a^^~^^^ which involves the w's is oB,,(<^j, Wg). 

§ 14. The transformation of the parameters of the double Bernoullian function is 
given by the relation 



gS;, ( a 



^1 ^3 



p~lq~l 



t t gS,, (a + 



^==0 Z,s=0 






0) 



, ft)g) 



+ pg^B^Hi (^^i. ^2) -^ 2-B: 



ft)i tO: 



;/4 1 



i> ' '/ " 



as we proceed to prove. 
Let 



f{a) 



S S 0^5 /^ (a "4^ 4^ • 



0)j, CO,), 



then we obtain at once 
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q~-l 
1^0 



I 



Wr 



SJ a + -- + ., 



S J a + 



'&)<: 



^ 



O ( CI ~p 



/ft).-, 
9. 



% I + lB;. + l(0J2) 



si a 



CO, 



9 



+.^".(?)(-— — "^-^ 



Thus f{a) is an algebraic solution of the difference equation satisfied by 



0/|( Cv 






The two solutions can then only differ by a constant, and thus 



2^«> Cfr 



ft)l fi)c? 









where R,, is independent of a. 

To determine this constant, let us integrate between and ~-. We find 



^ 



ft) I ft)Q 



c^a 






«=0 







'i^?4 b 'I 



/ft); 



CO 






_p 



and therefore by § 6 



ft) 



^1 %v I 

^ '-+11 .r? -~. I "T 



1 "D 

P \P '1 



P / ^3\ 

1 •*•*« + 3 

'?^ + 1 ' 



^^1 sB^+i (<^i. %) + S 



g-isS-M+if «i + 



ft) 



i^ 



9P/. 



7ft).^' 



<+li 



ft)l 



1=0 



n -f- 1. 



^/ — /_ _L. , j[^^^ j-^Y b 1 ^. 

p 



t> / \ I / -^ ^Nl-I^^w+S (%) _j_ ^ i-t-%+2 (ft>g) _| *^l-p 



S' 



p 



And therefore 



(** Gamma Function,'' | 18). 






ft) 






so that 



R» = M.B„,,K-,)-,B,.,,(^\'^^ 



On substituting this value we have the theorem enunciated. 
As a corollary we have on making a = 0. 



V 'S^ Q I 1 .... {. , — ? 



^=o?=o \ i? 



^ 



0) 



.? ^3) — s-^« 



mi co^ 



S-^« + l .. ' 



P 9. 



pq aB,,+j (ft)i, m.^), 



§ 15. We now proceed to prove the expansion of fundamental importance in the 
theory of double BernouUian functions : — 



ze ^^^ 



(1 ~ fr'«i^) (1 -- 6~^a^") "^ .<s^ 



^S,<^)(a) + — V^ ^^ 
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In this expansion the double BernouUian functions have co^ and o)^ for parameters, 
and the expansion is vaKd provided 1 2; [ is less than the smaller of the tAVO quantities 



2nTi 


• 


27n 


'"i 



For within a circle whose radius is less than this quantity, the function 



.^2Q-az 



(1 — e"''^) (1 — er^'^ 

has no poles, and hence it is expansible in a Taylor s series of powers of z. 
Thus we may assume, for all finite values of a^ 

,,-^- - "^; ^ = ^^ - Ao(a) +^^.z + .. .+ (_)»-i.^i^f),« + . . . . 

(1 — G""^^') (1 — e-'^^) z ^^ ^ ' 1 ! ' » \ / ^^ t ' 

where it is obvious that the functions of a which enter as coefiicients are all alge- 
braical polynomials. 

Change now a into a + ^^i, and subtract the expansion so obtained from the one 
just written. 

We find 



/C 



.-«^ _ A,(a) - A,(. + .,) ^^^^^ ^ ^^^^^ ^ ^^^ ^ m^±^ , + . _ 









But in the '' Theory of the Gamma Function," § 20, we obtained the expansion which 
may be written 



S 






1 -_ e-^2^ 

Equating coefficients in these two expansions, we obtain 

k^{a -f 0)^) — Ai(a) ^ 

Kj^a + ^1) — Aq(c6) =• '^^^^{a\(ii^ 



««*asBs«9eaeeeeno* 



and it is obvious that a similar set of equations hold in which (o^ and coo are inter- 
changed. 

Hence A^{a) is a constant whose value, from the first term of the expansion, is 



^1^3 



3 



S^^^\a\co^,(o^) by §4. 



Again Ao(a) ~ gS/^^(a]a)^, m^), for these two expressions can only differ by a 
constant which by § 4 and the actual expansion is at once seen to vanish. 
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Finally for all positive integral values of n, fu{ci) can only differ from 2S'«(a | cd^, cd^) 
by a constant,— ~^.^. 5 let us say, 

Differentiate the expansion thus obtained with respect to a, and we find 

(1 — e"^i^) (1 — e-"^0 "^ ^ ^ ^ 1 ! • \ / ^^ I » 

and hence we have finally (§6) 



(1 — 6~"i') (1 — e~"2') ^ xJ 1 \ / I J^ I I • • • I V / ^^ I 

which is the expansion required. 

This expansion may be used to define the double BernouUian numbers, and all their 
properties may be deduced from it. A procedure analogous to the one here suggested 
will be the one employed in the general theory of multiple BernouUian functions. 

§ 16. Several expansions of constant occurrence may be deduced from the one just 
obtained. 

In the first place, note that we may write the expansion in the form 

ze~^^ 1 . ft). + ft>Q <^ , oS'i (a) 



(1 -— ^~'^^^)(1 — e""*«^22) CO^CO^Z 2ft)^ft)2 ft)^ft)2 1 1 

^ ^ nl 

Put now a = 0, and we have by definition of the double BernouUian numbers, 



(1 --6-'^i^)(l — e-'^'^') co^co^z • 2ft)ift>3 



^ ^ (71 -- 1) ! 



We thus have (§11) 



1 , «! + ^3 , sl^l (^1' ^3) . I sl^S (^1' ^2)^3 



+ -~^-^ + ^-^^- ^ + ^-^Vf-- '^' + • • . 



(1 _ ^-a>i.)(i _ g-cu,.) ft)ift).32: ' 2ft)ift)3 0! 2! 



(2?i) ! 



"^„.t/ ^ 2. (2m)! 

And the last series may be written 



1 

2 



z z 11 

1 ^ g - <0i^ « J _ ^ - 0,2^ ^^ ^^ 



Hence we find 

2 o 2 
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Z -, z z z 



as the expansion from which the odd double BernouUian numbers may be derived. 

Finally if we integrate the -fundamental expansion of § 15 with respect to a 
between and a, we obtain 

„ '■ -,V- , = ^'--^^ -''■""■> - [.S',(c.) - .S',(o)l + &?'^ + . . . 



(-)»-\S„(a) 

n ■ 



or, as we may write it. 






as the expansion from which the double BernouUian functions themselves are at once 
obtained. 

All these expansions are valid within the circle whose radius is the smaller of the 

two quantities 



tTVl 



(jd^ 



and 



mi 



^3 



§ 17. Hitherto we have considered the double BernouUian function as defined by 
one of two difference equations, each of Avhich involves the simple BernouUian 
function. 

We proceed now to prove that, to a linear function of a, oS,,(a) is the only rational 
integral algebraic function of a satisfying the difference equation 

In the first place it is at once evident that ^^^^(a) does satisfy this equation. 
Again the difference of any two solutions is a solution of 

f{a. + cui + %) -— j^(a ~\~ o)|) — f{ci + 0J3) + f{c^) = 0. 
Putting /(a + cui) -- f{a) = (^){a), 

we have </> (a + w^) ■— dy (a) = 0» 

Hence, if /(a) is the difference of two algebraic solutions of the original equation, 
^(a) will be an algebraic simply periodic function, and therefore a constant. 
And thus we shall have 

f{ci + ^1) "^ /(^O ^ constant, 
so that if f(ct) is to be an algebraic polynomial, it must be of the form 

where X and ju are constant with respect to a. 
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Thus the difference of any two rational integral algebraic solutions of 

f{a + o>i + ^2) -/(^^ + ^1) - /(^ + %) +/(«) = ^' 
is of the form \a + /^> Whence the theorem in question. 

[Dr. E. W. HoBSON has kindly pointed out to me that the analysis of the 
preceding paragraphs would be much simplified by starting from the direct definition 
of the double BernouUian function in § 17. 

We should thus define the double BernouUian function by the expansion 



1 — e"^^ a 



So (a) + . . . + (-)"-V'^^'^^-^^^-^ + . . 



(1 — e-^i^) (1— c-^2^') «i«o^ '^ uv^/ « I \ / ^^! 

On differentiating with respect to a, we get the expression of § 15. 
From the relation 

— — — _ ^ 



(1— e"~^i^)(l-^~''^«^) (1— g-'^i^) (1— 6~^3^') 1^0-^^^ 

we obtain the fundamental difierence relations for the double BernouUian function. 
The result of § 10 follows from the identity 

and that of § 14 from 

t t\i - e-^"" V -^ i > \ = pq - ^ -Z^' L^- 

Inasmuch as theoretically the properties of an algebraical polynomial should not 
be derived from consideration of the coefficients of an infinite series, the original 
investigation has been retained. I had already proposed to myself to work out the 
theory of multiple BenouUian functions by a method closely allied to that suggested 
by Dr. Hobson. — Note added Jidy 3, 1900.] 



The Double Gamma Function V^Jya\o)^,(xi^ and its Elementary Properties, 

§ 18. In the elementary consideration of the simple gamma function it was found 
to be necessary to rely on two algebraical limit theorems : — 

(1) Euler's theorem Li^[l + f + • • • "1 l<^g"'^] = 7- 



(2) Stirling's theorem lit 



n = CO 

n ! 



n 



n = CO 



^n+^^-n 



\ 



/2ir. 
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In an analogous treatment of the double gamma function we may expect that 
similar limit theorems will be required. This, in fact, is the case ; but for our 
present purpose it is sufficient to take particular cases of the asymptotic expansion 
for log r(^). 

To make use of this approximation we need only remember that (''Theory of the 
Gamma Function," § 39) if 2; and w be any finite complex quantities, and n a positive 
integer. 



n 



log n {z + mo)) =z log r^ [z + (n + l)a) ] co] — log V^ {z \ o)). 



m = 



We suppose that such values of the logarithms are chosen that additive terms 
involving 2776 do not enter. In other words, we shall say that the logarithms have 
their absolute values, the formula just written being merely a convenient way of 
writing the identity 

(0 = ■'•IW 

On differentiating this identity with respect to z we have 



n 



w 



t - — : = ^^^^\z -\- in ■\- l)(i}\oy\ -- \\sS'^Xz 1 o)) 



with the notation of § 2 of the '' Theory of the Gamma Function." 

§ 19. The double gamma function of z with parameters w^ and o)^ we write 

r^ {z 1 0)^, 0).^), 

When there is no doubt as to their presence the parameters are omitted. From 
this function we form the subsidiary system 

^3^^^ (^ 1 ^1, ^2) = J^ log V,X^ 1 co^, co^), 

'P2^^ (^ I ^i. ^2) = ^o log r^{z 1 0)^, CO.3). 

and so on. 
As a definition we assume 



CO 00 



xIj.^^\z I coi, CO.) = — 2 S t j-~-- — -, r. . 

m,=o m,=o(^ + '?%««i -f m^oy^y 

This double series is, by Eisenstein's Theorem, "^ absolutely convergent, provided 
the ratio of co^ to co^ is not real and negative. This limitation on the parameters 
holds throughout the whole theory of the double gamma functions. It corresponds 
to the limitation in Weierstrass' theory of elliptic functions that V must not be a 
real quantity. 

* V, Forsyth, '* Theory of Functions/^ § 56. 
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We shall show that by successive integrations we may determine T^iz \ w^^, co^) as a 
function symmetrical in co^ and co^ such that 

where /3i(co) = ^/["IttIo)) {" Theory of the Gamma Function/' § 31), and m and m' are 
integers (unity or zero), to be determined in accordance with the detailed theory 
which we proceed to give. 

And the function so determined will be unique, provided 

§ 20. We readily see that the function 



CO 00 



where Q = nij^o)^ + ^^2<^2 satisfies the two difference equations 

where here, as always, we suppress the parameters of the functions ^<i'^\^) 
(r = 1, 2, . . . ) and T2{z) when these parameters are supposed to exist in perfectly 
general form. 

For we have at once from the definition-series 



^^(3)(, + ^^) _ ^^(S)(,) = _ 2 



CO 



-t 



1 



n 



(2 + m.M^y 



(^^ Gamma Function," § 2) 



Next, we may show that the function 
satisfies the two difierence equations 



CO 00 



z' 






1 



1 



(z + ny n^ _ 



whatever be the value of the constant y^ii^i^ ^2)- 

For the series for (f>2^^\z) is absolutely convergent so long as 

1 ^ J^ 
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is regarded as one term, and we may subtract two absolutely convergent series by a 
term-by-term process. Hence we have immediately 



^2^^^ (^ H" ^l) "^ ^2''^^ (^) 



(2) /^\ ™. 



CO 



1 






m., ^ 



and therefore ('' Theory of the Gamma Function/' § 2) 



Similarly, 



t/./> {Z + CU,) - ^,<^) (.) 



K' ( 
V'^ ( 



^2/' 
^l)- 



§21, It may now be shown that the function 



-X) CO 



satisfies the two difference relations 



z A- ft 



1 _l" 



for certain values of the numbers m and m\ provided 



27/1 



TTt 



COc 



2m' 



ITL 



ft)i 



731 («1> o^i) 



Lt 



n n 



1 

03 



1 






loP^ n 



log C0| -- log ft) 



2 



the principal values of the logarithms being taken, 
AVe may write xfj^ ^^\z) in the form 



Lt 



■^ mi = W2-O [■^ "T ^-^^ ^^ ^-^ 



n 






n 



1 



it tt 



t t 



, Ci>i 



and now we obtain at once 
^/) {z + cO - ^P \z) 

= — ^1721 (^n ^2) + L?^ _ ^ .. ^ ^ ,. 

Hence we may take 

provided 

^i72i(^i. ^2) ™ ^/^H^h2) — 2 



'/?z-7r6 



COc 



mTTi 



+ L^ 









51 «. 



' '^^ 1 



"^.'co I ,;,,, ^' ~ -t- 7%ft)3 ,n, - '^ + ^'^^gcwg + {n + 1) &>i ^^ ^ ??.3 - ft^J 
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or, utilising § 1 9 Corollary, provided 



(^irA^h^ ^s) + 2 



mrri 



Lt 



(Oc 



3 11 = CO 



i|;^(l) 2; -j- (n, -j~ 1) oj 



9, 



0)o 



^~ ^]^^'^ [^' '^~ {'^'^ ^^^ ^1 



OJ 



3„i 



^1^^'^['^ + (^^ + ^ ) • (^h "^^ ^2) 



U ';l 



0). 



2. 



OJ 






But (''Gamma Fimctiou," Part IV.) we know that when \z\ is very large, and 



ft) 



not real and negative. 



log r^ (2; + a 1 0)) == 



.t; 4- a 



ft) 



2 



1CV 
og - - 

ft) 



— I"- l0£ 






CO 



- + I log 

ft) "ft) 



+ terms which vanish when | z \ becomes infinite. 

In every case the principal value of the logarithm is to be taken, ie., that value 
whose amplitude lies between — tt and tt. 



Now 



log ^ ^ log ^ z= log 2 
ft) 






in all cases except when z lies in the region formed by lines from 
points —0) and ■— 1 (shaded in the figure). 



3 origin to the 




(&} 



When z does lie within this region, we readily see that 

log ^^ + log ft) =:= log z -]- 2m 
ft) 

if I (ft)), the imaginary part of w, is positive, and 



r4 



^ ft) ''^ 



if I (ft)), the imaginary part of cw, is negative. 



i nus 



ft) 



i) U^g^ + 2/{;7rt} 



log Vi{z '•\' a\m) '=^ 

+ terms which vanish when | z \ is infinite, 



+ I log 



ft) 



where A = 0, unless 2: lie within the region between the axes to — 1 and 
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and where k =^ :£; 1, the upper or lower sign being taken accordijig as 1(a)) is positive 



or negative, wdien z does lie within this regioiL 



On differentiating, we have the derived expansion 



7V 



" at 

+ terms which vanish when'U| becomes hilinite, 

the ])rinci])al value of the logarithm being again taken, and k being determined as 
3fore. 
Substituting in the expression for y^^ii^i, w^), we have 



IJlTTl 



^iTsi i^i^ ^-2) + 2 7; ~ L/^ 

Wo n = Gc 



1 , 1 , 

- - log Ud)^ + - log noj^ 



*«« .V 



ft) 



log ?l . (oj,^ + oJo) — ^ 2:v -^-j +" 



V/i, ■-:-. 7/i., --- 



1 7 n 1 «^i does ] (o), + con) does 

where /v == 0, unless / ^ , x , , Vor ^ S ^ 

(o), + ojo) does not J a)) does not 



lie in the region bounded bylines from the origin to — t and — o^.^. [It is understood, 
of coiu^se, that the principal values of the logarithms are to be taken. 

CO, does 



When, as in the figures 

(ojj + 0)3) does not 



lie within the region of exception, k 



== i 1, the upper or lower sign being taken according as 1(^0) is positive or nega- 
tive. 



Cfi)^+fo)^ 





From the diagrams, we see at once that it is impossible that ^ 1 ' -^ 

0)^ should not 

lie within the region bounded by the lines from the origin to —1 and (/jo. 

Take now m =. k, that is to say, let m be such that we have m :== 0, unless 

' '' > lie in the region of exception, and m ™ i 1 according as I(a)o) 

(oij + ^2) does not J 

is positive or negative, when this exceptional circumstance takes place. 

^' According to M. Poincare, we niciy not in general difii3rentiate an asymptotic expression. The one 
in question, however, may be readily established by the methods employed for log Ti {z + a \ w). 
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Then 



1 n n 1 

— [log no)^ + log ncDi — log n{o)j + ^^g)} "~ ^ ^ 77o 



But this expression is symmetrical in wj and o)^ ; and we must therefore have the 
analogous relation 



^2")(2 + u>,) ~ t/v"e^) = - ^i''\z\w,) + 



Ml 



0) 



does 



where m ::= 0, unless ^ ' V lie within the region bounded by the axes 

(o)^ + wg) does not J 

from the origin to ■— w^ and —1, in which case m =: ± 1, the upper or lower sign 

being taken as I {o){) is positive or negative. 

Provided therefore that 



1 '' '"' 1 1 

logti --tt^ --, + "— {logwj; + loga)^ — log(a)j + oj,} 



w^e have, with the assigned values of m and m\ the two difference relations 

The function ygu (^1? ^^s) ^^^ propose to call the first double gamma modular form. 
It will subsequently be expressed in terms of the function D(t) introduced into the 
theory of the functions G {z\t) {'' Grenesis of the Double Gamma Function/' § 4). 

It will be seen later that the algebra of the double gamma function would have 
been slightly simplified had a modified value been taken for this function 72i(*^i? ^2)^ 
and the analogue shortly to be considered, y^-ii^i, 0)2)' ^ ^^^^ not observe this fact 
until the theory had been completely developed, and the matter is scarcely of suffi- 
cient importance to demand the labour which such a change would entail. 

Covollm'y, — Notice that it has been proved incidentally that 

n n 1 

% ^' ^^ 

is infinite, when n is infinite, like log n, 

§ 22. As the numbers in m and 'jn! enter constantly into the analysis, it is necessary 
to consider their properties. 

Suppose that the functions log z^ log ^z^ log ,^z are natural logarithms (with e as 
base), which are real when z is real and positive, and which are rendered uniform by 
cross-cuts along the axes of —1, — ojj^ and — w^ respectively. 

2 p 2 
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hen it is readily seen, that 

log,3, (^i + (^h) - %o.,^^i 



log (o)^ + oj.^) 

log («| + ^1)2) 
log (ft)^ -f 0)^) 



log 0)1 - 


zmm 


log «3 


2??iVi 


log 0)3 - 


■ 2m'TTi 


log OJj 


2'rmri. 



^ log ((1)1 + Ct)Q) - 

By inspection of a diagram we see that m and m both vanish if the difference ot 
the amplitudes of m^ and a)o is less than ir, these amplitudes being measured between 
and ± IT positively or negatively from the positive half of the real axis. In 
particular -when the real parts of co^ and mo^ are both positive, m and m' both vanish. 
Not only so, but in all cases either m or vi^ must vanish. 

Again, if the difference of the amplitudes of m^ and o):^ is greater than it, m and m 
cannot l)oth vanish, in fact^ in this case we have the important relation 



ra 



rri 



riz ^ 



the upper or lower sign being taken accordhig as I( ') is negative or positive. This 
result is intuitive geometrical!}^ ; in the figure, for instance, two cases are indicated 



in which ].( -^ ) is negative. 



fi) 



1/ 






&i: 






/, 



/-^' 

m^ 




For corresyjondinff to the unaccented value of o)^. 



m. 



m 







and corresponding to the accented value of cwg^ 



m 



■m 







1 



Thus in both cases m— m ™ L 

No such simple expression can be given for m + m^ a number which is of constant 
occurrence in the higher theory. 

However, from the values for m and rri previously given, we see that when the 
axes of a)| and &)o include the axis of — 1 within an angle less than two right angles, 
the values of {in + qu) are given by the table 
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('/1I + m-) 


I(03X) 


I((Dl + 0)1 ) 


1 


+ m 


+ ve 


^1 


+ ve 


~ ve 


_^1 


- ve 


+ m 


1 


- ve 


- ve 



and therefore w + ^^^' = i 1? the upper or lower signs being taken according m 
l((jt)i + w^) and I(o>|) have the same or opposite signs. 

§ 23. It may now be shown that, if C have any arbitrary vahie, the function 



",•% 



Fg \z\mi, %) =z C e- 



.^ 'L. y2i(oji, wg) + jyaaCwi, Wg) 



00 CO 



. z X n n' 






L\ 






where ft = m^coi + '^^^ooj^, will satisfy the two diiference relations 






AC'- 1 <»3) 



V 



27r 

Cl)0 



e' 



•2'r/itrt 



(s _ 1 I » 



ztt 



V 



2/ 3 



ft)l 



where m and m' are the numbers previously specified, and 



n = m 



^ J;, 1 6>T + ct>Q 1 

^^ ^tOj^tOo 



1 






lOff 1 + 



n 4- 1 



tt>i 



©o 



fc5 



Ci>i 



Ct>i 



log 1 + - 



Wo 



+ "i—-"'^^ l\og{0Jy-\- W3) 



Jct)jCt>c) 



log 6)| »— log OJg} 



the principal values of the logarithms being taken. 

Observe that with the notation previously introduced {"Theory of the Gamma 
Function," §§16 and 31) we may write these difference relations in the form 

The i)roof, to which we now proceed, is exactly analogous to the one just given. 
We have 



, o Hz + mo) 



Tilcaj2 ^ y.2\^(ai + y-jgWi 



e-^ 



9. 






■;?- = OD 






4- (Wli + l)ft)j 4- MgCWjj 

2za)i "^ ^1^ 



X e 






20^ 



where y^i and y^o are understood, as always, to mean joj^i, ^0) and y^^io^j:^ ^2). 
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Substitute now the value of 721(^1, (o^) obtained in § 21, and we find 






\jJ%. O. 



iuj (jt/ T (a) :> 

A -I- <^j 



^ ~h (jt)j 



X L/; 



■/?- = 00 



exj 



3. S :: 



2( 



log ^i 



1 n 11^ [^ + ( m, + l ) + m ,a),^„,g] 



exp. 



^^^^^ Y^ ^ {log a)| + log 0)2 ™- log (6)2 + ^>,2)} + W|yc 






23 



, r^ (z coo) 



/\ .1, jc- 



-/l = QD 



i^-X.. L), 



i-J/i/ 






« '/I 



og?^ 



,6);, 



i^ 1 



J:\[,^ + Oi. + l).(a), + %)h3] 



1 



~;^ + {11 ~h 1) 6>i I Wq] I'lf;^' + (n + l)ft>3 I Wg] 



by the employment of the identity of § 18, Their principal values must throughout 
be assigned to the logarithms. 

But, as has been seen in §21, from the formula obtained in the '' Theory of the 
Gamma Function/' §41, we have when \z\ is large and z not in the vicinity of the 
axis of — cti, 

'z H- a 



logV^{z + a\a)): 



(i) 



\ ■ z 

2/ {^^S ^" + 2/{;m} — -+ log pi(6)) 
+ terms which vanish when % \ is infinite. 



where ^ = 0, unless z lies within the region between the axes of —1 and -co, in 
which case l = -^ 1, the upper or lower sign being taken as I{oi) is positive or 
negative. The principal value of log z is to be taken, and the prescription to be 
given to log V^{;^z + a| co) is left indeterminate : it is obvious that we only get additive 
terms involving Itti, which vanish in the sequel 
Inasmuch as when n is large, none of the points 



% 



+ {n + l)(wj + 0}%), ^- + {n + 1) wj, and z + {n + 1) oj^ 



lie in the vicinity of the negative direction of the axis of wg, we may substitute the 



values given by the asymptotic expansion in_ tlie expression for 
We shall find 

r,- (.) r, {z 1 0),) ^ 



rrX^) 



exp. L^ 



% = CXi 



^^^^ {% ^^^2+ log 0), - log((yiH" c^g)} + ^iJn 



ft);) 



z 4- C^?' + 1) coo 



4- 1 loff r^oi 



^ 



o 



'^ 4- (n[+ 1) ft)^ \ 7i(6o + a»o) 

— 2" /log ^'^Ct)] + ^^ '"'' 



ft)o 



2 f^ "^ i^L±-l)^(^L±_„^^^^^ 

\ «3 



Wo 

CD.i 



\ 



i h'^'' 
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In this expression their principal values are throughout to be assigned to the 
losrarithms, and the numbers h are to be such that 

k^ ■=! 0, unless (oj^ 4~ ^i) has within the region bounded by axes to -— 1 and — wo^ in 

which case 
h^ =z :i: 1, the upper or lower sign being taken as I(ft>o) is positive or negative, 

while 
k,-^ ■= 0, unless Wj lies within the region bounded by axes to -- 1 and •— co.,, in which 

case 
Jco^ = dz 1, the upper or lower sign being taken as l(ft)o) is positive or negative. 

On reduction we now see that 



r.-'(z)i\(z\co:,) 



exp. Jjt 



(^>yisi<^i, <^i) + tt' "^ + {n + 1) "' - '^-'log«{coj +Wo) 

"^-^ ("^i 

(2n -j~ 1) 6>9 — 6>i . (2'ib -f- 1) (Oi — (jOo 1 1 / \ 

^ - — ~ log nw^ — — — ^ ' log tift)! — log p](wo) 









z 4- i^n + l)(fW| 4- tt)^) 



CO 3 



"o 1 Zi/C^TTt 



z 4- 0^ + .l)ri>x 



(Wo 



I ) 2korrL 



We must consider the three possible cases in which /;, and k^ do not both vanish. 

-I 

(1) Firstly, when i/' / x , \ li^ within the region bounded by the 

\.\(^i ~f~ 0).)) does uot J 

axes to — • 1 and — o)^* 
In this case i ,^ \ the upper or lower sign, being taken as I(wo) is positive 

L A/o '^^^^ -f- 1 J 

or negative. 
And we have 



r 



r. 



3 ^{Z + ft))) / ^'^ T , 



CO 



V90 



/6 /4 1 



2 2' „- + (« + i) ""' ^^''log »(«i + wn) 

^^ COjCOo 



log nw^ ^;7^ - - log ncoi 



2<:o,&)o 

{V, + 1) TTt 1 -^ ,-, / S 



^ft)o 



TTt 



the upper or lower sign being taken as 1(0)2) ^^ positive or negative. 
But in this case m =: ± 1? the signs being chosen in the same way. 
If then we take 

73^(^1. ^^s) = S S^ ^ - {n + I) ~^- - log ti(o), + CO,) + '- -^:: :7 log noj. 







II 



ft), ft) .9 



2&)ift}o 



, (2^2, + 1) ft>, — &)o -, , ,^ (n 4" 1) TT^. 

2ft>,, ox ^^ 
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we shall have 



(2) Secondly, when \ / , . . [ lie within the region bounded by axes to — 1 

L (wj + wg) does I ^ ■^ 



(1 •— 
















111 this 


case 


Ay .7 


- + 


1. J 
1, 


m ^ 


0. 


We 


shall h^ 


bve 


then 


the 


same 


relation 












p^ 


rX^ + 


•^i) 



provided we take 






•n n 



/ t / L 1\^ "+1^3 1 / I \ ! (^''^' + <^;3 "~ ^^11 



r33(«i> wg) =: S ^ J, - {n + i)— ;7;-- logn{o}y + co^) + ' ■'^--.- ~^log nw^ 

(2'ft + 1) O), — 0)3 , _^ „0l + 1) TTt 

the upper or lower sign being taken according as I(ft),) is positive or negative. 

(3) The third case, when \ .^ , . ^ I lie within the refifion bounded by axes 

[{ojj + w^) does J "^ 

— 1 and — oio, is easily seen to be impossible. 

In all other cases we shall have the relation (1), provided 



729.(^15 ^2/ ^ S S ~ — (fh + t) -- - log W (oil + CO2) + ;. 1 



'^■^ 0)tO},^ /i 



ju/^ *-/ft)iC!>o 



log na).2 



(2/1 + 1) o)i — . 

+ ^; log ^^^^^1 ' ' ^ (2). 



Suppose now that we had investigated similarly the quotient 
we should have obtained the difference equation 

where y^i (^n ^2) ^^^^^ ''^^^^ value D, let us say, given by equation (2), except in two 

cases. 

f /y does 1 

(1)^ When i , ^ , . . ^ i lie within the region bounded by the axes to — 1 

[ [(0^+ (Wgjdoes not J 

and •— • o)i^ in which case 

(n + 1) Tit 



722 (^1? %) = D i 2 



®i 



the upper or lower sign being taken according as I(cO]l) i^ positive or negative. 
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r 1 

(2)' When ^ '^ . ^ [^ lie within the region bounded by the axes to 

a)|5 in which case 



1 and 



72-3 (^1. Wg) rr= D =f 2 



{n + V)iTi 



COo 



But the cases (1) and (2)^ are precisely the same, as an inspection of the figure 
shows at once, except that [(ri)^) is positive when l(wi) is negative, and vice versa. 




And, similarly, the cases (2) and (1)^ are the same, with a similar change. 
Hence the values which must be assigned to 722(^1? ^2)? ^^^ order that the equa- 
tions 



/- 



Fg 1 {z -f 0)^) r^ (z 1 0)3) 



J 



L r,-^ (.) 



^/27r/a)3 



e 



e 



2miTi 



■2m'Tn 






may co -exist, are precisely the same. 

We shall have then these two equations, provided. 

722 (^i? ^2) =^ ^^ D, 



n — 00 



where D stands for 



n n 1 



tt' 





{n + i) 



1\ «! + W 3 



«1«2 



log [W (Wj + ft)^)] 



+ (2!L±J>iZL»i log ^,^ + (2_'i+ 1) -1 



9 



^1^3 



2a)2C()^ 



(Da 1 

— log ?^a>|^ 



(the principal values of the logarithms being taken), except in two cases, 

(1) When {»' ^r . , , I 1- «*- «« -g- bounded by the axe, to 

^ ' 1(0)1 4- wg) does not J ° •^ 

and — Wn, in which case 



1 



722 («i. ^-2) = Li 



D +2 



m{n + l)7r^ 



a)o 



0)2 does 



(2) When ^ 7^ """"" \ t Mi^ within the refijion bounded by the axes to - 

' (o)^ 4~ ^2) does not ^ 



1 



and — co^, in which case 



722(^1. ^2) = L^ 



n ~ 10 



D + 2 



')vJ(n + l)7rt 



0) 
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Since m or m always vanishes, we see, on combining these results, that in all cases 



73-2(^1. ^2) — ^i 



But 



n = CO 



n 



^t -~ {n -f- I) --7--— log n{a), + w,) 

'^^ 6i)|ft)Q 



(2^/1 4 1) coo — Wi 



2&)iCt)o 



1 , (271 -I- 1) O), 

log nm 4~ - ^y^---'- 



f».> 



r/l. . Qll^^ 



■^ log' 'Hoij + 2 - + — ) (n + 1) 



TTt 



-^ 



'y? -: CO 



0^^ 



CO1 4 cOo . 

■ ,^ " log n 

1 



(ti +1) ^^^ (lc>g{^i + ^^%} "^ l(>g ^1 "~" 2m7rt 



CO: 



1 



/ 



(n + 1) -"-^ {log (w]^ + C03) — log c^3 — 2m Vi.} 



^ r. I log (oj{ + ^-2/ ""~ ^Og c?ij 






OgWg} 



/ ##1 

log io)^ + COg) — log OJj — 2?}l'7r6 := log 1 + "~^ 

log (o)i + ^2) -^ log o}^_^ —- 2m in = log f 1 4- ^ | 



the principal values of the logarithms being always taken. Hence 



722(^1? C02) == -Lt 



n = 00 



■■ n n 1 
0^^ 






log ^l 



^^ 4- 1 T / _. . oj.> 



ft)c 



iog I 1 -p 



71 + 1 T / 



CO, 



log ( 1 + - ^ 

COj + CO3 .-. , V 

-f" -— — — — v^og (^o>| ~|-~ cogj '"- 



J- .<j 



log Wy 



As a corollary, we see that, when n is very large. 



n 



n 






m^ = r% = "i^liOy^ ~f '^^%^3 



r « ^r^'j 1*1 COi-rCOo 

- IS innnite like -i - - io.o" n 



iCOiCOo 



o 



'?^ 



COr 



loF 



to 



log COg} 



1 "-f" — 1 

1 lop' ( ] 4- ^ ) 



24. We now determhie the constant C in the expression 



Ce 



^ 721 + -yn 



00 



'7 



CO 

r-r/ 



■J/li = 55*0 = 



1 **-p -pr ) e> 






by the condition assigned in § 19, that 

l^t zV^ ( 

g = 

This at once gives us C = 1. 

It is evident that, with the conditions of § 19, one and only one fimction can be 
constructed, and this is the double gamma function V^{z\m^^ c^g), which is such that 



2^ I 0)|5 COo) 
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r.3 ^{z\a)^^,a).^) = e'' 



'.) Vs] + -yas 



CO 



z 



. n n' 1 + 






a: 



6 






where y^2i{^i^ ^2) ^^^^^ y^J,o)i, oj,;^ are two constants, which we call the first and second 
double gamma modular functions of the parameters o)^ and cocy. 
These constants are given by the relations 



731 (^1^ %) ^ L/; 



n = CO 



1 



ft)l (Do 



log n 



^ 






m, = Jii., — 



4 J/ 



1 



+ ;;7, (log ^^i + log ^^3 "- log («i + 'ys)) 



72.3 (^1. %) = L^ 



5i = 00 



71 4- 1 



0): 



n 









- log ?l •' log ( 1 + 

At 



i(jOi(i)rt 



G>i, 



SD 



a> 



1/ 



log 



1 + ^f J j 4. ~|^-^^ {log (oj, + 0)3) ^ log m^ - log oij 



where the logarithms are such that their principal values must always be taken. 

And the theory is the natural extension of that of the simple gamma function 
r^(^z\a)i), which is such that 



CO 



T{^^{z\o)^) — e~y^\ z, n 



m^ = 1 



\ 



l_l_ 



Z 



mifDi^ 



z 



e 



riiiMi 



J 



where the constant yn is given by the relation 



TriK) =^ L^ 



n =: 00 



n 



1 






CO , 



logno)^ 



and the principal value of the logarithm must again be taken, 

I 25. We may now see at once that 









mm 



— imri 



For we have seen that 



where m is either zero or unity according to the determination of 
But Jjt [zT^{z)] z=z 1, 

and U[zV^{z\(^,)'\=l, 



Zi Zii 



2 = 



as we see immediately from the product expression for 

Ti {z I o>g). 



2 Q 2 
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Hence, making i^ — 0, we have 



rii'iTi 



whAch is one of the relations required. 

We thus see that T^{o)j) is independent of o)j, and not only so, but its value is 
substantially a quantity that appeared several times in the theory of the simple 
<^amma function. 

Thus we saw [§ 4 cor. " Gamma Function "] that 



that (§ 8) 



n — 1 

II r 



"CO 

logF 



n 



OJ r-r n 



2^/co)]^^^^J 



i:^ 



CO 



)dz 



a) log ^(27T/a)) 



ind the most ^s:eneral form of Stirling's theorem, was seen to be 



iD 



log ri [a "f m,oj) ■=: p[n\ogii ■— n] + ?^{S/^'^(a + o})j)o}logp(.o} 

viy := 

+ ll + 8i{a)] logn — logi\(a) + log y(2'/r/a)) 



+ S/(a + co)logpa) + S 



I. (-)"'"' Sm(a 4- o)) + ,Ii, 



4- 1 



■»i — 1 



'??171 



"iii 



{]m) 



m 



We now see an additional reason why it was proposed to write 

and to call yii(oj) and pi(oj) the two simple gamma modular forms, the latter being 
sometimes called the simple Stirling modular form, We shall see that there exist 
three double gamma modular forms 

of exactly analogous nature. 

§26. We proceed now to connect the function r.2i^z\a)^^ oi^) with Alexeiewsky's 
function G{z\r), some of whose properties were investigated in '' The Genesis of the 
Double Gamma Functions." 

In the first place, we take r = <^2/^i> ^^^^ then we have 



G 



z 



CO 



z z 

a — +b - — g 

(3 W2 l(0./ 



z 



00 



OD 



— . n n' 

ft>2 mi = m^ ~ 






1 + ^ )e «'^ s iP 



where U =: mj_co^ + "^^h^^, 



and wherein 



a " - los: 2ttt + i los: ^° 
z ^ ' ^ ^ 



C(t), 



6 



loi 



9 o 



T 



6 



t^D(t) 
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We also have 



Ff^z) -e^^"^ 






00 CO 



' . z . n n 

Trli ~ 1112 »= 






so that on comparison of the two products we find 



Tf^{2) = Ke^^''^'G' 



z 



\«i 



where a, ^, and K are suitable functions of Wj and coo. 
Now G{z I t) sa.tisfies the difterence equation 



z 



G(2|r) = r-i(~ )G(.?+ 1|t), 



and hence G 



;<; 



00 



T I satisfies the relation 



i-._ 



f{z + 0)i) =: 0)3 -. Ti (^ I C03)/(2^) . 



Hence a solution of the difference equation 



IS 



f{z + 0) ^) = ^(0)^/2^) . Ti (;^ I o)^)f{z) er^"^^^^ " "^> 

s'^ 2 g _«_ / ^ \ 



And it is evident that the coefficient of 2mTn in the last exponential may be 
written ,So(. I a>,,a,,). 

The general solution of the difference equation is 

where |9(^^| coj is a function of 2; simply periodic of period Wp 



J — LCxXOt? 



Jiut 



^2 (^) 

' T(0) 



P (^ 1 ^i)- 






has been seen to be an expression of the form Ke*"^'"^^^^, where K, a, and 



/3 are independent of z. We thus have 



so that a = and /3 



2r 



iri 



(i)i 



^a (s + Mif + ^ (.t; 4- (0]) __ y3a%^ + ^z 

where r is some integer. 



Hence 



and since 



2rirrts 



r^-X^) = Ke "a (w3e-'""'')^^«^^>(27r)S^G[-| T 



• # » » 



. (IX 



ijt izT^(zy 



3=0 



jjt 

K=0 








g(^ 


')' 



0) 



2> 



we have at once 



K 



(O^. 
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Now G 






0)1 



T j satisfies the equation^^ 



1—1 



r(-)/(^-) 



T to 



l^ ^^' ^ (27r) 2 



where r "^i^^ ™e ^"'■^i ^; ^^^^ ^^^^^ j-].^^ principal value of log r is to be taken. 
(The same remark, of course, applies to every many-valued expression of this nature 
which occurs in the course of the investigation.) 

Employing the relation (1) in conjunction with this equation and the equation 



O 3 



we obtain 






M'-' 



^^"^i^ "^ ^3) ™- A (^1^1) 



2(0 



+ ^ 



1 e 



{2rr) 



T 



which reduces to 



y)2{m~-rii') -m 8/(2 | ojj) _^ .;, 'Irmr + S'i(a I to,) [log oj,, — log W| -■ log r] 



But we have seen (§ 22) that 



log 



0) 



2 



loi 



CO 



1 -"-^s 



log T =^ 2 (/n — m') TTi 



for 771 — m^ = 0, unless the difference of the amplitudes of o)^ and ft)^ is greater than tt, 
in which case m — m^ = ih 1? according as — I M is positive or negative. 

We thus find r ™ 0, and incidentally we obtain a valuable verification of our 



results. 

And now finally 



Vf^ {z) ~ 0}.^ (27r) ^-^ (oi^e"-'^'^"^ Y^'"^ ( 



o; 



OJi 



T 



the relation between the two forms of double gamma functions. 

§ 27. From the relation just found we may at once express the gamma modular 
constants C(t) and D(t) of the former theory in term_s of 722(^1' ^-2) ^^^^ 72i(^i? ^2) 
respectively. 

For we have 



G 



\«i 



00 00 



(f-,^^%^ . -^~ . n 11' 



0>0 '/ilT=0 )/lo=0 



o3 n 






1 + -■ c « * '^"' 



where 



a 



T 



- log 27rT + ^ log r -\-' yr — 0{j)^ 

27rT^ 



6 = — • T log T 



6 



"■^D(t), 



and also 



CO CO 



Tf\z) = 62^='+'^^^ . z . n n 



^7ijrrO 7%=0 



2; 2;^ 



2; \ _.- + 
\ !_ — - J ^> O 203 



■^ '' Genesis of the Double Gamma Functions/' § 10. 
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Substituting in the relation 



r2-^(z) = (27r)~2«'.W3 . (wjje--""^')'-''"^ . G 



M 



rU 



we 



find 



77721 («i» f'^a) + 2=722(^1, w.) 












^ 



2ft> 



log 27r + ( --'^--~- — f- - \, + 1 lo 



g^2 



2m'7ri 



^2 







9 

/^.V 



2ft>^ft>2 2ft>| 






ft>c 



O 2 ' 



And hence equating coefficients of z and 2;^ we find 

721 (^1' ^2^ __ 



2 



1 & 

-— {log % — 2m7nj 4 



722 (^1. ^2) ^ 



2caTft>o 

J. ^ 

1 



9 2 ? 



'2 



1 



{log 27roj2 — 2w2m} — „^" {log co^ + 2m7rt] 4' 



a 

0)g 



Thus 



# 



D (t) = - oj.^rri (^^1. ^2) + "f (log C02 - log r} - 






6 



ft) 



3 



C(t) 



CO 



iJn {^i^ ^2) -{h + 2;) {% ^2 - log T --" 2m7n) + 7^ 



which are the relations required. 



Since 



log 0).^ — log T — 2mTn = log w^ — 2m^m, 
we may evidently write these relations in the form 

D(t) == - (o^jYu {(o^, 0)^^) + ~^' {log oji - 2)11 wl] 



'TTo 



c* ? 



C (r) - - 0)1722 (coi, C0.3) ^ / 1 + ^~-ij {log (o^ -' 2m Vt} + y. 
§28. We will now show that, when the parameters oj^ and o)^ are equal, we have 



721 (^. ^) 



1 



60' 



log CO — 1 — y 






and 



722 (^. ^) = ^ [7 



1 

2 



log CO^ 



By the definition formula of y^iicoj, o)^), we at once have 

1 



^^21 (^. ^) 



L?5 

-11=00 



n n 



m, = Om,=o{m,i + W2) 



2 - l^g 



tlG) 



Group together all terms for which m^ + m^ = e^ and we have 



n n 



« e + 1 , 



7^ 



+ ;^"4+.-. + 



^ ^' Genesis of the Double Gamma Functions/' § 6. 
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for we may suppose that the terms are represented by the corners of the small 
squares mto which the positive quadrant is divided — -in the new grouping we take 
together all terms lying on a line equally inclined to the two axes. 













' 












\ 


..... 


... — _ 


-- 


V 

\ 


> 


\ 








\ 


L 







rrs 



rhus 



n 71 



(^?h + ^^^•'>) 



So that when n is very large 



n n 



1 



11 - 



1 



3 









^ (n + If ^ {n + 2)' 






1 



- .-4- 4- 



(71 -f- Qiy 

n — Qi 
(n + ny ° 



(w2i H- mo)^ 



^ ' ^ * 6 Jo(l + ^^ 



+ terms which vanish when n becomes infinite 



log ^i -|™ y + 7.^ + 



6 



[ ^^ 1q(t 2 + similar terms. 



Hence 



721 (^^. ^) 



1 



ft) 



r 



log Oi 



6 



which is the first relation, 

In the second place we have, when o)^^ and cy^ are equal. 



^722(^5 ^) ==" ^^ 



n n 



1 



and by the same method of reasoning as before 



2n log 2 — log 2no) 



n n 



1 



^^h + ^^^2 



^ e -f- 1 , n 
^^i € n + 1 



. n — 1 



7^ 4- 2 



I 



1 

7^ -4- n 



n I I 

S - + n + " 7 + 






1 



6 a e 



r 



n "f ^2, 



, / 1 1 1 



n 
n 



=: n -{- log 7?. + 7 + log 2 + 2n log 2 — 
neglecting terms which vanish when n is infinite, and thus 



1 

2 



n. 



1 



722 K ^) ==^,t]7 



1 

2 



log Ce)]. 
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It is an interesting piece of work to show that these results accord with those 
previously obtained for the G function of parameter unity. 

§ 29. We proceed now to write down the expansion of log Vo{z) and the first few 
of the derived functions t|;./'^(r^) in the vicinity of 2; = 0. 

We have, by definition, 



CO 00 



» • e 



\Ij:^^\z I o)i, CO..) = — 2 S X r —-^ J where O = m^m^ -f m.>oj,., 

Since the series on the right-hand side is absolutely convergent, we may expand 
in the form 

2 f «>,, 1 ^ , Sz "^ o . 4: z^ 

^3 ^' (^) = "— "i; — 2 -^ S ^ ^o •— S S ^i + S S :j ^ Yy^ . 

Hence, integrating, 

1 ^ Z ^•' f Z'^ ^ t ^^ 

^/'^ (^) = ;r — ysi (^i^ %) ~ ^ 2 S' ^3 + ^ S S ^-^ — 4 tt ^^ + . . . • 

the constant being determined by making 2; = o. 

Thus integrating again, and determining the constant in the same manner, 

1 00 ^ CO 0«1 

,/,3<l> (2) = — -^ — y^i («^, ft)^) Z - 733 (ftJi, W3) - S S' ^3 + S S' ^, - . . . . 



and finally, 



A-.W, 



^A/ 00 ly^ CO ^* 00 2;^ 



log r^ {z) = - log z - zy,, - ^-f -tt' ^^ + tt' ^ - tt' ^ . 
the expansion holding good within a circle of radius just less than the least value of 



^^h ^i + ^-i (^z 



nil = 0, 1, 2, . . . c/D 1 

m:, =: 0, 1, 2, . . . CXD J 



excluded. 



We note that by Eisenbtein's Theorem each coefficient in the series is an absolutely 



convergent series. 



§ 30. We proceed now to the expressions for the double gamma functions as simply 
infinite products of simple gamma functions. 
Consider the product 



m^ — 1 






The typical term may be written 



63! 4! 



00 

and the series t ylfi^''\m.2(i)o\o)i) are absolutely convergent when r>3. The product 



iH.y = 1 



is therefore in general absolutely convergent. 
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Agairij it has no finite zeros, for its zeros would be those of Tj{z + moO).2.\^i) ^^^ 



■111: 



0, 1, ... 00 . And its poles are given by 



z t\- m^o)y + rrkj^ i)^ = 






\J 9t I * ♦ ♦ 

0, 



A ~», (^ fl « A 



CO 



Thus --|v bas no zeros or poles in the finite part of the plane. 



Change now z into z + w^, and we have 






2; 



00 

n 



i -J- 



6 «^2'»>2 



'?7loft)o 



<« <v / 



00 



/S. 11 



^//?2 = 1 L 






;?]i.>0)., . e 



Wl,l//j(l)(m2t02 I Wi) »- JL_t: J^-!,'/|(2)(in2W2 1 Wi) + 



''^'2^3 



Now the product last written must be convergent, for all otlier terms of the 
identity are finite for finite values of \z\^ and this product is evidently of the form 
gi'~ + a^ where j) and q are functions of Wj^ and w^j only; 

Hence Ave must have 



Now we have proved that 



r, {z + 0,,) = r, {z) v{-^ (. I o>,) 



'^TT \-^ ^2hnn 



ft)., 






Thus if we put 

we shall have 

and similarly we shall have 



f{z 4- Wi) 









__™, ^a^.i + ^, 



/■)(^^>^ 1" PJr 



Hence /(a;) is a doubly periodic function of the third kind, with no finite zeros or 
poles. 

Thus we must have 

/'/.A _« /1., A,;-:i X lb 

for in Heemitk's expression of such a function the cr functions are each associated 
with a finite non-congruent zero."^^ 

To determine C we put ;4^ = 0, and obtain 



C = Lt 

z = 



Pi {^ I Wl)_ 



Bifierentlating logarithmically, the identity 

Fo {z^ I CO, cuo) ■==■ P (^) e'^'' ^' 



B; 






Forsyth, 'Theory of Euuctiona/ § 142. 
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we find on making z = 0, 



1 
Jm (^i» ^2) = B + - (log w^ — y), 



«i 



and therefore B = — ■ y^s (^1? ^2) 4" 



G)| (Oi 



Differentiating again, we have 



00 



y\>^^'\z\o,,, 0,,) = 2A + ti»/2)(2lco,) + t [t///^>(^ + ^isco^lc,) - ^^''{m,<»,\oy,)] 



riio ~ 1 



Again, making z = 0, we find 

or A=: — iy^iico^, o)o). 

Finally, then, we have 

T^[z\(i)^, (0.2) = e - ^ ^1 ^^^XT{z\o)y) 



00 



X n 



r. (^' + rn,,0):, I ft)-, ) - #1 (^) (r/l2C02 I COi) - f- Xfji (2) (»?oC02 I COi) 

_.-.™_^ ___„„™, y> ^ 

. I\ (moO). i 0)^) 



This formula is equivalent to the one obtained in the " Genesis of the Double 
Gamma Functions,^^ § 2. 

It is an interesting verification to actually transform the one formula into the 
other, making use of the relations established between 721(^1' ^2)? 722(^1? ^2)' ^K'^) 
and D(t). 

On account of the symmetry of the present functions, the formula corresponding 
to that given in § 8 of the '' Genesis of the Double Gamma Functions " may be 
written 

Tcy {z 1 wi, 0)2) = e ^ * '^a ^=i ^ . Vy{z I W3) 



00 
X II 

mi — 1 



' V^{Z + m^ft)^|ft)o ) - #1 (I) {m,o,y I C02) - 1 1//1 '2) (.wicoi I C02)' 



The product formulae just obtained correspond to the expression of the cr function 
as an infinite product of circular functions. 

Such a circumstance, of course, at once prompts us to try and find a formula 
corresponding to the expression of the a function as a .st^m of exponential functions. 
But it is readily seen that such is an impossibility. We cannot express the double 
gamma function as a sum of an infinite series of simple gamma functions of varying 
arguments. 

It is this fact, combined with the absence ot any quasi-addition theorem for the 
double gamma functions, which precludes the possibility of any collection of formula) 
rivalling in number and elegance those of the doubly periodic functions. 
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§ 31. We pi'oceed now to express Weteiistrass' elliptic functions in terais of double 
ffamma functions. 

In Weierstrass' notation of elliptic functions we have 



00 



F'(2) = -2t t 



I 



Q'lii m.2 \^ 1 ^^) 



8 ' 



00 



where O = ^^i^i + ^^^2^3? ^^^ a)^ and o)^ are the periods of ^(2^). 
Now by definition 



00 a> 



Ml, = Ma = V' T i-i) 

Representing the various terms by the corners of the parallelograms of the figure, 
we readily see that 



^2® (^ I «^i> ^2) + ^^2^^ (2 



<^l> «2) + V'3^^' (2 I Oil, ~ CO,} + l//3«> (2 
: - 2 



"*— COj^p — ft) 

1 



■2) 



00 1 «3 1 



rt 

/^.'I 




Hence^ nsing the natural summation %\{j^^^\z j J^ ft)|, ± Wo) to express the left-hand 

side of this relation, we have 

2 
txfj^^^^ {z \±a)i,± (0.2) — f\z) + t^^f'^z I d: ^i) + tf^'^ {z I ± ft)^) + :^ . 

Ay 

and therefore, on integration, 

F (z) = St/./) (^ 1 ± 0),, ± (.3) - St/^i® (z I ± 0.1) - S^i® (2 1 ± CO,) + i + ^, 

where j^ is constant with respect to z. 
Evidently 



00 1 CO j 



(^zwj^y^ 1 ('^^fw.O^' * 



Now 
721 (^u ^2) == T^'^ 



'vl = 00 



1 



n n' 






[logncuo + log l^ft)l — log n {mj -j- wo)] — S ^ 



1 



r,H - m, - ('^h^l + ^^h^if J ' 



where the principal values of the logarithms are to be taken. 
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And hence 



V 



JUtt 







NoWj as may be readily seen by examining the different possible cases in a 
diagram, 

log(a>i + o)o) + log[— (o)^ + %)] — log (6)1 — %) — log[— (o)^ •— 0)2) 

-, 6)1 + ft)n 

= 2 log — — — '^ 5 



where that value of Wi^^cwo is to be taken which is on the same side of the real 

axis as w^ + %. 
With this proviso. 



p 



^J ^ + :^ log ""^^^ 






ft)i 






the infinities in the double summation being equal in absolute magnitude. 

I 32. We may now express Wriebstrass' ^ function in terms of derivatives of 
simple and double gamma functions. For on int%rating the relation obtained in the 
previous paragraph we find. 



remembering that 






i m 



9{z);^ 



— C(z) = St|/i'''(«| ± wj, ± Wo) —trli^^\z\± oj|) - :S»|'/i'(zi ± Wg) — ~ + vz-\-ix, 

where /x is constant with respect to z. 
Making then ^ = 0, we find 






= - ty^^{± 0)1, ± m^) - ~ {log(ftiO - y} + ;7log{- Oil) - y] 

(Ml Cut 

1 1 

- -- {log 6)3 -y} +- |log ( - 0)3) - y} + /x. 



Now by § 23, 

Sy32(±«l, ± Wo) 



6>i -f- 6>a 

1 j<^ 






1ii » ii fa wMi 



« 



\"^1 



log (cOi + cog) -- log 0)1 - log 0)3 

- - % [~ (^i + %)] + log ( - 0)3) + log ( - 0)o) 

ZL«J J ^^g (^3 - ^1) - % ( - ^i) - log 0)3 

«1«3 1 "- log (0)3 — 0)0) + log 0)^ + log (— 0)2). 

"f 9 ~){l0g(0)i + %) ■— l0g[— (0)i + 0)o)]} 



Z6>i jijCWo 



fi> 



1 



+ --[log( 

ft), ^ ^ 



«^i)-l"g«^i} +- {log(- %) 

ft>o 



loo: 0).^ 1 



I _ 97" H log (0)3 -- 0)^ - log(o)j - 0),)}. 



* JoEDAN, ' Coiirs cFAnalyse/ 2iid edition, p. 347. Note that Jordan uses 2(ui and 2ct>2 instead of 
o>x and Wo of this paper. 
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IJeiice 



1 



1 



/^ 



2« 



"t~ o )[^'^s('"i ~t~ ^''■-i) — ^^9^1 — (^i ~i~ ^3) 



1 



2ft) 



2ft), 



1^ 



• [log («o - Wl) - log (ft>; 



%)> 



(O), + «,)] = ± ... 



Now log (cdj^ + ^2) "^ l<)g 

according as I(a)| -|- co^) is positive or negative, and 

log (cOg — Oi^) — log (<>)j — %) = ^ TTt, 

according as I(oj| — co^) is positive or negative. Therefore the valnes of fx are given 



in tbe following table : 




negative 



I ((01 + (02) 

positive. 



ih ^ 



TTt 



(Oi 



fJL = 



Wl, 



0,).> 



I (Wi + Wo) 

negative. 



//. = 


TT/, 




0)0 


/^ =^ 


TTt 



Wl 



In other words, 

if 1 IToj^) I > I I(co^) I , ju. = ± TTi/ojy, the upper or lower sign being taken according as 
I(ft)j) is positive or negative, and 

if |I(w2)| > |I(o)i)j, /x = db 7rt/o>2, the upper or lower sign being taken according 
as l{o)^) is positive or negative. 

§ 33. The expression for cr (2;) in terms of simple and double gamma functions is 
now immediate. 

For on integrating the result of the previous paragraph 



vz'^' 



- log o-(z) = p ^ fiz -^ -=^ - \ogz + t log ^(2; ] ± wj, ±^3) - t log r^{z I ±«l) 

— SlogI\(r^|dbwo),, 
p being the constant of integration. 

Make now ,^ = 0, and we at once see that p -=- 0, 



jLJLt^xlvvv? 






HFg ^(-^1 i «!, i^o) 



and in this expression 



nr3~-i(0| ± cw^jTnvH^l ± ^3) ' 



i^ 



_^ w-t/ X f* ^ Cii)i ~p" (jdz 



■ 00 



\L' 






«e)Tft)o ft>i ^^*^ ft): 



the infinities in the double summation being equal in absolute magnitude, and that 
valne of m^ -^ Wo being taken which is on the same side of the real axis as m^ + ^% ; 
while 
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/x = 4: ~~ if I I(<^j)l > I i(<^'>)|j the upper or lower sign being taken accord- 
ing as I(w|) is positive or negative, 

and jLi = ± "" if I I(<^3) | > I I(<^i) I ? "th^^ upper or lower sign being taken accord- 

(J 

ing as I(<Wo) is positive or negative. 

§ 34. By means ot the preceding paragraphs we may now at once prove 
Weierstrass' relation'^ 

the upper or lower sign being taken according as I ( "~^ I is positive or negative, 
where -q^ and t/o are determined by the relations 

Take the expression for t^{z) given in § 32, and we find by use of the formulae 

of § 22 that 

*^ . 

l{z) — i{z + o)^) = ^- (m — m^ — m^ + mg) + z/^i, 

0)3 

where 

o)^ does 

mi = 0, unless i ^ r lie within the res^ion bounded by axes from 

^ ' o)^ — - o)^ does not J » j 

the origin to — 1 and co^, in which case m^ = i 1, according as •— l(o)^) is positive 

or negative, and m.2 and mg are obtained by changing the signs of (i) both o)^ and 

ft)j3 and (ii) o^i respectively in this formula. 

ihus -" = S 2 7c. — log "^ — " m •— m, — mo + mo , 

the infinities in the double summation being equal in absolute magnitude, and that 
value of o)j -^ 0)3 being taken which is on the same side of the real axis as [o)^ + ^^2)* 
And, similarly, 

where -^ = S 2 r^o — — log -^ h I""" ^^^ + ^^^1 + ^^h — ^^3 ^ 

where m! has its usual meaning, and m/ is obtained from m^ just as is m^ from m. 
Hence ^^ — . ^ z= - — [— (m — m') + (m^ — m,') + (^^^o — m^') — (m3 — mg')]. 

J, /-^ J- ^ 

^ JoimAN, ' Coiirs d' Analyse,' p. 351; Fousyth, 'Theory of Functions,' §129. Again notice that 
each of the quantities 1/ and w is double that usually taken. 
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But we have seen in § 22 that, if the difference of ampKtudes of co^ and Wo be 

greater than tt, 

m — m' = + 1 , 

^ — a— J 

-^ is positive or negative. 
Similarly it may be proved that if the difterence of amplitude of Wj^ and ■— w., 

is > TT, 

m^ — m/ = i 1, according as l(~"j is -^ve or —ve. 
And , . . , . , , "-- a)| and co., is > tt. 





m^ - m/ 


while, finally, 






m3 — 7n.^' 



ft)o \ . 



jt 1 • • ' . If ^j is -\-ve or — ^^^ ; 



^1/ 



0) 



• R ff f « • e 



and — Wo is > tt. 



1 ^^^^^ "^2 



ft) 3 



^ 1 s . , . If '- is -+^^6 or —ve. 



^L 



In all other cases the differences between corresponding nis vanish. 
But the difference of amplitude of one, and of only one, pair of the set 

can be greater than tt. 

Hence -- (m — m) + (^^i — ^lO + {"^h "^ "^h) "^ {"^^h "~ ^^3) 

must always equal Jb 1. 

And it is easy to see, by taking the particular cases which can arise, that the 
upper or lower sign must be taken according as l(o)Jct}^) is positive or negative. 

We have then finally 

™"' ' — • db 5 

and therefore rj^co., — tjc^co^ ^=^ Jt 2776, 



the upper or lower sign being taken according as If-" j is positive or negative. 

\ ]/ 

§ 35. For brevity, we merely indicate the relation of the formulie which we liave 
Ibund to the known relations : 



a 



O" 



(z 4- oi^) = — e' ^^ (T (z) 



By § 33 



TT 






W^'Czl + ft),, ± ft)o) 



« n^ www 
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and by the previous paragraph, the expression on the right-hand side of this relation 
will evidently reduce to the form e'^^^'^^\ where S^ is some constant whose value may 
be readily seen to be given by 

^1 ~" i^i^i I / \/l|l\, / \/l 1\ 

This value simplifies on detailed consideration. 
When jI(coj)| > i 1(6)2) I, 

according as I(coi) is positive or negative. 

When 1 1(^2) I > I I(coi) I , there are four subsidiary possibilities ;— 



(a) When co^ lies within the axes to — ■ 1 and — cog, 

O) J. ^1 ,, ,. ,, — 1 and C02, 

(r) '' — ^1 .. n „ — 1 and — C02, 

and (8) ,, — co^ ,, ,, • ,, — 1 and oi^. 

In cases (a) and (8) 8^ — |- 771 co^ = J^ Trt/cuj, the upper or lower sign being taken 
according as 1(^3) is positive or negative ; and in cases (y8) and (y) the upper and 
lower signs are interchanged. 

We thus see that in all cases 

so that we have the required equation 

cr (2; + wj) = — e"^'^' + ^'^^^ o- (a;). 
Similarly we find cr (2; + ^2) = -- ^'''^'^ ^ '"^'^ o- (2). 

The verification of these results affords substantial proof of the general correctness 
of the signs which are involved in the work. 

§ 36. It is interesting finally to notice that just as the gamma functions do not 
exist when r = o}J(o^ is real and negative, so the elliptic functions demand that r 
shall not be real. 

The condition that r must not be real and negative arose explicitly at several 
stages, and might have been predicted d priori. 

For, when oi^joi is real and negative, it is obvious that 



O = m^o)! + ^^^2^2 ' 



^ ' ? • • 5 I excluded, 

Lm^ = 0, 1, . . , 00 OJ 

will have a zero value at least once. 
And thus the function 

Tg ^ (2;) = e ^ . z . li n 



^^ o "^ 202 



1 +-^)e 

Till =0 Wg = L \ ^*/ 

will be infinite independently of z ; that is to say, it ceases to exist 
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For all other values of m^ and co^ the function r^""^ (^|^i? ^2) exists. But the 
product 

and consequently or (^), will not exist when either M^/cf^i or — w^/cu^ is real and 
negative ; that is, when r is real. The criterion for the existence of multiple gamma 
functions (^2-ple where n is greater than 2) is more intricate, and, as we know, 
n-ply periodic functions (n > 2) do not exist. "^^ 

Part ITT 

Contour Integrals connected 'with the Double Gamma Function. The Double 

Riemann Zeta Function. 

§ 37. In the theory of the simple gamma function it vfas shown that the interven- 
tion of a definite integral, coupled with the theory of asymptotic approximations, it 
was possible to obtain contour and line integrals to express Eulek's constant y, and 
the logarithm of the simple gamma function and its derivatives. We now proceed 
to show that it is possible to extend the method thus previously employed so as to 
obtain expressions as contour and line integrals for the gamma modular constants 
72^ and y^^, and the logarithm of the double function and its derivatives. It will be 
noticed that when the real parts of coj^ and w^ are positive, the numbers m and m' 
which intervened in Part II. vanish, and there is consequently a noteworthy simplifi- 
cation of the formulae obtained. This simplification extends also to the definite 
integral expressions, and consequently we shall first investigate the theory in this 
simple case, proceeding subsequently to contour integrals of greater complexity. 
Finally we make use of an extension of Mellin's method of defining the simple 
C function by a series instead of a contour integral, and we show that there is 
complete agreement between the formulae obtained in the different ways. 

§ 38. When the real parts of co^ and w^ are positive, and when in addition the real 
part of a is positive, we define the double Riemann ^ function 

for all values, real or complex of s, by the integral 

iT ( 1 - s) f ^jz^^i- ^^ zlM^^-^^^ 
lit J(i- g'^i^xi - g-'^'^^o ' 

wherein (—^)^~^ = e^'"^^^°^^-'^, log (—2;) being real when % is negative and being 
rendered uniform by a cut along the positive direction of the real axis. The 
integral is to be taken along a contour enclosing the origin (but no other pole of the 

"^ The existence-criteria for functions which are substantially multiple gamma functions have been 
discussed bj^ Ceani, * Batt. Gior.,' vol. 29. 
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subject of integration), and the positive half of the real axis ; and extending from 
+ 00 to + 00 as in the figure. 




Under the limitations specified the integral is^ in general, finite^ Moreover, by a 
theorem previously obtained,"^'' we have under such limitations 



^r(l -s) 



9, 



IT 



j / — zY~^ e ~ ^^* "^ "^^'^^ "*" ^^*2"2>^ dz ~ 



1 



(a + 'i>\'^i + ^3^3^ ' 



the latter expression having its principal value. 
And therefore 



mi = m,^Q (^ -f %ft>i + ft>3G>3> 



tV (1 - s) 
2^ 



X — g-' ^v'l'^ + i)<^i* 
1 — ^-'^i'^ 



^,-{qn + \)(a.^z 



e~''^(— z)' ^d: 



z 



r I \ ^r(l-s) 



1 _- er^^ 



(1 _ ^-«i'^)(l — C-«/^) 

4- — >- L I — --^^ — -- ^--- / '7V"^c^2J 

27r J (1 — e-^'i') (1 — (?-«^^) ^ "^ ' 



- (---^)^-y, 






all the integrals being taken along the fundamental contour. 

When n is a large positive integer, we proceed to throw the two integrals last 
written into the form of asymptotic series. 

For this purpose consider the expansion obtained in § 16, which may be written 



(1 - 6-"!^) (I - e-'-^'O «if03^' 



1! 



?l 1 



We showed that this expansion is valid provided | z \ was less than the smaller of 



the two quantities 



2in 



(Oi 



9. 



TTl 



(Ocy 



Outside this circle the series diverges. But within the region bounded by lines 
going to infinity from the poles of 



<;/^ — (« + 0Ji)3 

(1 — 6-^1'^) (1 — " e-""-^') 



1 



CO I CO cy^^ 



^ '' Theory of the Gamma Function," §§ 22, 33 and 34. 
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TTt 



the series will, in the language of M. Borel,"^' be summable, that is to say, within 
this region (which is conveniently bounded by straight lines from the |)oints i 

jh- — which pass to infinity through the remaining poles) it is possible from the 

values of terms of the series at any point to obtain, by the employment of 
intermediate functions, a magnitude, independent of these particular intermediate 
functions, which is the value of the function at the point. 

If then on any term auZ'^ of such a series we perform the operation which is 
expressed by 

^r(i-s) 



2tt 



CL^j/C ( """^ /6 ) Ct/O., 



we shall expect to obtain a quantity which is the r^th term of a possibly, and even 
probably, divergent sequence, which in turn is, by suitable operations, summable to 
the value which results from the performance of the fundamental operation on the 
function of which the original series is the expression. 

Such considerations being understood to underlie the operations, we have 



^r (1 - s) 



" 0" [Ph^ + 1)^1 + c/]z \ ^- [q{m + l)wa + alj 



9 



TT 



(1 — e""^i^) (1 — 6-"^^") 



(— zy hlz^ 



tr(i-H)f(-~~ zy-^ 



9. 



TT 



«1^2 



'^-pni»iZ I g~aft«2^"] ^l. 



z 






iV(l ~~s) 



9 



lilt 









-|-. -^ -^ / — £\S-^^-~(l>ii^-£i J 1 3 _|_ 3 (^^.^^ 

2w" J [ 2ft)j^coQ ^1% 



+ ^% ^ t V^--^^^ {^\n{ci + oy^er^'^-^ \ ,S',(a + o,^ e'^^''-^^\dz 



1 



1 



(5 - 1) (s - 2) 0)10)3 l(i>^a)i)--^ "^ Im^'^l 



1 0)^+003 + 20- 



s — - 1 



2o)j^O)c 



1 



1 1 



ww..i *w «» 



1 



1 



S — 1 [o)i " (pnft)|) 



1 11 



o)o (qnw:,y 



s—\ 






9>l=l 



JL « *^ * • « I f V 



Now it may be readily deduced from the results obtained in Part I., that 



t 



S;/(a) =: w[jjS«„|(a + a>,) + pB,,, 



•^ BOREL, * LioLiville,' 5 Ser,, vol. 2, pp. 103 et seq. * Annales de Ecole Normale Superieure,' 3 Ser., 
vol 6; pp> 1 et sea. 
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Hence we have 



tr(i — s) 
1 



(1 _ e-™!^) (1 — e-'"^^) 



1 






+ 






(S-1)(S- 2)ft)ift)2 L(p9t«iy-2 (^«,a,3>'-3j 



(— z)* ^(iz 






2ti + fflj + (Bj 



S- 1 • 



0^10)3 



1 



.(F^«i) 



s-l 



1 

•< 



1 



a)l(p^cOl) 



^ + 



(o^iqnoy^ 



s-l 



>S /(fX 4- ft)i) _^ 38/(^6 + 6)3)1 






3S,;,(<rj^ + 0)3) + oB 



Wl-fl 



(^'?ia)3) 



m-\-s 



In an exactly similar manner we find, since 



that 



^P n __ A r ^-[a + (im+l)(Oi + (fi/H-l)wa].-; 



2 



(1 — e-'^i^) (1 — 6~*^2.') 

1 1 



2a 4- ft)i + ft).7 



(s — 1) (5 — 2)ft)j^ft)3 * (p^ft)i + q^^o)^^ ^ (s ~ l)2c»^ft)3 * (P^^^i + ^'^^^y ^ 



3S/(a 4- ft>3 + 0)3) , ^ / y^^ fs + m - 1\ ^S^(a + co^ + co^) 4- gB^^+i 



(^7^6)3 4- g'^^Wg)' m=l 



m 



(pnco^ 4- ^^0)3^ 



If now we group together all the results which have been obtained, we find the 
asymptotic quality, true for all values of s, and for values of a, co^, and co^, whose real 
parts are positive, 



1 



p7i qn 

^ ^ — 

wi=o m,=o {^^ 4 micoi 4- nio^co^y 



= U{^, a \ 0)^,(0.2) 



WM**WM 



1 



(S- 1)(S~ 2)0)10)3 

1 1 



1 



(pnco^ 4- Q'^^o)^y ^ (pnco^y 



_o 



(qnco^y 2 



2(re 4-0)1-1-0)3 

2(s + l)o)iO)3 i(imo)i + ^7^o)3)^-l (F^ft>i) 



5-1 



(qnco^) 



s-l 



S 



~ 1 






.S/(^ 4- o)i 4- 0)3) 2S/(a 4- ft)i) gSiX^ + «»3) 



.«i(P^«i) 



r-l + 



<»3(2'^«2)' ^ 



(pncoi 4- g'^^s)^ 



(jmo)!^ 



(^^0)3)^ 



00 






(->'^ /'«i 4- 5 - l\j 



m=l ^^^^'^^ \ ^^^ / 



38,^ (g 4- o)i 4- 6)3) 4- 3B 



W4-1 2^?»(^ 4- ^i) 4" gB.;^^! 

(|?0)i)*+^ 

2S^(a 4- 0)3) 4- 3B 



»<+l 



(q^i) 



■m+s 



r • * (xL) 



It may readily be seen, just as for the case of a single parameter, that the series 

. 1 . . 

proceeding by powers of is a series of powers of a real variable, whose line of con- 

it 

vergency is of zero length. 
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This series is summable by an almost evident modification of the appHcation of 
Borel's ideas, which was employed in the ^'Theory of the Gamma Function/' § 38. 
It is thus the asymptotic equivalent of the sum 

2rii qn 1 

^ ^ - --.-- . 

and it satisfies Poincare's^^ criterion for asymptotic equality that the difference 
between this sum and the first m terms of the series has its absolute value less than 

1 
a quantity of order t^^^^.,^ • 

§ 39. The function t,^{s, a\o)^, ojg) has been defined, and the asymptotic equality (A) 
has been deduced only for the case in which the real parts of a, oj^, and co^ are all 
positive. 

It is natural to try and use the equality to define the function for all values of 
a, W|, and Wg. 

In the first place, it is evident that when U(w|) and ^{(cog) are both positive, the 
equality (A) holds for all values of a, for the various terms of the sum and the equiva- 
lent asymptotic series are continuous for all but an enumerable number of values of 
a, (0^, and (o^. Hence in this case 'C^i^^^l^i^ ^2) i^^^J be defined as the term inde- 
pendent of n in the equality. 

So also when s is a real positive or negative integer, the equality will hold for all 
values of a, co^^, and o)^. 

But when s is not an integer, the various terms involving s in their index are 
multiform functions, and to ensure uniformity we have to assign definite cross-cuts 
to the logarithms which arise in the equivalent exponentials. When, as under the 
limitations for which the equality (A) has been established, these cross-cuts are 
formed by a line outside the smaller angle between the axes of co^ and oj^, the expan- 
sion is perfectly valid ; but when the common cross-cut lies within this angle, terms 
arise similar to those which occurred in Part 11. of this paper, which are multiples of 
27n, and involve n. 



Cro&3-cut 




/J 



Cross -cut 

O 




m^ 



And, therefore, if we attempted in this case (see the second figure) to define 
^c^{s,a\o)^, 0)^) as the absolute term in an asymptotic equality such as (A) § 38, where 
for complex values of s the principal value of each term is taken, we should ultimately 

■^ PoiNCAKE, ' Acta Mathematica,' vol. 8, pp. 295-344; ' Mecaiiique.Celeste,' vol. 2, pp. 12-14, 
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find that Czi^n, ajcoj, ojg) as so defined would involve n. In other words, we should 
have made an assumption which could not be justified. 

If we wish to obtain an expansion valid for all values of co^ and oj^, we must con- 
sider as our starting point the integral 



.r(-.) 



^-az ( _ zy-hlz 



TV 



Jl(1 — er^^^^ 



) (1 — e-"2^) ' 



where ( — zf ^ = e^' ^^ ^""^ ^ '^ , the logarithm being rendered uniform by a cut along 
an axis L, coinciding with the bisector of the smaller angle between the axes of 



to 
1 



1 



and — ^, and where the integral is taken along a contour having this axis L for axis 

(as in the figure), and enclosing the origin, but no other possible pole of the subject 
of integration. That value of log(— 7^ is to be taken which is such that the imagin- 
ary part of the initial value of log (— L) lies between and — 27rt. 



&), 




This integral of course is only valid when a lies between the smaller area bounded 
by the axes of w^ and cog, or, as we may say, when a is positive with respect to cj^ 
and cog. We notice that the line L is uniquely defined, since the ratio cog/oii cannot 
be real and negative. The definition of the integral is not complete when 0)^ and ca^ 
include and are equally inclined to the axis of — 1 ; in this case we may take L to 
be a line nearly coinciding with this axis. 

We now define the double Riemann ^ function, when the variable a is positive 
with respect to the w's, and 8, w^, and co^ have any complex values, by means of the 
equality 



4 {^^ ^ 



CO I, 6)2) 



iV{\ -s) 



2 



e 



2Msiri 



TT 



f g-^^ ( — z ) ^-1 dz 
J L(l — e-**'!^) (1 — e-"2^) 
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where M = 0, unless the axes of -- and (coi ~\~ co^) include the axis of — 1, in 

E 

which case 

M = ^ 1, the upper or lower sign being taken as I^co^ + ^2) ^^ positive or 
negative. 

§ 40. Let us take now the integral which has just been defined, and apply to it 
the procedure of § 38. 

We shall evidently have to consider integrals of the type 



2mr 



I er'''{-- z)''~hk, 



where the axis of n lies within the smaller angle between the axes of o)^ and co^. We 



can at once see that this integral 






1 e"'"-^( -z)''^hk 



n 



where u = 0, unless the axes of L and — embrace the axis of — 1, in which case 

' n 

jLt = i 1 , the upper or lower sign being taken as I(n) is positive or negative. 

Let n = r&\ co^ = ae'"-, o)^ = be'^, where 0, a and ^ are measured between and 
2tt by rotation in the positive direction from the positive half of the real axis. The 
axis of L proceeds from the origin to the point e""2^'' + ^^ and therefore where z is at a 
distance p along the axis of L, 

nz = rp 0'(^»-~^)> 

This quantity has its real part positive when 

_ ^- ^ /9 _ ^^ ±A ^ 5 

2 ^ 2 ^ 2 ' 

a relation which is satisfied when 6 lies between a and ^, and the difference between 
a and /3 is less than tt. It is also satisfied when the axis of L proceeds from the origin 

to the point e~'("'2~~"^ V ; where € is a quantity less than half the excess of tt over 

a ^^ ^, We see then that the axis of - lies within a range of a right angle on 

either side of the axis of L. Hence by the propositions previously proved ['' Theory 
of the Gamma Function," §§ 33, 34], 



z 
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unless the axes of L and embrace the nep:ative half of the real axis. In this latter 
case, since the imaginary part of the initial value of log ( — ■ L) lies between and 



27rt, as also does that of log 



1). 

n r 



we are giving a different prescription to the 



many- valued function which occurs in the subject of integration. 
We therefore have 



"^^^ '^^\^e-'^\- zY~\lz 



'TV 



~"^~?"r^-' e^''"'|^e-«(- zy-'dz 



tlT 



n 



where ft has the values which have been assigned to it. 



Now 



'^ ^^ '^ ^1 ^"" ( - ^y~' <^^ == ^' '"' ^'"" '^: 



'}. 



TT 



n 



the logarithm having its principal value {'' Theory of the Gamma Function/' p. 107). 



J LtJXXvytJ 



ir(i - s) 



'iTT 



L 



e ""^ — z)' ^dz = e 



■s [log ti + 2ix^Tn] 



1 



where /a' = 0, unless n and - embrace the axis of — 1, in which case /a' 

XJ 

upper or lower sign being taken as I( -^ ) is positive or negative. 
Finally then 



± h the 



27r 



e '^' ( — 2:) ' ^ dz == 



1 



n" 



where the latter function = e "^"s^^ where log 11 has a cross-cut along the axis of 
— Y, and is real when n is real and positive. In other words, log n has its principal 



L 



1 



value with respect to the axis of - , . 

XJ 

§ 41. If now we apply to the integral 



cVil — s) C e-~«'^( -- zy-hh 



9, 



TT JL(1 - e"''^i^~)(l— ^-'^sO 

the procedure of § 38, we shall, for all values of 5, a, cd^, and co^^ such that a is positive 
with respect to the ms, obtain the asymptotic equality 



2Jn qn 



1 



iV {1 — s) 



coi = m, = (« + ^/^i6>i + ^%W3> 



9 



IT 



6~«^(— z)'~hlz 
L(l — ^-'^1^) (1 - ^-'^a^) 



wii I HUM 



1 



n'~^ (s — 1) (s — 2) w^wj 



. iP^i + ^^2)' ^ (P^i) 



5-2 



2a 4- Wj H- W3 

2(s — 1) w^w^ [(p^^i + ^^'^'Wg)^"'^ (P^<^i) 



1 



5-1 



($^^^3)' ' 



[ovee] 
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1 



I 



1 



S — - 1 [ 0).^ (i^^'^^i) 



A— 1 






1 



o)^{qri(o^' 



if + 



gS^^C^-*^ + ^1 + ^^-i) 2^\(^^ '^ ^-^l) ^^\(^ + ^^2) 



{jmar^ + (Z^^^g)' 



a. 



-f- S 



, /iv'/M -\- s 
111 ^ 1 /^ 






V <v \ JL c^ /O 



y/i + 1 






(f/lwg)^ 



;;i + 1 



{V<^i + ^^3)"' 






(^i«2 



ill + A' 



wherein all the many -valued functions with .s' as index have their principal values 

1 



with respect to the axis of 



L* 



It proves convenient to consider these functions as 



having their princijjal value with respect to the axis of — [o)^ + co^). In order that 
this may be the case, we must multiply the integral by 



6 



where, as in § 39, jj. ■=: 0, unless (w^ + ^%) and ^ includes the axis of 



1 

L 



L in which 



case M ::=:: i 1? ^s I(co^ + oj^ is negative or positive, 

liemembering the definition of t,^{s, a|(i>|, ojo) given at the end of § 39, we see that 
we obtain for our fundamental asymptotic equality an expression which in form is 
identical with (A) § 38, but in which the many-valued functions with .s' as index have 
their principal values with respect to the axis of — (tt)^ + ^2)' '^^ ^^ evident that 
the equality will hold for all values of cf, and will thus serve to define l^c^i^s^ a\(i)^^^ o)^ 
for all values of s^ a, co^, and ojc^. 

§ 42. We proceed now to take such particular cases of the generaJ asymptotic 
equality which has just been obtained as lead to expressions for the logarithm of the 
double gamma function and its derivatives. 

Suppose that s is a positive integer greater than 2 ; then, making n infinite in the 
general asymptotic equality, we see that 



( — y 



a 



0) 



1. ^2): 



where 



d/ 



^2^'^ (a 1 0)1, o)^) =- , dog J\{a I to,, 0)^) 



This relation is true for all values of a^ co^ and co^ ; it is the first of a series 
connecting the double zeta and double gamma functionSe 

Let us next put 5 + ^ for s^ where e is a small real quantity and s is, as before, a 
positive integer greater than 2, Then, provided a is positive with respect to oi, 
and coc 



'3? 



^S + 6, 



a 



0J|, (i) 



).- 



a^Xi 



S 



e) 



2 



e 



./^MTrL{s+e) 



TT 



L (1 ^e~^i'){l ~ e~^^^) 



the integral being taken along the L-contour, and M being the integer defined at the 
end of S 39. 
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Hence if log (— z) is real when z is negative, and is rendered uniform by a cut 
along the positive direction of the axis L, 

""■^ "^27r' (5-1)!' IJt, (l--e--^0(l-^~^^^^)r ^ '^^ ' ' ' 



X \ 1 



1 1 

€ ( ^7 + ... -f " ~ 7 

1 S — 1 



1 



. . . U 1 + ^ ^og ( — 2; ) + • 



X il + 2M7ne + . . . 1- 



for 



( - y-' 






1) 



r(6) 



(S - 1) ! 6 



1 — ye \- 



(e - 5 + 1) (e — 5 + 2) ... (e — 1) 



1 1 \ 

1 -|-- € ( ~ p . . . -f" 7 ) 



o 






• • • ]• 



ow when s is an integer greater than 2, 



JL {1 — e- ''^') (1 — e~ ''•^^') 



0. 



for the integral may be reduced to two line integrals which destroy one another 
and an integral round a small circle enclosing the origin whose value is zero. 
We have then, on making € -^ 0, 



C{s,a\a);, C03) 



1 1 

t (-y-i r 6-«^(-«)»-i|-log(-a)+ 2MTn + j + . . . + — -^ - yj 



27r (s - 1) ! 



IJ 



(1 — e~^l^) (1 — 6-^^22) 



<i2;. 



Butj when 5 is an integer greater than 2, 



e-«^(— ^)-^~"^J2M7r/^ + -.- + 



I 



1 



. -f- 



1 



1 



7 



L 



(1 — ^-^1') (1 - e-""^') 



dz 



vanishes for the reason just assigned. 

We see, then, that when a is positive with respect to o)^ and co^, and .s" is a positive 
integer greater than 2, 

1 r e"'''(— zy-^iog{— z) 



^'% o^ 



^i^(^^)={'^y 



277 (s " 1) ! J L (1 - e-^^') (1 - e-^^^') 



-dz^ 



and therefore under the same conditions, 



^a^'-(^ 1 ^2. ^^2) 






cfo. 



43. Put now in this result 5 rr. 3 ; then with the assigned limitations 



^%^^\^A^i^ CO,) 



27rJ 



L (l-e~^»i^)(l - ^-^^0 



cfo. 
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and it is obvious that we may replace log(— z) in the subject of integration by 
log ( — 2;) + y without altering the value of the integral Integrate successively 
with respect to a, and we obtain 



log T^{a\ 0)^,0)^)1-:= 



27r 



L (1 — e"'''^'){l — c^'^-^^') 



dz + (Ij ^) 



where the coefficients of the additive quadratic form are constants with respect to a, 
Remember that 



LP* ■ / 

6 Clli^ ict 



then we evidently have 

log r.2(a 1 0)^,(^3) 



a = 

2nT 



^1^ ^2). 






,-«..(^,.)-i|log(-™,^)+^} 






^- «*■(__ 0)-"i{log (^ z)+ j} 



dz + Xj_a + \^a 



2 



L 



(1 — e^^i^)(l—e""^2^) 



+ log a 



We now define the third double gamma modular form p^{coj^ co^) by the relation 

logp^{a)j, tOa) =:: ^ 2Mm ^S'-^(0 I 0)1, 0)3) 



—'" Jut 

a = 



and we proceed to show that the constants X^ and X^ are such that 
log t^-^^1^1^ = 2So(a) (M + m + m')2^t + 2^'i{o)2M7n 



-|- 



^TT J 



' ^•"- -(-^rHiQg("-^)+7} 



dz^ 



where the numbers m and m' have the values assigned in Part II. 



If this relation is true we shall have 






— • B\{a I %) (M + 071 + mf)27n + ■; 



e r c-^^^X — ;s)"'i{log(— ^) + 7} 



27rJL 



1^0-^2^ 



XAjZ, 



§ 44. Let us now consider this integral 

It is to be taken along a contour embracing the axis L, which we take to be the 

,11 
bisector of the smaller angle between the axes of -™ and — , unless such bisector 



ft)i 



should be the axis of — 1^ in which case we take it to be nearly coincident with this 
line. And in the subject of integration that value of log (— 2;) is to be taken wliich 
is real when z is real and negative, and is limited by a cross-cut along the axis of L. 
Let us consider the relation of this integral to the integral 

^ r 0~az{^_ 2;)~~^{log(— z) H- 7} 



2ir 



Wg 



{l—er^'^') 



dz. 



which is defined in the same way with reference to the axis of l/rjo^. 
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We proceed to show that 



c ' 



2 



TT 



L 



(On 



where jjl' = 0, unless the axes L and 1/c*)^ do not enclose the axis of -— 1, in which 

case 

fji' =z rfc 1, according as I(ft>g) is positive or negative. 

For take the integral along the contour embracing the axis L, and suppose the 
contour to expand so that it embraces also the axis of l/cog. 

Then since a is positive with respect to o)^ and cog, and since the angle between 
the axes of L and l/o^^ is less than Jtt, the value of the integral will be unaltered, 
for its value along the part of the great circle at infinity between the axes of L and 



l/wg IS zero. 

Suppose now that the contour is taken to lie on the infinite-sheeted Neumann 
sphere, whose sheets intersect in the cross-cut from to oo, on which the subject of 
integration of the integral is uniform. We may, without altering the value of the 
integral, deform the cross-cut so as to take up a position along the axis of l/co^, 
instead of along the axis L, provided that in doing so we do not give a new specifica- 
tion to the logarithm. The latter phenomenon will occur when l/w^ and L embrace 
the axis of — 1, in which case we take the first contour in a sheet in which z can 
assume real values, while the second is taken in one in which log {— z) for real 
negative values of 2; is equal to ± 2wl 

After deformation of the cross-cut we may compress the contour so that it 
embraces the axis of l/co^. It is easy to see by this repetition of the argument 
previously employed in the ^^ Theory of the Gamma Function," that we have 

c f er'''(-z)-^ {lQg(- ^ ) + 7} ^^^ 



27r JL 1 — e-'^^^ 

I r g-a2(— ;^)-l jlog(— z) 4. 2/iVt + 7} 

1 



- 1 — e.-^^"" 

«2 



27r 

i_ f e--«^ (— z)-'^ {log (— z ) + 7} 
27r J- 1 — 6-*^2^ 

«2 



dz 



dz — 2/xVt S'l (a I 0)2) 



where fi has the value previously given. 

§ 45. The assumption made in § 43 for the values of the constants Xj and X^ will 
therefore lead to the relation 

^°S -T,- JP^ = ^^g ^i? ^'^ (« h.) {M + m + m' + /.'} 27rt, 
for (" Theory of the Gamma Function," § 37) 



T^{a\(o) I [ e-«^(_2)-i{log(-^) + 7} 



log ( \ — c, 



I 1 - 6-"^ 



dz. 
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Now, by considering the various cases which can arise, it may be readily seen that 

~-|- lib "j— /x — - 



These constants all vanish unless oj^ and ojo embrace the axis of — 1. When 
this take place, suppose that cd^ lies above, and o)^ below, the real axis. 



Then 



[J. 



ft 



=: 0" 
zz 1 



when (oj, + coo) lies above and 1/L below the real axis. 



m 



- ■>when (co^ + oi.,) and 1/L both lie below the real axis 



>when [(jd^ + ojo) and 1/L both lie above the real axis. 



IL 



:==, 1 
:==: 



m ::=: >when (ojj + co^) lies below and 1/L above the real axis. 



:rr 1 

We get similar sets of values when the imaginary parts of m^^ and o)^ have opposite 
signs to those just assumed. 
In all cases 



and, therefore, with the values assumed for X, and Xg? 






Vf' (^^) 



Pi (^^ 



Similarly we should find 



To"! (a + CO,) 



Pl(^l) 



But these are identically the fundamental formulae for the double gamma function, 
found in S 23. 

The values assumed for X. and X^ are therefore correct. 

We have, then, the two important formula3 



log 



P3 («!' ^3) 



2 



So (a) (M + m + m) 2nL + gS'i (0) 2Mm 






f (r^^' (— ,a;)-i {log(-- 'Z) -h 7} 



o 



/j7r 



T 
JU 



(1 — 6"^!^) (1 — g-^3^) 



^&, 



and log p2 (co-^, ojg) =^ — 2M7rt gS'i (0 | co^^ o)^) 



Jut 

= 



9. 



IT 



L (1 — e~"i^) (1 — e™^2') ^ 
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which express as contour integrals the double gamma function, and the third double 
gamma modular form. The first relation only holds when a is positive with respect 
to the oi's. The second is valid for all values of o)^ and cog, subject of course to the 
dominant condition that c^i^/w, is not real and negative. 

It is worth noticing that the first formula) may also be written 

l^g Mii») ^ ^s^ (^) (,,, 4. ,,^') 2in + ,S\ (a) {2M7rt + 7} 



P^i^v ^3) 



■h c 



'X 



7r J 



L (1 — e~ "1^") (1 - e~^--^) 



\A/4/t 



§ 46. Subject to the condition that the real parts of a and L are positive, we may 
now express our contour-integrals as line-integrals. 



Consider the integral 



27r 



e ^^^ ( — z)~'^ lo£? (— z) , 



By hypothesis the logarithm has a cross-cut along the axis of L, the initial value 
of its imaginary part lying between and — 27n. Hence if the contour of the 
integral be reduced to a straight line from 00 to €, where e is a point on the axis of L 



very near the origin, a circle of small radius 
line from e back again to -f^ 00, we shall have 



round the origin, and a straight 



t f e «^(— 0) '^log('-z) 



fco 



2w J L (1 — e~ ""y- {l — e" '^^) 



dz 



e^^'\l 



e~«^(— zY^d. 



e-^i')(l — e-'^s^') 



W«ILrAiM<W 



1 



r27r 



c - ^^^^^^ { log e + i{e -tt )] ^0 



27r .0 (1 -^"'^0(1 



e 






) 



The logarithm in this second integral, which results from the small circular contour 
surrounding the origin, has its principal value. The integral itself is evidently equal 
to (§ 15) 



1 



r2ir 



9. 



TT Jo 



d^[log e -^ i{0 — tt) _ 



g-2.0 






,^^''K») .. - .e j^ i>jM_ 



e 



€ 



jL. * 



+ terms involving positive powers of e 



2S\(a)log 



1 S ^'^Ua) 

I ■+• ^-^-— ^ -- J^ terms which vanish with 1 e 



'Iw^oi^e^ € 



Thus 



2 



e 



^ - «^ ( — ;>;)" ^ log ( — z) 



TT J 



L(l 



— w 



'') (1 



e 



to.) 



') 



dz 



X) 



6-«^(— z)~'dz 



g(-L)(l __ e-<^^'){l — (J-^^^) 



+ gS'i (a) log e 



1 



O 2 



4' 



Sf> (a) 



+ terms which vanish with lei 
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Now, as has been seen in the *' Theory of the Gamma Function," § 28 cor., when 
15. (L) is positive, 



/'OO 



(L) 



e 



ch 

X 

z 



log e — y + terms which vanish with 



And evidently 



/»CX3 7 



y«00 



(L) 



cfe 



2; 



2e^ 



The integral under consideration is thus equal to 



7 

(L) 



6 



-fls; 



,^ [ (1 •— g - ^^i^) (1 — e ~ ^^^) CO^Wg^' 



1 +^-^-^-^-e-,S\(a)^~y,S',(a) 






+ terms which vanish with 



If now we make e coincide with the origin, the integral last written remains finite 
and we have 



'm 



6 ^^(~ z) ^ log { — z)dz 
L(l — e-*^^') (1 — ^-^2^j 



.00 



(L) 

Jo 



dz 



e 



~~a& 



Z 



\{l — e-""^') (1 — e-^'^^) 



(Oj^moZ 



1-^ + A^^K^O 



^ 



e ''2S'i(a)|- rzS'iH- 



This equality may equally be written 

]^a;^M^i^ - „So(a) (M + «i + m') 2^1 + oS'i(o) 2Mm 



/. CO „_, 



2ir J 



z 



t ' e ^^^(— ;a;)~^{log (— z) + 7} 7 
L (1 — ^-^1^) (1— g~"«^) 

dz 



f (L) 



e 



'rt;^; 



•'0 



Z [(1— 6-'^l«)(l— £J-^2^) 



.Si^!M + ^^li^ 



4/ 



e-%^\{a) 



J 



under the assigned conditions that a is positive with respect to the a>'s, and that the 
real part of L is positive. We thus express the logarithm of the double gamma 
function as a line- integral. 

In order to obtain a line integral for log p^ (w^, co^), we notice that we have 



.00 



log a = (e"^^ 



e 



'a. 



. dz 
A - 



and therefore 



27r 



e~~^^(— 2^)~^ {log (— z) 4- j}dz 



,oc 



= (L) 



dz 






L (1 — 6-^i~)(l — e. 

(1 _ ^~'^i^)(l — ^--^a^) 



Wgg 



) 



e 



•fl^ 



-|- log €1 






+ e-il-,S»j 
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Therefore on making a = 0, we have by § 45 
log p{o>^, cog) = — 2M7rt 2Bi(«i, %) 



'00 



(L) 



dz 



1 



z [(1 — e~'^i^)(l- e-""'^') 



'^■"'<"-' + A^"^-"-' + «-=[l-,S',(a)] 



Z" 



On differentiating the formulsG which express the logarithm of the double gamma 
function as line and contour integrals, we obtain 



xfj^^{a\(o^,(o,^ = 



27rj 



e "-^ log ( ~ z) dz 



.00 



= l.(L)'fia 



L (1 — (!~"i-)(l— er"2^) 



+ 3S'o(«) 



7 + (M + m -f- ^^') • TTt 



g-C02S) (1— ^-'^■2^) 






e-^^ gS'o (^^) r + 2S'o (^0 [(M + ^?^ + m)27ri\. 



Similarly, again differentiating, 



i^%^\<^\(^i. (^^) 



2iT J L 



f. «« (— ^)|log (-- z)}dz 



6-coi^j(l™ g-co,^) 



- + .S7'(«) 



y + 2771 (M -f- m -f- m' 



00 



(L)zdz 



e' 



-az 



(1— e-'^i^)(l— e-"'^') 






(M + m + m') . 27rt 



And, if ,s be greater than 2, 



\\)^^\a\ (x)^^ (Oc^) = 



t r e-^^X— ^)^_~iiog(- ^)cZ, 



.00 



(L) 



^"«^(— ;2;)^~^ 



0^2;. 



27r J L (1 — 6"'^!') (1 — e"'^2^) J ^ -^ (1 — e-^i^) (1 — g-'^^^) 

Notice that, when we have the more narrow restrictions, the real parts of a, 
co^, and co^ are all positive, the constants m, in, and M all vanish, and there is a 
substantial simplification in the formulae. 

§ 47. We may now deduce expressions as line and contour integrals for the first 
and second double gamma modular forms 



731 K. ^3) a^nd 700(^1,^2). 



We have seen (§22) that 



1 



xjj^^'^ia I c^i, 0)3) :=:z 723 (a)^, co^ + ay^^{o)^, CO.^) + ;; 



00 00 

"T"" 

and, therefore, on making a :=r. 0, 



so that, by the last paragraph, 



1 



+ 



a 



a -4- m^coi + '^^^^3 7?i^ft)3 + '>?^2^2 (^^^i<Wi + "^^h^^y 



xjj^^'\a\a)^,(o.)^- ^ 



? 



a^:::0 



00 



J 



VOL. CXCVI. — A. 



1 



(1 — e-^^'){l— e-"-'^^') 



— 1 



A^'Ko) ^^l!^g^(2)(^>^ 



z^ 






gS'o (0) 27ri(M + m + m'). 



2 u 
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The additive term will of course vanish when the real parts of 0)^ and co^ are positive, 
Similarly, we have 



721 ( 



0)^, O)^) 



xlif^\a\o)^,ox^ 



1 



a' 



a~0 






2 



ttJ 



G~-- (-~~z){ log ( - z) + 7 } J 1 ^ 

L (1 — <3~'^J^)(1 — e""^-^') c(? 



' gSo^^^^(a)27rt(M + m + ni), 



so that the first douhle gamma modular form is expressed as a line integral by the 
formula 



721 ( 



^00 



O)^, Wg) 



J 



1 



^0 



(— Z)dz\ 77 ,7/ 

^ ^ [(1- C"'^i;)(l 



g~0>2^j 



1 



e 



3 2Si^^^(0) 



It will be noted that for the modular forms 



^2(^1. ^2)5 72i(<^]^ ^2)3 722(^^1^ %)3 

we have, by making a vanish, obtained, line integrals which are in general finite, 
although in our fundamental formulae the restriction was made that the real part of 
a should be positive. 

This restriction was necessary to ensure that the contour integral should be finite 
at infinity. It is clear from the mode of generation of the line integrals, that the 
process which has been carried out is perfectly valid, since by the introduction of the 

terms log a, — 



05 

ay a 



allowance has been made for the manner in which the contour 



integral tends to an infinite value as a tends to zero. 

§ 48. At the beginning of § 43 we entered on the investigation which has just been 
given by integrating with respect to a under the sign of contour integration, and in 
this way we deduced the contour integral for log Vj^a) from that for \jj^^'^^ (a). 

We now proceed to show how the contour integral for log r2(a) may be obtained, 
without the employment of this process. 

For this purpose we take the fundamental asymptotic equality of § 38, valid for 
all values of s^ a, w^ and w.., the many -valued functions with s as index having their 
principal values with respect to the axis of — (co^ + ^2)^ 



ji9l q^l 1 



4 ('^? ^ I ^n %) + 



1 



{s -^ 1) {s 



2)a)2&)o 



1 1 



(j^ncoiy '^ (qnoy^Y"'^ 



2a + ri)j 4- ft) 



3 



2 (5 ■— 1) (OiO)^ 



1 



1 



1 



1 



(pnco^ + q^co^y 



s~l 



{pnw^y 1 {cpiw^y 1 



s — 1 



1 






1 



[over] 
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,Bj (a + (0^ + a)g) gS/ (^-f ft)i) qSiX^ + ^%} 



(pn(Oi + T^fi^^y 



(pnay^y 



(qnco^y 






1\ r3Sm(<X + ft)| + M^) -|- gB^^i+i 2S,,i(a 4- Ct>i) + 2-^»i+l 



(1^0)1 + qmy"-' 



(pco.y^'- 



>pin (" + Wj) + jB^nj 



(««3) 



W + ,? 



where, if a is positive with respect to the co's, C (•5? f* [ ^^i? (^^) niay be expressed by 
the integral 



2lT 



e2M..i r(l -s). 



L (1 — e-^i^) (1 — 6-^2^) 



(§39.) 



Make now 5 = 6, where e is very small ; expand the various terms of the 
asymptotic equality in powers of e, neglecting those higher than the first, and we 
obtain, if the real part of a is positive. 



pn qn 



1 



{pn + 1) {q7i + 1) "— € log 11 n (a + ^i^i + '^^^^g) 

QHl = /"% = 

p — dZ f _____ iy\ 1 

■zz^^ (1 + clog(-- z)] {I + ye] {I + 2M€7rt} dz 



6' 



'■) 






- e log (pw^n + ^«2«)] 
{npw^fll - elog|J72«J - (•ng«2.)^2[l 



logqnoy.^} 



4" 3^1® (a 4" Wj -f Wg) (1 + e) 



n {p^i + 2'«2)[1 — £ iog (P'^'^i + T'^^i)^ 



«pwj [1 - elogpwwj] - nqoy.^\l - eloggwco^] 



+ (1 + e) [p'^ (1 — e logpnwi) + qn (1 — e log §«((><})] 

+ aS/ (a + wi + W2) [1 - e log {pnoi^ + qnoi^)] - ^S/ (a + w^) [1 - e log^HO) J 

- 38/ (a + wg) [1 - elog qno>;_ 



00 



-p 2 



(~y^ 6 



m 



:= 1 ?i'" m 



>^ni (Cf^ 4- Wj + O^g) + 3B 



2/ T^ 3-^^?«+l 2 



S„, (a + w^) + oB 



m+l 



op-if^i^CC + Ct>2) + 2^^* + 3 



(ptt)i + ^0)3)^^ 



(_p«l) 



m 



k^i) 



m 



This equality will hold for all values of s, a, co^, w^ ^f the integral be replaced by 
^2(6, a|a)j, 0^2)? the logarithms having their principal value with respect to — (w|4-^2)- 

Equate now the absolute terms in this asymptotic equality, and we find, if a is 
positive with respect to the cos, 



e~«^(-— 0)"~i 



27r JL (1 — e-^^') (1 — e^^'^) 

But we have seen (§17) that 



dz +38/ (a + to^ + co^) -~ gS/ (66 + (o^)—Si (« + %)• 



>S/ (a + o>i + C02) -" gS/ (a + 6)1) — gS/ {a + oig) 

2 u 2 



,Si(a) + 1 



OO ju 
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Hence, provided a is positive with respect to the oj's. 



2^1 H ^ 

and for all values of a, co. and co^ 



2ir 



e~az (^^ z)-^i 



L (1 -^ e-''-'^) (1 - e-~^^ 



ch 



% 



Oil J U)^)\ 



§ 49. Secondly, equate terms involving the first power of e 
We obtain the asymptotic equality 



p» qm, 

log IT n (a + ^1^1 4" ^%%) 

ri\ — •?/io = 



■a 



,. ^s, a 1 Wj, «,) 4-«3 (|- + -|)_pg - n2 log n ^S/^' (0) [pa)i + qoj^f 

- w\ S^(3) (0) 1 [(pa>i + qoj.f log (pwi + gwj) ~- {po),Y log (poji) - (qox^Hog qm^] 
+ ?^ (_^9 + q) -n log n [gS/^^ (a + «i + wg) (pw^ + gwg) -• gS/^* (a + co^i^wi 

- n [gS/^^ (a + wi + 0),) (jowj + gwg) log (pw^ + qa.) ~ oSi® (a + «]) 25«i log pcuj 

" gSi®' (a + wg) g«2 log gcug 

- [jS^' (a + wj + (Ma) log (pwj + qcoo) — 28/ (a + ^1) logpwj 

— nS/ (a + o).,) log gwo 

- log n [Si (a + wi + 03) - gS/ (a + w^) -- gS/ (a + Wg)] 



00 



i-r 






m 



cj ^'m+i (a + (Oi + ft) g) 3^'»^-f:i(^'^ + ^1) _ 3^'^H-i (^^^^ Hl^g) 
_ (jp^i 4- ^«3)'- (pa)if (pm^>^ 



valid for all values of 5, a, 0)^^ and cog, provided the logarithms have their principal 
values with respect to the axis of ■— (coj + Wg). 

In order that the labour of writing down cumbrous formulae like the one just 
obtained may be diminished as much as possible, we propose to introduce a symbolic 
notation suggested by Cayley's notation of matrices. "^^ 

If f{z) be any function of z^ we shall represent symbolically 

by Fg [y(^ + ^)]5 the suffix 2 denoting that we are dealing with two parameters. 
Thus the difference equation for double BernouUian numbers (§ 17) is written 

i\ [Su {z + a>)]= - oS, (z) + z^ 
Similarly F2[3S/^^(a + a))pa} log pai] denotes the function 



* Cayley, * Collected Works/ vol. 2. The corresponding theory for multiple gamma functions will be 
developed by employing a symbolic notation ab initio. 
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- 2^1^^^ (^ + oy^) 5^2 log gcog. 

[The analogy with the matrix notation would be more complete if p and (j were 
replaced by pi and jo^. The convention adopted here is, however, quite natural. 
And now our asymptotic equality may be written 

pn qn 

log II 11 {a + m^o)i + ^^2^2) 
W] = W2 — 



■a 



3^ 4 {s, « I ^1? ^^2) 



s = 



+ P3 [^^ log n •— ^" (x + |-) ] + (p + q) [n log ?i — n 



+ F3 [2S1' (« + «) log {poi) ] 4- log w [1 - Si ('') ] 



»n 



1 mti 



{pay) 



for, as may be readily proved, 

When the variable a is positive with respect to the oi's, we note that the part of 
the absolute term in this asymptotic equality which is equal to 



"3 

3^ 



b2 v*^? ^ 



(x>l, OJo) 



s = 



may be written 






IT 



L 



g-ff^ (.— ;?;^- 1 |fy 4. log(-- Z) } d0 
(1 — 6 - *^1") (1 — ^ - ^'^a^) 



2M7n2S/ (a) 



which is the expression which has been proved equal to ^Sq (a) 2 (m + m') ttc 



log 



To (a) 



P% («i. wg) 



by the process of differentiation under the sign of contour 



integration. 



§ 50. But if we take the expression for log Fg (a) which has been obtained in § 24, 



log r^ (a) - - 731 (wi, W2) + ay22 (wj, cog) + log a 



00 00 



and write it 



m^ = W2 = 



log (a + O) — log n 



a 



log Fg (a) = ,^ y^i (a>i , Ci)^) + ay^^ {oy^, w^) + log a 



2 



+ Li S 2' 



log {a -^ D) — log 11 



12 ^ 2a2 



+ 202 



we may obtain this expressiou independently. 
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For putting a = in the asymptotic expansion of 



pn qn 

log 11 11 (a + a) 

mi = m.) ~ 



we 



find 



pn qn 



log IT 11 {"in^o)^ + ^^2/^2) 

niy = m2 ~ 



^ 



"3 

as 



{3 (s, a I o}y, 6J2) + log a + pq [n^ log » — n^ (|- + -g-) 



-f (p + <?')["' log n — ■^^] + ^*^F3[3S/^* (ft))pwlog|>&>] + ?^F2 [gS^®' (&>)_2:)a> log pco 
+ Fa [3S/ (c) log^^c] + log n\l~- 38/ (0) ] 



00 ( V'^~^ 



m ~ 1 



mn 






\ • • 



t a s 



o e e a e a 8 



(1). 



This is the extension of Stikling's Theorem to two parameters. If for all values of 
o)^ and ^2, ^® P^^ 



logp2K.^2) =^ — L^ 



a = 



a 



s 



l{s, 



a 



ojj^5 ct)^) + log a 



we may call p^ (coj, ojg) the double Stirling function of the parameters o)^ and a)^. 
It is the same as the third double gamma modular form previously defined. For we 
have by § 43, 



logp2(oJi, W,) 



a ^-- 



6 r 6-- (-S)-Hlog (-2) + 7]'fe 1 ^,_, ,; 



TT 



s = 



(1 -» ^-"'^'^) (1 — e-''^') 



9 



TfigS/ (0 I Oi^, Oig) 



a 



^^15 %) 



log a by § 42, 



We now see the exact analogy between the function p^ (co^, a)^) and the simple 

Stirling function p^ (oj) ™ Y^(27r/co) as defined in § 31 of the ''Theory of the Gamma 
Function/' 

For a brief inspection shows us that the result of § 30 of that theory may be 
written 



pn 

log 11 (m^ co) r=i _p {■n log n — n) + ^^[S/^' ((i>) paylogpa) 



■/)h — I 



+ [1 + S'l (0)] log n + log p^ (oj) + S/ (o>) log po) 

which is the complete form of Stieling's Theorem for a single parameter 
The analogy between this asymptotic expansion and 
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log n n {m^co^ + moco^) = pq [n^ log n — n^ ( f + i)] + ^^^'F^ [2Si^^^(w) po) log pw 

+ (P + 'Z) ['^ ^og ^' - ^^] + «Fo [3Si<''>(«) po) logpo] 

+ [i - gS'i (o)] log II + log p3 («i, 6)3) + Fa [gS'i (w) log jjw^ 

-|~- 2v ■ ^ iJ o S — — -------- — ■ . V. 

is so evident as to determine the nomenclature. 

Note in the second place that the fundamental asymptotic expansion (A) of 
may, with the symbolic notation subsequently introduced, be written 



38 



■pil qn 1 

t t 7 - 



^.3(5, a|cOi,COo) + 



^ Tp 






1 -T-I foS/^^ C^'. -I- ft) 

— - — - \j\^ -- -^-^ — - 

(s — 1) . n' ^ [ (^jwy-i 
(—■)"' fm + s — 1 



1 



+ '■"■" Jj 



oS'5 (^(^ + oj)" 






t?l 



^ 1 (jm)"''^' 



§ 51. If in this equality, true for all values of 6*, a, co^ and w^, we make 6' 
obtain, since 



~ 1, we 



pn qn 

t t 



1 






S/2)(a)- 



4(s, a] 0)5,0)3) + ^y_— 






+ «F3{2S/^>((«)2>a) logi^ojj- - 3S/->(a)log n + F2{2Sj<->{a + o>) \ogpo>\ 



+ S -^^:rr F^ 



00 



J?l == 1 



7V 



ni+1 



o&;«(r6 -f- CO) 4- 2.^">»i!41 



which is equivalent to the result of diflerentiating the asymptotic expansion of § 49. 
If now we make a = 0, we obtain 



pn qn 



1 



\jih 



s= 1 



C(s, «l«i, «s)4- 



,Si(2) (a) 1 



1 



ffi 



+ ^^F,{,Si^^>(co)pa;log po)} -~- ^S/'H^^) l^g ^^ + F3{,S/'^^(a;) log po>} 






^ . 



J 



We now see that 



T I 

Ijo 



n = 00 m 



1 



'i' I' z;-^"f--^+ 2S/^'(o)logn.™Fl,S/'^>(c.)logpa,} 



— nF^{^Si^^^{o))pwlog po)] 



Jljo 

a ~ 
= 






which is a quantity independent of p and q. 
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But Avhen p :r:r. q = 1, the expression on the left-hand side may be written 



Lt 



n •— Qo 



n n 1 
0^^ 



} 



o— -' log '^ + 9:,7r 1 iog(«i + «2) - log «i -- log wj 

/ 1 1 \ -y^ + 1 ^l + 1 
- (7^ + 1) ( ' + " ) lc>g (^1 + ^2) + ' l^g ^2 + l^^S ^1 



\^i 



CO 



3/ 



0). 



ft)o 



Since the principal values of the logarithms with res]3ect to the axis of — (a)| + <>>3) 
are in all cases to be taken, we see from § 23 that this expression is equal to 

For denoting by a capital letter the logarithm, which has its principal value with 
respect to the axis of — 1, we have 

log (coi + OJ,) :^ Log (o)^ + ai^), 
log co^ === .Log 0)^ + 2mTn, 



log a>2 = Log 0)2 + ^niTTL 



x4nd therefore 



L^ 

a=:0 
s = 1 






1; 

a 



Lt 



n ^ CO 



~ n n 1 

o 

'^^ 



.^^ log n + "^;";- log («i + (^,) - log w^ 



log 0)2 



J 



{n + 1) ( " + J log (^1 + ^^^2) + ^ ^ log <^2 + --77"- log o>i 



ft) 



COo 



COi 



CO: 






11/ -\~ 1 111 "V' 1 

-^"— — 2m7r6 + ~ """ 2in'7Tt 



COi 



COo 



ft)] -f ft)o 

2ft)|rjL)o 



2('??i + m)Tnri 



721 (o>i, 0J2) + 38/'^ (0) 2(m + m')7r6. 



We notice that this formula agrees with our previous results. For by § 47 



721(^1. ^2) + 2S\(o)[(m + m')27n] 



Lt 



i C er^''^ [log {■— 76) --V ^]dz 
2it] L ( 1 — er''^) (1 — e-'^'^') 



1 



a 



2M7rtoSi(o), 



and by § 42 this last expression 



™ L^ 

s = 1 



4('^', a] 0)1, 0)2) 4 



,S/2)(f^) 1- 



s — 1 



a 



It is worth noticing that incidentally 702 (^1? ^2) ^-^^^ been obtained as a limit in the 



more general form 
722(^1. ^2) = L^ 

n — 00 



' 2^n qn 1 
0^ 



ft)i 4" ft)o T 

"^ log ?1 
ZcOiCOo 

Ct)i -"T C^o ri / J \ 1 

+ -o 1log(pa)i + 50)2)-- logpo)i 



2a)ift)3 



log (^0)2} 
[over] 
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+ ^ log ( jJwj + go),,) + __ ^ (log go), + 2mm) 



0):y 



(jdx 






0)o 

qn 4- 1 



w, 



(log (^cog + 2m7n) 



where log {loo>i + g^o)^) has its principal value with respect to the axis of — (co^ -^-o).^. 
§ 52. Let us next put 6' = 2 in the final equality of § 50. 
We shall obtain, in the limit when n is infinite. 



nil 






Jut 



4(^, a\o)i, o).i) 



,Si(«) (a)- 



+ (1 + log n).2^fio) 



the logarithms having their principal values with respect to the axis of — (co^ + coo). 



rn 



Thus 



L^ 



4(5, a\o)j^, C02) — 



,8^(3) (0) 



9 



1 



a^ 






1 



1 



t t' 

('%«! + %«2)^ 



1 



[log n - log (pwi + r/cojj) + log pwi + log gwg]. 



Now the left-hand side of this relation is independent of p and q, and therefore we 
see, by § 22, on putting p ■= q z= \^ that 



Li 



ct= 



4(6^ a 



0>|, O) 



2/ "" 



,Si(3) (0) 



1 



a' 



+ 



1 



«1«3 



= - y3i('^i> wj) - 2{m + m') 77128^^(0). 



This formula again agrees with results which can be deduced from the integral 
formulae. Incidentally 73i(ft>i, Wj) has been expressed as the more general limit 



y2i(«^i.<«3) - 



loff n 






1 



ft)ift)o 



V V -..-J - ^ 

[log (pa)i + qo).^ --. log pcoi - log qco^^ 



where log {2^(0 1 + qco^) has its principal value with respect to the axis of — (co^ + ^2). 
§ 53. We can now finish the investigation indicated at the beginning of § 50, and 

r (a) 
obtain the expression for log --^^^^ - without integrating under the sign of contour 

integration. We have 



log To {a I Wi, o),) = - 721 + ctyo^ + log a 



2 



where O == "^^^i^i 4" ^^^;2<^i- 

VOL. OXCVl, — A. 



+ Lt t t 

H = oo 



Zi A. 



log (a + O) 



lo^ O 
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Utilising the results of §§ 49-52, we find, when a is positive with respect to the w's, 

-loff To ia) = 1^77- t"" 7i wfi" :r <^'« " 2M7noSi (a) 

+ TiFJ [2S/->(a + w) - 2S/^'(a>)]pwlogjjw} + F3{2So(a + «)logp«} 

+ log n {38/ (0) - 38/ (a)} - log /)3 (o^, og) 

+ a [- 38/''^ (0) 2 (to + «^') TTi + o^f> (0) log *i - 'F.{J^^' (w) logpo)} 

— nFgfaS*'^ (a))pa) log |)co} ] 






4- 2 [-,S,» (0) 2 (m + TO') TT. -f 281* (0) log n^ - F3{3S/^> (a.) logpo.} ] 
tlie logarithms having their principal values with respect to — {oj^ + ^2)- 



> 



? 






Now 38/ (a) - 38/ (0) + a3S/^> (0) + ;^ 38/^) (0), 

for the higher terms in Taylor's expansion all vanish. Hence 



1 To (a) ^ f 



dz 



p, (0)1 , o).) 277 J L (1 - ^"-^■^) (1 - e-'^^) 

+ ,2^0 (a) 2 (M + m + m') ttl + ^S/ (0) 2M7rt , 

which is the expression which has been previously obtained in § 45. 

§ 54. We next proceed to consider the values of the double Riemann C function, 
{o(5j a 1 o)|, ^3)5 when s is a negative integer. 

By the definition of § 39, 

a being positive with respect to the co's. 
Now in § 15 it has been proved that 





0(1 






(1 — 6-^^ 


— (?" 


-W2A 


Therefore 








i% (^s', a 


<^i? 


^2) 


ri:: 



e s a 






tr (1 - s) 






+'-5^.... +(-).-?--;f ..+ ...} 



the latter integral being taken along a small circle surrounding the origin. [Since s 
is an integer, the two line integrals in the standard reduction of § 46 destroy one 
another.] 
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Hence if ^" =: — m, when m Is a positive integer greater than or equal to zero. 

{0(5, a 1 0)1, ft)o) = r (1 + ^^^) 7-'-~4-r- . 
Thus, when 5 is a negative integer, 

4 [s, a wj, 0)^) = ^fzj' ■ 



But by §12 corollary, 



and therefore finally, when 5 is a negative integer, 

4 (s, a I ^1, 0)^) = gS., (a) + ^Bi^, (c^^ , ojo). 

When a is not positive with respect to the oj'b this formula continues to hold, as is 
immediately evident hy the theory of the function 4C"? <^K^i? 0)2) > '^^ which we shall 
shortly proceed. 

§ 55. We proceed next to find the values of 4G^? ^<^ j ^u o)..) for jjositive integral 
(including zero) values of .9. 

We have seen, in § 48, that when a is positive with respect to the oj's, 



27r, 



(L) 



(1 ■- {?~'^i^) (1 -- ^-'^a^) 



2 



S\ (a), 



so that 



4(0, a I CO J, 0)3) = 2S\(a). 



Differentiate with regard to a, and we find 



27r, 



(1j) 



e-'^'d^ 



(l-~ e~^'l^)(l~- (?-'*>2^) 



,SV^>{a). 



Now by § 39, when e is a small quantity, 



4(1 



e,a 



0)l^ 0):yj 



ir(e) 



e 



■SMcTTt ! 



^-«^(— ;a;) 



— e 



27r 



Jl 



^' -' 27r ■' JL (1— e-"i^)(l-e-"2=^) 

e-"-^{l- e log (- a) + . . . }{1- ye + . . . }{1- 2M67r t . . .} 

e(l— c-''i-)(l— e-"2^) 



(?: 



Z 



ch 



SO that, neglecting powers of e above the first, 



e4(l-— €, c/.|a)i, (Oo) == 2S/^^(a){l — ye — 2Me7rt} 



f 



log( — ;^)6^"«-J,^ 



But we immediately deduce from § 53 that 



27r JL (1— e-'^i^)(l — e-'^a^) * 



log r^ (a) = ^\^ ^z^.^f^[zr-^ - \ M^ + «^ + ^^ ) + 7 } A''' («)■ 



da 



3-\r O 



o 
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Thus 



Lt 

s = 1 



Us, a 



^15 ^^/ 



.) + 



3S,«(«y 



s 



1 



VA(^0 + 2(w + m')m5S/^^(a). 



Again difFerentiating, we have 



e-^'^— z)ch 



27r . L (1 — e"'^^') (1 — e"""^^) 







Sj®(a). 



Also, if e be small, we have 



4(2 



e, a 



W], Wo) 






27r 



L (1 ~™ c"~''^i'^)(l- ^.-"2^) 



and r(€ — 1) 



Thus 



r(6) 

€-1 



(l-~ 7^ . . . ) (1+ 6 + . . . ) 



\ 



1 +(i-r)€ . 



9 9 I « 



/ 



4(2 — e, a|wj, 0)2) 



1 7^1 + 2M7ri 



27rJL (1._fl--.^)(i_fr-'«) " 



z. 



so that, since 

t/;,<^>(al«„ CO,) = 2^fL7^)(Jz^!|' + [2(M + m + 7n')^. + y]2S/«)(a) , 
we have 






4(5, a I oj], oj^) 



(s -- 2) a)ia)3_ 



^/^^(^O 



1 



ft)Tft)o 



2(m + m')7n2Sj<^')(a). 



Tabulating our results, we see that 



4 {s, a 



«1> ^3) = 


- 5-*--i,, *'•'("). 


when ,9 > 2 ; 




- f ""^ + «'(«) - 


— 2 (^n + 




- sm") ., „,,„ 


J J_ / 'iio »JL /no » /^ 



= 2 



s-- 1 






5 = ] 



.9=0; 
s < 0, 



These formulae hold for all values of the variable a, though we have only established 
them for the case when a is positive with respect to the co's. They evidently agree 
with the results established for the case c^i = w^ in ''The Theory of the G Func- 
tion/' § 34. 

§ 56. In a note appended to the '' Theory of the Gamma Function," it was stated 
that a theory of the simple Riemann ^ function had been developed by Mellin.'^' It 



^ Mellin, ' Acta Societatis Fennicae/ vol. 24, No. 10, 1899. 
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wa8, however, published after my paper had been sent to press, and I was therefore 
ignorant of his elegant results. Expressed in the notation which I have adopted, his 
method is as follows. 

It is evident from the expression of the function ^(5, a, co) as an asymptotic limit, 
that its importance lies in the fact that it is a solution of the difference equation 

1 

when s has any value, real or complex. [From this result we see at once that w^e 
should expect that, when s is a positive integer, the simple ^ function should be sub- 
stantially a derivative of the gamma function, and when 5 is a negative integer, a 
BernouUian function.] Now when '^(5) > 1, the simplest solution of our difference 
equation is evidently 

This solution becomes nugatory when ^(5) j> 1, but Melltn has succeeded in finding 
a modified solution by the following ingenious modification of Mittag-Leffleb/s 
process. 

We construct the function 

when 1 > %{s) < — h, by writing »— s in place of m in the m}^^ simple BernouUian 
function 

and taking the sum of the first h -^ 1 terms, k being of course a positive integer. 

Thus c;n / ! X ^f^"' C^~' .^'^V"~"A -.-r-p / \ 

And. now ^(^, co^ co) is defined by the relation 

C {^> ^-. w) 

= •— S^^;, {a\(o) + t I — --^ — -7„ — S^^,;^. {a + (n-jr l)o)\o)} + B^sjc {ci' + no) | co) 

We readily see that this function formally satisfies our fundamental difference 
equation, and we may at once prove that the series does define a function existing 
over the whole plane. 

For when s = 0, — 1, . . ., — • k, 

it is evident that B__s^k {(^\o)) = S^^ {a | w) -f- constant^ 

and therefore B^-s,ji:{c('\o)) — • B^s,k{^) •— a"-^ = 0. 
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When a is large and s has any value, the left-hand side of this relation may be 
expanded in the form . 

Po (s) P, (s) P,. (s) 

where the P's are integral polynomials of 5 of degree indicated by their suffices. 
As they vanish for (^ + 1) values of s, we mnst have 

S ; /// 4- ^^ — S J (a) -- cr' — -^^±iii^ ■+ ?-^-±^^^ 4- 

CO "i 

Thns i(s, a, oj) is convers^ent with t ; — — — ^^^ ^^j^^ jg therefore convergent, 

provided 

%{s + k+l)> 1, or f . (s) > -k 

which is the condition with which we started. 

There is one point which does not arise in the work of Meltjn, wlio takes the case 
0) =z I, It is that throughout we must w^ofk with many-valued functions with s as 
index, which have their principal values with respect to the axis of — o). For In 
expanding 8„s,k{(^ + ^) "^ S^sjc{^) — <^i'^% where a is large, we have tacitly 

assumed that log (a -^ co) •=: log a + log ( 1 H — ) ? which, unless m is real, is not the 

case when a is large and nearly real and negative, so that a and a + (^ he on oppo- 
site sides of the axis of -— 1 , if this axis is the axis of the logarithms. 

§ 57. It is now possible to construct the double Riemann ^ function by extending 
the previous analysis. The function so constructed might be made fundamental in 
the theory of double gamma functions and double Bernoullian numbers, these 
functions arising for particular values of the variable s. We will indicate the 
development of the theory from this point of view, for brevity establishing only the 
principal results, or those which, as in §§ 54, 55, have been established only over part 
of the a plane. 

The double Riemann t, function ^^{s, a\(Oi, cog) is the simplest solution of the differ- 

ence equation 

1 
f{a + 0)1 + %) -/{^ + «i) -/(a + %) +/(^0 = 7^^ 

a, s, (w^ and (j)^ having any complex values such that o^^/^i is not real and negative. 
The determination of a^"-^ will appear in the course of the investigation. 

In the first place, it is evident that when ^ (s) > 2, a solution is given by the 
sei les 

2 ^ ^ — ^ — ^ 

mi=0 m,:^0 (^^ + Wift>l + ^%^3,y ' 

which will then, by Eisenstein's theorem, be convergent. 
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When 2 > |l (s) > — • h, h being some positive integer, we form a modified solution 
by the introduction of the function ^S^^,/, (a|co^, co^), formed as follows. We take the 
mih double BernouUian function, ^^,i^{a\(o^,o).2J, write — s in place of m, and then 
take the sum of the first (h -{- 2) terms. 

We thus have 



1 W;^ -f Wg 3B3 ^^ .^ /'~'A ^'^>' 



+ 1 



^-sjc {a I coi, co^) - (j^ _ ^y(2 -^ s) 0)10)3 a-2 2 (1 ~ s) o)ioJ. c^-^ + V + ,^0 \ r ] a^'r ' 



and then C2 (^) ^^^ 1 ^1? ^2) ^^ given by 



CO CO 



— gS^,,^. {a + {m^ + 1) CO^ + ^^r2<^2} 



— gS^,,^ [a + m,^a;i + (m^ + 1) co^} + ^S^,,^ {a + mico^ + m^co^} 
an ex|)ression which for shortness we shall sometimes write 

00 GO 

It is at once evident that the function so defined formally satisfies the fundamental 
difference equation, and we may readily prove that the series is in general con- 
vergent. 

For when s = 0, — 1, — 2, . . ., — A:, — ■ {k + 1), obviously 

where X and /x are constants. 
And therefore 

x(^\s,k) vanishes when 6' = 0, •— 1, — 2, .,.,—/(*, — (^' + 1). 
When z is large this expression admits of expansion in the form 



V 



where Pq (s), P^ (s), . . . are integral polynomials in s of degree indicated by their 
suffixes, provided that the logarithms which intervene in defining the many-valued 
functions with s as index are such that, when e is small compared with z, 

log ^ + log ri + -J = log {z + e). 

If the axes of coi and o)^ include the axis of — 1, this will not be the case for 
terms of the double series, for which the numbers m^ and nic^ in the term 
z =: a -}- m^coy + m^oJo are large, unless the logarithms have their principal value 
with respect to some line between the axes of — co^ and — • w^. We take this line 
to be — (coj + coo). 
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And now, since x (^ | ^\ ^) vanishes for h -\- 2 values of s, we see that its ex|)ansion 



when z is large must be 






c^s+k+2 







The series for 4('^, ajcoj, co^) is therefore convergent with 



00' 00 



1 



It is then convergent when 5 + ^ + ^ > 2 or 11 (6*) > — k, \ s | being finite. 
We have then obtained a solution of the difference equation 



f[a -f- ^1 + ^2) — f{^^ "1" ^1) "^ f{f^ + ^2) "^ fy^) 



a' 



where a"^ has its principal value with respect to the axis of •— (o>j + ^2)' which is 
valid for all values of s^ a, m^ and w^. 

§ 58. The identity of the function 4(^5^|^iJ ^2) j^^^ defined with that previously 
employed is easily seen. 

From the mode of formation of c^^__sjc{<^^\^{^ (^=1) it is evident that when a is positive 
with respect to the co's, we have 

. 1 
c>^-sjc {(^ 1 o)i, ^2) ~ "th^ sum of the first (k + 2) terms In the expansion in powers of - of 



iV (1 - 

27r 



s) 



e 



2M.7ns 



and therefore 



f 



L (1 -^ e-co,.^)^ _ ^^co,.) "' 



2S«-^,a5 (<^l<^iJ ^2) 



^^ (^ ^) ^2M?rt6' 



g-«.*; (^ _ ^^y-1 



^ 



TT 



L (1 — e"-""^) (1 — e-''^) 



i^ * 



Therefore when n is large we have the asymptotic expansion 



2S-...00 (^ + {jpn + l-)t^i 1 ^1. ^2! = 



1 



1 



2a + ft), + ft);, 



(5 — 1) (s — 2) ft)3_ft). (Fift)i) 






s — 1 ' 2ft)j^ft)n (p^ooi)^"^ 






1 



S — 1 , (Oi (p?^^i) 






'HI— 1 



•yyi i 



a fornmla which may be proved to be true for all values of a by a term-by~terni 

expansion of the series for nS_,.^oo{c^ + (j^^<^ + O^^il^u %}• 

Now from the expression for t^{s,a\oi^, co.^ given in § 57, we see on taking the 
{pn + 1) first values of m^ and the {qn + 1) first values of m2, that 



^o {^'5 c^ 1 (^1(^2) ^^ ^^^' 



n = <x> 



pn q',i 






c^S^s,k{(^ + (p^'^ + l-)<^i + (5'^^+ i)<^^2^f 
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Putting ^ = 00 , and employing the asymjDtotic expansion just obtained for 
S-s^oo {<^<^ + (p^^ +1} <^i] we obtain ^^ (s, a | w^^, co^) from the same asymptotic equahty 
as that by which it has previously been defined in §§39 and 41. 

§ 59. As an example of the way in which we should proceed in a theory based 
on the double gamma function as defined in § 57, we will prove the relation 






C^is, a\o)^, 0)^) + 



s - 1 



.So' (a) 



= *— xfj^ {a\o)i^ 0)^) + 2('m + m')7ri3SQ (a) 



established in § 55 for the case when a is positive with respect to the co's. 
In the first place, when s =^ 1 -\- e and e is small, we see that 



oS_,^;,(a| ft),, OJ3) 



1 CO^ -\~ COo ^ 

- 4" log <^ j^^oi^) ~" """ + higher powers of e. 



Therefore taking ^g (6', a \ o)^, 0)3) as the limit, when n is infinite, of the sum obtained 
by taking the first (pn + 1 , gn + 1 ) terms of the double series 



^2(1 + e,a) + -~Si^{a) 



n — 00 






1 



a 



^m, = ^2=0 <^f^ + '?/^ift>i + %ft)3 CO^Wg 



-- aSoT^^ + (p^ + 1)^1 + (?^^ + 1)^2]%^^ + (F^ + 1)^1 + (^i^^ + 1)^2] 
+ a^oT^^ + (?^^ + 1) ^2] log [^ + (^/^^ +1)^2] 



the logarithms having their principal values with respect to the axis of — • (co^ + ^2)- 

- q =1, we find 



Putting a =: Oi, p 






1 1 



«. 



'ii 



1 



J- J. f^ t^ 

7b -f 1 



- log^. — {n + 1) ( ;;- + 77) log [coi + ^2) 



^/> 



COiCO.2 



\0)j 



^0. 



r 



0)n 



log (Oy + '' ;_ log «3 + ^:;^_;^' {log (Wi + Wo) - log W; - log OJ.2} 



co^ 



«1«3 



722(^1. ^2) 5 27r6(m + m ). 



ft)lO)o 



And now, in the limit where n is infinite, we find 
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4(5, a\o}i, mq) + 



1 



1^ 



n 



1 



So' (a) 






^1^3 



9 



^1^3 



log n{o}^ + wg) 






1 ^i_^^ 

f'' 



0)^0)2 



a)n + a)ol T fa + (?^ + l)a), o). + (Wcl i 



a + o 



1 a 

n ^ D^ 



+ ayo] + 722 + 2(m + ^^^)7rt . 



2a — a)i — 0)0 



9 



o). a)c) 



XJJ^^^'^ (a I CD;l> <^2) + 2^0^ (^^) 2 (^'^ 4" ^^') ?^^ 



where O ™ m-^o)^ + ???2^^-2? 



the complete form of the result established in § 55 for the case where a is positive 
Avith respect to the co^s. 

We may establish the other results of that paragraph in a similar manner. 

§ 60. Finally we will briefly consider the reduction of the second form of the double 
C function to the double gamma function in the case when s ::rz 0. 

If we put .9 " e, where [e] is very small, we obtain 



a^ j e 



20)^C0:^ 



>{1 — e log a + • 



a (ft), + tOg) 



1^ X "f- 6 -|~ . , . ] { 1 — € log a + » . . } 



+ ol^j (*^i? ^2) { 1 — € log Cl^ + ■> ' * } 



2S/(a)[l — eloga] + e 






4 



a 



fjy ^ Q) 



&)]^a)o 



9 



(Of^O).)^ 



+ higher powers of e. 



And therefore, by the second definition of the double { function 



^.(e, €^j w^, 0)^ =z 2S/(a)[l — elogaj + 



4.0)1 «n 

J. .^ 



a 



0)3 + O)^ 
2o)iO)o 



CO 



S S 2Si'(^^ + ^ + «^i + <^2)(-^ '"^ elog (a + ^ 4" ^1 + ^2) 

'iilj^ — O 7/22 = [_ 



9, 



2 



S|' (a + O + <^2){1 "" ^l^S (^^ + ^ +^2) ' • • i 



+ oS/ (a + ^)[1 — ^ log {a + O) . 



) 

) 






1 + ^ log (a + O) + IT 



wtiere .(2 = miojj + '^^2^2 



Thus 



^2(0, a \ 0)^,(0.) ^ oS/(^0 



and 






to (e, a I o)p 0)3) 



„s/ («)■ 



'/;'/?. (pi 



S S log (a + '?>^'i,^i + ^^^3^2) + 



/)lj^=:0 //?o-— 
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Ct C0| -f- Ce);) 






4cL)jft)3 2 






1 (pn + 1) {qn + 1) 



2 



in the limit wh6)n '/^ is made infinite. 

Substituting the value of ^{{ct) in the various terms, expanding the logarithms, and 
re-arranging the result in powers of n, we find, with the symbolic notation of § 49, 






^cj (s, a I ft)i, cDg) •— gSj ^(a) 






log n n (a + O) + P9\j^ Jlog ^'^ — ^^^(i + i)] + (P + <][) \j^ log ^^ — n 







O 



+ --F2[2Si^^''(o))(2JC0)^logJ9ft)] + '^'^"^%{,^i^K'^ + ^)i^<^logj)ct)" 
+ F3[3B/(a + '^)^c)gpw] + [1 — oSi (<^0] l^g '^^^ 

in the limit when n is made infinite. 

Since the left-hand side of this equality is by the definition of i^{s, a]^^, ca^) fin: 
unless a is at one of the points 



(wicoi + m^^ 



m 



mcy 



0. 1 



9 



• « « 



CO 



*Jm X, ^, • • tj ^■•^ 



we see that w^e have thus been led in a purely algebraical manner to the determina- 
tion of the dominant terms of the fundamental expansion of § 49. 
If we make a = 0, and remember the definition of § 50, viz. :— 



l^g P^ i^v %) == T^^ 



s ^-- 



^3 (s, a 1 03p 0)3) - gS/Ca) 



t*? 



+ log a 



we arrive at the dominant teiins of the extension of Stirling's theorem to two 
par am e u e r s . 

If we utilise this result in conjunction with the one just obtained, we find 



Lt 

s-^ 



■^sM -2Si'(a)" 



s 



log P3(«l. ^2) 



ni] 





log a + log n n' O — log II 11^ (a + O) + ~~ — F^ [pw log j9oi] 



^1% 



+ F3 [{2S0 {a+ m)- 3S0 (w)} log^ft)] - sSq (a) log n, 

and by the definition of the double gamma function the expression last written 
reduces to 

2 Y 2 
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log r^ (a) - 



«- 



C0^ (Ocy 



— a 



fjOi + COc, 



9, 



.1 rV 



2(/ni + m) TTL 



m 



Thus 



exp 



. i 



Lt 



& ('%(^') - 38/ (a) 



s 



.^ ^'3 ^^Ji^l'_.^3) ^-^380 (ft) 2 {ui -V m') ni 



We may utilise this formula to establish the fundamental difference equations for 



the double gamma function. 



By the definition of the double ^ function of § 57 



4 (a + wi) — 4 (a) = Ijt 



n = T> 



ob_^5^ CC- -j" o)^ -f" y^ ■^" 1 ^■">) "^ ^i^^s^/c ^^ "4" ('^^ "4" 1) 



0) 



9: 



V, 



1 



,n,^o(ct + m^wgy 



Therefore, in the limit when 9i is infinite, 



Lt 

H = 



'?3 («' + ^1) - 3^^/ (^^ -^ «l) ^3 (^0 - A' («')' 



.9 S 

= S log (a + "^^^^2) — 2^/ [s^ + 0)^ + (n + 1) <^^-3] ^^^8' I ^^ + ^1 + (^^ + 1) ^^p 

'I1I2 :=■ 

+ 2,8/ [c(' + ('/I + 1) W9,] log [a + (n + 1) co.,] + ^ '^ + .7 [^^' "1™ (^^' ^" ^-^ 



COj + (0^ 



2ft)o 



On reduction we obtain 



log 



1 3 (f/^ + (Oi) ^^_ .3 ^,^.^ ^ .^^^/^ ^^ ^c5^/(^^, I ^^) 



A («) 



— J_jf 



n ™ CO 



S log [a + ^^2^2) 

7n., = 



f/^ + (/// + 1) f'^a 



Ot); 



|r j log nC02 + il' 



This latter exj)ression is, by the expansion, obtained in the ^^ Theory of the Gamma 
Function," § 30, equal to 

— Jog "^^y^-^' ^ 2m m b^ {a o)^). 
Pi (^3) 

[The term 2m in {71 +1) which arises is absorbed by the identities which change 
log {a + ^n^oig) i^to log ^>^2^^^ + log ( 1 + — - ) ^ The prescription of the absolute 

\ '^%^3/ 

logarithms has been throughout left indeterminate/J 
We thus have 

one of the fundamental formulae of 5 23. 

Sufficient indication has perhaps now been given of the alternative development of 
t]ie theory of the double ^ function. 
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o 



Note, — The asymptotic expansions of this part were obtained in my original cast of 
this theory, to which reference has already been made in the note which follows the 
" Theory of the Gamma Function," by the assumption that they would involve merely 
powers of n and log n coupled with inductive processes. Such a method, though 
long, is, could the fundamental assumption be justified, probably the most elementary 
way of obtaining these results. 

Additional Note added July 5, 1900.— ^Br. Hobson has pointed out that, if we 
admit the validity of the application of the calculus of operators to a parameter in 
the subject of integration of a contour integral, the theory from | 57 onwards may be 
developed in a very eleo:ant manner. 

J. Jo 

We take the formula 



f{a + o)^ + o>o) —f{a + o>i) --fa + ox^ + f{a) = 



a/ 



and the known theorem 



1 ^r (1 - s) , _ 



a/ 



2 



IT 



and deduce 



f\a) = 



^''(l - s) 



1 



2lT 



d 



{e 'da — l)(e ''da 

J (r- 






zf-'d^ 



z 



2tlT 



(1 — e-""^) (1 — e-^^) 



; dz. 



Part IV. 

The Multiplication, Transformation, and Integral Formulm for the 

Double Gamma Function, 

§ 61. After the developments of Part III., we now return to the pure theory of 
the double gamma function. As regards the multiplication and transformation 
theories, two distinct courses are open to us. We may either proceed entirely alge- 
braically, making use of the limit theorems which have been established, and so deduce 
the required results without the intervention of contour integrals at any stage, or we 
may directly utilise these latter to obtain the formulae in question. The former 
course is, on abstract grounds, preferable r^ we ought to deduce algebraical results 
by algebraical processes. But it is open to the fatal objection of leading to very 
lengthy algebra. We will employ the two methods, side by side, to deduce the 
multiplication formulae, and it will be observed that the second method is both more 
elegant and more speedy. For the sake of brevity, the results of the transformation 

* In the first sketch of this theory, before the discovery of the contour integral expressions, all the 
results were obtained in this way. 
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theory are obtained solely by this course. Inasmuch as the function Vo(mz\(j}j, coo) 
can be expressed in terms of the function vAz\ '' ? '" ), the multiplication theory can 



be deduced from the theory of transformation. As tlie work of obtaining the new 
expression is in every case almost equal to that of obtaining the results ah initio.^ we 
adopt the latter course. 



MultijjUcMtion Theory. 



§ 62. We have from the definition (§ 19 



00 w I 

^b:j^ imz) = ^2t t 7 - ^ '^^ ^ - ',. 



and tlierefore 



CO a:- 






n}?^bf^hnz) = - 2 S t r -— 

\ in m I 

ill — I ii(, — 1 CO .00 - "I 

. « 9 'V "V V '^ 

,■ = „ = ™, = „■,, ., / ,, ^. ;^-«, ^ .s«2 ^ ,^^^^^^^ V 

\ '111 '111 / 

the parameters being understood to be co^ and rx).) wlieii not explicitly written. 
Integrate with respect to z and we obtain 

w^i^f^ {mz) :=^i ' '% ' ^sf^ iz + ^'""-^ ^ -^''^) + r' . , . , , (i.. ), 



where r is constant with respect to 



'-;•' 



Now , (•^) / \ / \ I 'k L V V^ [ -^' ."^- 



^3^'^ (^) = "- 7-21 (^o f^^) + ,a + S^^' ] (i^"+7D^ 



A A'^ 



Substitute from this relation in the identity (l)^ take the same number of terms 

1 1 

involving z on both, sides of the equality, and j-emember that y---^-^—~ ^ ~_r^ is 



always to be regarded as a single entity. We find that, in tlie limit when n is 
infinite, we have 



/ j lib 2L^ 2^ ^-\0 —— /ll/ JLf ^ C\Q 

i^"' *^" 

nil m^ Till m.T, 



where 11 =::= "i^H^i + '^^^sX^^a- 

Now we have seen that (§ 22), in the limit when n m infinite. 
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i>. il 






' 1 



1 



1 



D} 






log 0)j + log COo "— log ((^1 + COo)] 



the principal values of the logarithms being taken. 



Therefore 
so that we have 



ni'- 



1 



^1^3 



og m. 



m^s.f^'^ (mz) 



Til—l m~l 






ri'b 



m^ 






log m , 



Integrate again with respect to z and we find 



m\pQ^^^ {vi z) 






pCDi + ^2^3 



7??,. 



7)ih 



WjCOr, 



log m + -^ 



. (iL), 



where s is constant with respect to z. 

Now we have, in the limit when n is infinite. 



1 



n n 




1 



1 



Hence we find, on equating the irresoluble terms involving z in the same way on 
both sides of the identity (ii.), that, in the limit when n is infinite, 



m - 






L 



723 {^1. %) 







«ll 9».2 



and, therefore, when n is infinite, 



r«t)i+m— 1 1 



72.3 (^1. ^^2) 



Now in § 23 we have seen that 



'iih m.^ 



■n 



1 



721 ""T ^ ^ 03 





7n, ^H. 



o^ 



P '1 



m^ i 



n 



./I 



^ ^'a '^' '^2:2(^1? ^2) 



+ m (w^ — 1) 



^1 -r «2 [n:^'%/ ^• 



9. 



^ -^^ 0?3 + 721 (^1? %) 










n 



Voo 



6>i + ft)o 



9 



0)1 0)0 



log 72. + (^^ + 1) 



1 1 



0) 



«c,, 



— log 0)i 

tt)n 



Ce)n 



log W, 



2m7rt 

0)o 






CO 



1 J 



0>1 + Wo 

20)^0)ry 

x ^ 



{log (a>] + (1)2) — log a>i 



log ft) J, 



and, therefore, after some reduction, we see that 



.s* = m 






.^ log m. 

2ct>ift>o 

.1. iV 
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We thus have 



m i/icj*'^ (mz) 



m—1 m~\ 



t t ^|J,''^z + 



V^l + 2^3 



^) = 5-0 



m 



m^z 



^1^3 



log m + m — ^— ---^ log m. 






Integrate again with respect to z^ determine the constant by making z = 0, 



and remember that L^ 






: we obtain the formula 



9/1 



m Tg (mz) 



n n Fo z + ~-^ — 

2j=o s=^ "_ \ J^_^ 

n n'' r, ^^ — ^ 



- e 



■j — :Sq ('"^-) log ul 



rr = s~0 



m 



the principal value of the logarithm being taken. 

§ 63. By means of the extension of Stirling's theorem to two parameters (§ 50)^ 

it is possible to express the form 11 fF Fg ( ^ -^ j which has thus arisen in terms 

of the double Stirling function p^{o)i, w^)- 
For, in the limit when ?^ is infinite, 



'i\-' (^ 



22 



1^ 'ih 'il 

^ 7-21 + ?'y2-2 -r-r- Y-|-7 

e^ . 2: . n n 







1 + 



Z 



a 



^ o '^ 202 



? 



and therefore, under the same condition, we obtain from the result of the previous 
paragraph 



m' 



vtn + m — I 

m?e^-^'"2~'"""^" n II' 








m?: . mh^ 



] J_ "£' ] ^ O +^~02' 



me 3^« ^'''"■^^ ^^^ 



7/i 



'/;y — 1 ■//?-—! 



11 11' Fr^- •-- ^ 



m ~ 1 rrt — 1 

n 11 

r := .s' = 



, ro). + SCO A y-^^^- 
z + — i---^ - - ) e 



/s + 



rcoj + S(02\2 



w?. 



1 



2 






•\ 11 11 

0^ 



J. ""!"" 



— ± —j^ \ 



IV1I, •!■ 



7 '00 1 •!- .s'co.^ 



'Ill 



o 



^ A 






02 



e 



Make now ;?; = 0, and we find 



'"' "" ^ ''' "/ ^ „-, (r(D^ + §6)3 



n n' F 



= .S :== 



2 



\ 



7)1 



m — 1 rii — 1 

n n 

r = &' = 



i/> 



ft) J -f tS'fo).^ 



•/7. 7? 



7yv. 



11 11' 1 1 + 





m, 'III-: 



771 



X Exp. 



rcoi + ^^^V' 



ra)^ + ^'Wg 



Ql 



., I 



m 





Till m.2 



n 






S t ^ — yo.g) + J ( — :::"^") (^ S ^3 + 721 



?/i 





nil ''^h 

r ' 



O} 



> 







where ( ' '^---^^~-"— 1 1 denotes that in the product the term for which ^ __ q 



m 



simul- 



taneously is to be excluded. 
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We have tlieiij in the limit when n is infinite, 






'inn -\- 'ill — ' I "I 

n n' -™ 

4i/ 

It 1 "-i 

n 11' - 

^-^ 



• m 



wi^^iS ™ 2mhi - wiS + 1 



■in- 









(m — l)(2m — 1) 



9 



il 



12 






(to — 1)^ 

r^ 



Ctlj^O)^ 



(ft)!^ + 0)^^) 



n 1 







Q} 



Utilise now the extension of STIRLIMG^s theorem and the limit formulae for 

We find that 



n 1 n 1 

^ki^ ,ib>i)! ^-^ CllXXvA. ^mi timtf g^ • 

A Ai A ui/ 



loff 



•j-r TT' T1 — 1 I 1 ' ^ I 

r=:0 5 = "" \ '^^^' / 



W%1^ 



2w% 



- m^ 4" 1 } log '^ 

+ (mn -|- m — 1)'^ log {mn -f m — 1) — f (mn + m — 1)^ 

+ ('^ii'^i' + wi — 1)^ Fg, [^Si^'^^ (ft>) ft) log o>] + 2 {mn + ??i — 1) log (mn + m — 1) 

— 2 (Qnn 4- m ~ 1) + (^i^ + ^'^'t — l) Fjj [gS^^^^ (g)) ft) log ft)" 

+ [1 — gS/ (o) ] log {mn 4- m — . i) -|- (i -. 99^2^ ^^^g ^^ ^^^ ^ ^^^ ^^ p^ |-^g^^ ^^^^ j^g ^j 



9 



ffi log tl 



3^3' 



m'^"n/F;> [qS/-^^ (ft)) ft) log ft)] •— 2m^ [n log n -— n.] 



m%i Fe [.3S 1^^^^ (ft)) fti log ft)] 



m'^ 



1 — S{ (0) ] log n — m^Fo [gS/ (ft)) log ft)~ 



m (m 



., X CDi + ft>Q 

- - -^ 



9, 



Bi^^^ (0) 2 (m + ^^') ^^ + '^^P^ [^Si^^^ (ft)) ft) log ft)] -— gSj^^^ (0) log n 



+ F, ISf' (c) log co] 



-X. 



(ft)| ~y" ft)r) / ""I"" 



12 
1 



4 



ft)lft)3 



X 4 ^ 2 (m + ^^')^^ Si^'' (o) + ;;;; log 7i ^ F, [,Si^s>(ft)) log m 

where the logarithms have their principal valnes with respect to the axis of 

' — (ft)| -j-* ^3/' 

The labour of reducing such an expression is evidently very great. It is diminished 
by observing that the result must be independent of w, so that we may neglect all 
terms which involve this letter ; but even then it is only after several steps that we 
prove that the right-hand side is equal to 



[1 — Si {0) ] log m + (1 -- m^) log Pa (wj, wg) 



9 



^ y , ^x I / g X (ft)-, + ft)n)^ , 2m^ — 3m + 1 w,^ + ci)<, , 

2[m -{- m)7n\ — (m^ — m)—'Y'~'^~~^ + "^^"~ """ ^-^^■•---~- -- . ~^~--™ — ^^ 4. 

LA 



m^" 



*Z'in> -f- J~ 



4ft)jft)jj ' 12 ^1^3 * ^ 

1 — gS/ {0) ] log m 4* (1 — m^) log p^ (ft)] , ft)^) + ('^^ "^ 1) 2 (m 4- ^^0 ^^ 2^1 (^) 
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We thus see that 



n n' r^ ' ' 

r =0 H^O 



'Uh 



[fi^{m^,o>^Y-^ . m^^'^^-\ e^^'^'^^^'^-^^'^'^^^K 



§ 64. Ill the case when m -■=: 2 the preceding result has an especial importance. 
It IS convenient to write 



6i>i 



l\U + ^2^ 



r.{^ + | 



<"i> ^3 =y, (z|ojj, Wg) 



OJj^, COg 



72 (2 



W|, cogK 



Tg I z 






(jOj + G)o 

2""' 



6) J, COg j =^ 73 (^' 1 60,^, OJq,)^ 



and in accordance with this notation we put 

To {z I 0)^, (0^) = 70 (^ 



OJ]^ %)• 



Tliese functions evidently correspond to the functions 

(r{z), (r^{z), 0-3(2;), (ro^{z), 

in Weierstrass' theory of elliptic fimctions. 
Omitting the zero argument, we take 



r 



9, 



9. 



^U ^2 



yj (a>j5 coo) and two similar equations. 



that 



'n 'n'r,('^''-t^) 



^0 .? = 



Ti Ts 73^ 



\ 



the parameters coj^ and o)^ being omitted. 

And now, from the result of the preceding paragraph. 



so that 






We thus express the double Stirling function of oj^ and coo in terms of the product 
of double garnma functions whose argument is a half quasi-period. 

We have previously seen in the theory of the simple gamma function that 






and the formula just obtained is the natural extension of this result. 

I 65. From the results of §§ 62 and 68 we see tha,t we may express the multiplica- 
tion formula for V2{z) in the form 
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}ll — 1 Til — 1 



n n r J s + 






T.,{mz) = 



r~Q .s = 



in 



p^'^ - 1 (a)p ct).) 



' ^ ~ g'^r' <jnz) log 9;t -1- (i/*/ -1)2 {m + »fc') irigSj' (o) 



We now proceed to obtain this result at once from the expression of losf -y^^^^ 

as a contour integral. 

We have seen (§45) that, when a is positive with respect to the w's, 



log 



To (cc) 

P2 («]' ^-i) 



= 3S0 (a) (M + ni + m) 2m + ^S/ (0) 2M7rt 



JL 






4„ -1-- I ~: ^^^ — 'iz.^..jl:z^j^ — ZjL-1 — ll ^^ 



and, therefore, m being a positive integer. 



fll — 1 Hi 

n r 

= s = 



"1 

log n 11 



1 J To ( '^■'^ + " ) ^ 



P3 («^P ^3) 



J 



(M + ??^ + m') 27rt S S So ^ + 



'»i — I m — I 






s~0 



m 



+ w^S/(o)2M7r6+ ,^ 



m—l m—l 

e~'^'(- r;^)~i{log(- z) + 7} ^ X e' 



/ 



r=0 s=0 



TT 






But 



and (§ 14) 



"^1 , , (m—l) to I 



(1 — CT"'^ (1 — 6- '^^^") 



1 — er"^^ 



1 — e z 



m 



m~l )/i— 1 



t t S() €1 + 



7'0>{ -j" 6YO0' 



r^O s=:0 



7)1 



So ' ^* 






m^Bi(coi, a)o) + 3B 






= gSo {ma I Wj^, w^) + (1 — m^) gBj (w^, coo). 



Therefore when a is positive with respect to the co's. 



)ll — l Yll~\ 



log n n 



•< 



rJa + 



fCOi -f 50)o 



r=0 s^Q 



m 



(M + ^'^ + 'in)2Tri Soh^^\^\^ ^2) + 2M7rtoS/(o) 



4" (^i + m')27n(l -— m'^) 36^(0;]^, Wo) + , 



L 



OiiZs 



L (}_-.e-'t){l-^i) 



Since m is a positive integer, the axis L defined with reference to the parameters 
^1, 0)^ is the same as that defined with reference to the parameters m^i, mw.^ ^^^ "the 
lines representing these two sets of parameters are coincident. If then we change 
z into mz, the integral last written becomes one which (§45) is equal to 



log 7"-' — , + oS/(ma){logm — 2M7n] 

the arithmetic value of log m being taken. 

2 z 2 



3 



^^{ma) (m + rn')27n, 
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i 



rjLlivl ullv3X i3Xv/x o 



lop' n n 



rm, 4- SO) A " 



.s=0 ^ 



ni 



p,{co„o>,) 



J 



{in + m^)2m(l — tii^) 08/(0) 

+ oS;^'{ma)logm + log 



Fg {ma) 



which is the result required. This result has of course only been proved by means 
of the contour integral, under the assujnption that a is positive with respect to the 
co's. To establish it in general we should appeal to the principle of continuity. 

§ 66. Before concluding the multiplication theory, we deduce expressions for the 
values of 



m—\ TO=1 



% i' xjJ^'H^'''^- ^^'], where ;■ = 1, 2, 3, . . 



We recall from § 29, that within a circle of sufficiently small radius surrounding 






the origin, we have 



logI\(^) 



1 



OQZ 






9 






* % % 



Again, from, the multiplication theorem of the preceding paragraph, we have, by a 
similar expansion, 



log V^^nz) 



{m + w')2trt(l — ^^'^'^')sS/00 + '^ '"'^^)^-^gP2(^i5 %) 



log m 



mrf 



oS/(o) + mz Si^{o) + -^-^" ■^i^{o) + log Tj^ 



Til — 1 Til—X 



+ log n n' Fg 



T^l + g«3 



m~l m—1 



L ^ 



m 















9. 



m 






Combine these two theorems, and equate coefficients of various powers of z in. the 
resulting identity. We find 



72^, 



o>i, m.2) 



m 



1)1 



Sq{o) log m 



1 



m 



"M — 1 Yil '=■ 1 



^ ^SL (1) /£^^L_±^I^ 



rail^i? ^a) 



'??i.-' 



■'/?■?/ 



1 



a' 



^Q^\o) log m 



•w 



1 Dt— 1 1)i— I 

— 1 2^ = g=0 



(-2) /i^^l + 



fE- 



and, when r is greater than 2. 



,,=0 (/=0 \ "'• / ),i| =0 iM.j= 



1 
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Transformation Theory. 

^67. We shall now consider the theory of the transformation of the parameters 
of the double gamma function. It must not be supposed that we intend to consider 
the general linear transformation. 

There exists no such theory for the present functions—at any rate, no theory 

having the simplicity and elegance which is characteristic of the elliptic functions, 

and the reason is obvious — ^the change of a)^ into coj^ + «^:3 makes no difierence of form 

in such a series as 

1 



'/Hi 



^,_oo m,=:--oo (^ + m^coj + m^coc^ 



but it makes a change of comparatively great complexity in such a series as 



CO 00 



ft,. 1=0 m., = 



The former series is the basis of those occurring in the theory of elliptic functions, 
the latter of those occurring in double gamma functions. We shall then limit our 
consideration to transformations which result from the change of o)^ and co^ ^i^to 
o)j/p and ojpjq respectively, p and q being positive integers. 

By definition we have 

CO 00 J_ 



//i,=o v».j=o ('^ + ^ly 



where O =: v/^x^i + w^cu.^. 



Hence i//^ 



(3) 






00 00 



Zj ^ ^ 



/ 



m,=0 m, = I ^ _!.. "^ ^^1 _j_ 



Z + 



P 






r=rO S = 



CO 



1 1 





TO), 



SCOo 



n, m, [ Z + J- 4- Y ■^" ^^H^l +" ^^^3^2/ 



V 



p-1 r/-l 



— 'S ^ i// ^3) f -/ 4- -- -^ 4-. 



S:=cO S = 



9 'ft) J 



5ft), 



S' 



it being understood that the parameters when not expressed are always Wi and ce>.> 
On integrating successively three times with respect to z, we shall find 



lo 



g 1\ (^ z 



ft)n ft)o 



i^ ^ 



X2(^-) + log 11 






2 {7v + 
s = \ 



•rft)| , sft)y 



P 






J. 



where Xs ('^) is an algebraical polynomial in z of order 2. 

As has been stated, Jt is possible to obtain ^s ('^) hy purely algebraical processes, 
use being made of the limit theorems previously established. We may, however, 
obtain its value as follows : — ^ 
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In the relation which has been obtained, change z into z + (f^Jp ^i^d subti'act the 
first result from the^one so formed. We find by § 23, 



r. 



z 



log: 



ft>o\ 



N / 



Pi 



ft). 



\^ 



2mp^(j m Sj 



ft)o 



X2(-«)"X3 2 + 



COi 



.p J 



+ log n \ 



\v,U + 



SCOo 



9 



^3 I L 



'/-I 



2m TTt 1/ b^ ( z + — ==■ 



Wo 



/^ 



ft) 



does 



where m has the value asigned in § 22, and where '?%„^ = 0, unless < 



(jO 



CO:, 



'pH/ 



lie between — 1 and — .— , in which case m 
Now (^^ Theory of the Gamma Function," § 7) 



1 as 



/ ft>o 



i 4^ ""^ does not 
is positive or negative 



log n r 



2: 



mjprn 



Sft);, 



log Tj ( ^ 



ft>r 



ft)., 



^ ^ + g log pi (o),) - log Pi I - 1 , 



and by § 1 8 of the same paper. 



Ti / 



K5i 1 /V 

\ 



ft)c> 



•^ . .-.i / 5'ft),T 



Of).-, 



We therefore have 



Similarly, 



Xs (^ + "^ ) - X2 (^) 

X2 (^ + ^] - X3 (^) 



27TL (fiipyq —- nt) S/ ( z 


ft)., 

/■J 

^2 


27r6 (wip,^ — ni) S/ [2: 


1^ 



where m has the value assigned in § 21, and m'^,,,^ differs from it in that pw^ + Qco^ 
must in the definition be substituted for Mj^ + o).,. 

Now we have seen that (§ 22) m — m = ^ 1, the upper or lower sign being taken 



ft)o 



o. \ ♦ 



as I( ^~ ) is negative or positive ; and, since p and q are positive integers, the same 

is true of m^.q — m'p.rf 
iiius 



"^np, 



'pn 



m 



m;,„ -- ni ^ /x^„ (say)^ 



and now ^^ {z) satisfies the two difference relations 



Vo f Z "4- 



ft)i 



P 



ft)o 



Xs (^) = 2m fXj^.q S/ (2; ^^) 



-j~ / 



Xs ^ T" ^~ ) 



^ 



Xs (^) ^ 27rt ^^,^ S/ ( r<; 



Oi 



i5 
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and, therefore, vsiiice ^2 (^) is a quadratic polynomial in z, we must have 



Ct)j Wc 



Xg {%) = 2776 fip^g 2^/ ( ^ I 7:^5 ^^ ) + constant. 
If we determine the constant by making z = 0, we have finally 



r\ 






2^' q.] 



•-0 .=0 ^\ 2^ qj 



p-i Q -1 

n 

r = s ^0 



11 1 o I "i 



^ '^'^^l^p,q 2^0 ( 






)■ 



3 



\ 2^ 



q 



From the values of w^,^ and m, it is readily seen that /x^,,^ = 0, unless the axes to 
(oij + ^2) and (qcoi + jyco^) include the axis of — 1, in which case 

f^p.q = — 1 if l(a)| + 0)^) is + ve and I{qo)i + j9co^) — i^e 
/^ii?i? = + 1 if l(a)^ + ^^s) is — ve and I(g'o>j^ + p^g) + ^^. 
I 68. The constant which enters into the transformation formula of the preceding 



paragraph can be expressed in terms of ^^(a)^, co^ and p 



(O, O): 



1 



? 



P Q/ 



|. For this pur- 



pose we consider the contour integral which represents the double gamma function. 
Since p is a positive integer, 



roii 



1 -« g-witt P- 1 p~a^^J 



i — e -J, 



r = 



e p 



And therefore if the integral, its contour, and the logarithm which occurs in the 
subject of integration, be defined as in Part III., we have 



I r e - «^ (~ 2:) - 1 {log (- 0) + 7} dz 



'Ztt J L 



Wl?- 



(1 - .-^•) (1 - e-i;) 



2) -1 ^- 1 



2ir 



e-"«-^ ( - .^)-i { log (- ,^) + 7} 2 2 6- 

r = ,s = 



( 



'. -f- 



L 



(1 



***i'^) (1 — (r "^^ 



dz 



for the bisector of the angle between the axes of m^j-p and oy^^/q is the same line as 
the bisector of the angle between the axes of I/cd^ and l/co^. We therefore have by 
§ 45, when a is positive with respect to the co's, 



i. \ (a, 



log "^ 



(Wi w. 



JL^Sl - ,S. a 



(Wi CO 



2*^0 



9. 



P9 \ J 



?^?) ^^^'" "^ ""''^'''^ 2m--oS/ (a 



^1 ^2 \ o 



p s/ 



p — 1 <j — 1 
: S 2 lose < 



'i-? 



Tola + -^, -^ 

P3 (^P ^3) 



> 
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.1 



Now by § 39 M ^0, unless the axes of - and (wj + wo) include the axis of 

which case 

M := ^ 1 as I{a)i + C0.3) is positive or negative. 



1, in 



Therefore M^;,,^= 0, unless the axes of -^ and (qm^ + pa)^) include the axis of 
in which case 



1 



M^,,^ = ^ 1, as I{qo)i + po)^) ^^ positive or negative. 



Again, by § 14, 



.SqI a 



«i 



2J 



(1/ 






0) 



3 \ 



^ 






b| (a>j5 coo) 



D /^l ^• 






We therefore have 



r, ( a 



log 



0)i 



I p 



9/ 









2 2 log 



T(jd 



5C0o\ 



I' 9 



P-i{o^l 



, Wg 



mA(«i> ^2) 



.B, ' "^ 



&)o 



p 2 



(m + '?'^^') 27T6 + 2^^( ^ 



+ o.Sn { a 



M 



I "^3 



i / 



r/ 



m —' m 



f 






M" 



27rt (M^,,,^ - M) 



But from, the values which have just been given, it is clear that 



XyjLqx A fj J^T.L 



l^p^q" 



We thus have, when a is positive with respect to the co's. 



log Fg ( a 



60, Ct)., 



\ i*9~lf/-l 



rct)i . ^(jdA 



P 9 1 



t t log rj a + — ^ + -'^ 1 + 27rt/X;,,,/2.S(} 
'=0 s=o ■ \ P 9 I 



a 



COi 0)0 \ 



? 9 



+ log p. 



\ / J J 






P 9 !. 



{971 + m') 27T'L 






\ 1^ ^ 



This result agrees with that of § 67, and on comparison of the absolute term, we 
see that 



log n n'rJ-J+-;^] 



P*? (log po(<y], Wg) — (m + m') 27ri gBj (w,, Wg)} 



r 



I % P3 I ■,, 



/CO;^ COc 



\l^ 9 



(TO + m' + M^,?) 2^t gB, i^j , - ji 
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This may be regarded as the transformation formula for the double Stirling 
function /02 (^i? ^2)- 



Notice that when p =: q =z m, fi^^,,^ = 0, g^i 
fore the preceding formula becomes 



^1 ^ 



,B^ (w|j 0)2), and there- 



i-5 — -- ^9 



log n II' 1 

r = s = 



m 



m^ log ^3 (to,, 0)3) - log ^0 ^ J , J 



+ (1 — ra^) 2ttl {m + '^') gBi (<^i, o)^). 



Comparing this with the result obtained in § 63, we find 



^^g P2 






log p^ (o),, 0)^) + [1 - Si (o)] log m. 



This result may also be obtained by the transformation of the line integral 
which expresses log ^^(ajj, oj^)- 

§ 69. We have still to consider the transformation of the first and second double 



gamma modular forms 



721 (^1. %) and y.g {co^, w^). 



With this object we write symbolically 



00 a? 1 

% %' y - 






and we write 



S« = t S't//3<"'^ ( '^ + -" , where m > 1 ; 



SO that S;;^ is a form analogous to the modular forms introduced into the theory of 
elliptic functions by Abel. 

By § 29 we know that within a circle of sufficiently small radius there exists the 
expansion 



^? 



00 






CO a4 



log Tg {%) r=r. ^ log Z - Zy^^ - ^yy^^ - ^ S^^^^ + 2 S' J^--, - 



2 



30' 



412^ 



Take now the formula of § 67 



log Tgi ;i; 



— , — p:r: 2 S I02; To ^ + + 



'I 



SCOo 



t s' log r, ^-^ + 



+ 'iin\Lp^q 2^0 '^ 



(Wi &)c 



i^ ^/ 



and expand in powers of %. 
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We evidently obtain 



log % 






iO^ 



z 



Ju 



\ 9 



Z-' 



V ' ^/2 



f> 



OL?" 



,^4 



log 2 — ygj (ft)), W3) z - 731 (wi wg)^ 



S3 



«■'' 






+ .s, + Is, + ... + .„,,„ 



2 



^OqI 






« < e 



Oft) 



:Ct) 



«i 






ft)o 



) + ^^ 2^0' ( ^ 



ft)i a>o 



1 V 






(2) 



+ 9 3^0 ( ^ 



^1 ^2 

1^ ' (Z 



+ 



where we put symbolically 



00 00 



1 



w.,=o 7/^,=o(^%^i 4- '^%&>3r 



ft) 



//i 



And now, equating coefficients of the various powers of z, 






rwi 



sftJcA 






/ ^1 



«, 



721^.^ 3 



V 9. 



731 (^1. %) 






juTTLIJjp^g QK^Q I 






P' Q 



(Oi co^ 



J /? 



the transformation equations for the first and second gamma modular forms, 
ote that we also have, where m > 2, 



Ti^i.^H'^^-/'^) = {^r(m^i)i[\ 



1 



r=0 s-0 



P 9. 



V' 



CO 



m 



If T = 0)^1 o>i we may put w^ y^i {m^, m^ = g^i (r) ; and now, putting p 
we have 



9^31 [;) - 921 (r) = - ^i ^^^^i^f)'^ 2^t/x,,, "^ 



1, q z^^ n 



And putting p = n, q ^= I 



n 



3 



:93i i^^^) - g%i {^) 



9 x:^ f (2) / 1 



?l 



!.^ 



We get analogous formulae by writing w^ y^^ (co^, a)2) = 9'23('^)* 

The analogy between these results and those obtained in the analysis of elliptic 
functions is obvious. We cannot obtain, however, results which will connect such 



expressions as 



g^i{T+ 1) and g^^ir). 
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Integral FormulcB. 

§ 70. In the general theory of multiple gamma functions the fundamental integral 
formula expresses the fact that the integral of the n-ple gamma function can be 
expressed in terms of (7^4- l)-ple gamma functions of specialised parameters. As we 
have not yet defined the treble gamma function, we cannot prove this theorem for the 
case when ii = 2. In the case when n = 1^ this proposition reduces to Alexeiewsky's 
theorem ('^Theory of the G Function/' § 13). We proceed first to translate this 
theorem into the notation of the present paper, and then to give an alternative proof 
capable of extension to the n-ple gamma functions. 

The G function is substantially the double gamma function with equal parameters, 
the two being connected by the relation 

r^"^ {z\0),(0) = G ( ] (27r)*~^2to co' 2co2 

[" Theory of the G Function," § 29.] 
By differentiating Alexeiewsky's theorem we obtain 

d d 

= 1 log 'Zn — z — a -^ ^ -i- (z -\- a — 1) -r log T (z -}- a) — ~ log G (2; + a), 



dz ^ V ' / d. 



and therefore, writing z for z -^ a, and substituting from the relation just quoted, we 
find 

= _ - + 1 + r^ _ lU^iogr K + cu^^iogia^h) + ^ - loga>. 



^ \ ... f z 



But log r -- ) = log Ti{z\o)) ■— — — 1 ) log 0), 



- ft) / \ ft) 



We thus have ^ i/// (2; | co) + co i//^' (^ + co | co) = S/ (2; | ctj). ...... (i.). 

On integration we have 






log Fj {a I co) da = a log V^ (a | co) — S^ (a | co) -f ^ log --^-^-"--pL ^ 



We may put z -{• a in place of z^ so that 

d d 

{z + <^) v log Vi{z -^ a\ci))'\- 0) Y log Fg {z -^ a '\- (t)\ci)) = S/ (;<^ + <^ | ^)- 

And now, on integration with respect to z between the limits and z^ we obtain the 
extended formula 

log Fj (2; + a j co) dz == (;^ + ^) log F| (2: + "^ | w) + ^ log ~^jt7 r V" 

J 1 2 (^^ -}- ft) I ft)} 

— Sj (;^ + ^ h^) + Si (a I co) — a log F^ (a | co). 

o A. ^ 
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§ 71. There are two alternative methods of obtaining this formula, or rather of 
obtaining the f imdamental relation ( 1 ). 

Firstly, we may directly transform the series which expresses i/// (z | m, co). 
We have, from the definition. 



^h (^ 



a),(o) = z y.^1 (ct), (o) + y^^ (&i, m) + 
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00 
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% 4- ('??^'i + ^?^o) m (r/ij + t?Lcj) ft) 






(i)%^ -f m^)^ ©^ 



Put now m^ + '?n^ = c, a change which is equivalent to grouping together terms 
corresponding to points on the cross lines of the figure. 




There we have 
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1 m 
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Now we have seen in S 28 that 



L. - --^ 

I q 

ft)" 
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rax («) 



We therefore have 



i/zg' (2 1 w) 
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72-3 («) = 7; [7 
ft) 



o 



1 J t/^l' (2 






gwa wg aa i i swssiMai 



6' 



loa: col 



ft) 



1 



ft) 



which is equivalent to the former relation. 

Secondly, we may make use of the contour-integral expression in the following 
manner. 

We have (" Theory of the Gamma Function," § 30), when a is positive with respect 

to CO, 
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Therefore, under the same limitation, 



a 



^— ti;/"^ (a\ 0)) 

00) ^ ' ^ 



9, 



TT 



0) 



— _ :i ^ i. — 1 z 1...- w^ 

(1 - C— )3 



^ 



^3^^^ (a + Oi I Oi), 



JNow, smce 



1 



m = (<^<^ 4- mm) 



it is evident that i/^/^^ (a | o)) is homogeneous of degree —• 2 in a and w. 
Therefore, by Euler's theorem, 



a 



so that 






ai|;j^^^ (a | co) + ^^j^^^^ (c^ + ^^ | ^) ^ 21^3^^^^ (a | m). 



On integrating twice with respect to a^ we obtain 

a^i {a I ct)) + (wi/zg^ (a"-f"<w I co) ::^ ;)(■{ (a | o)), 

where Xi (^ | ^) is a Hnear function of a. 

Changing a into a + <^, and subtracting, we see that Xi {^^ 1 ^) satisfies the difference 
equation characteristic of S/(a|w), so that it can only differ from this function by a 
constant, which will vanish, as we see by making ct = 0. 

We therefore obtain again the relation required. 

§ 72. We proceed now to the analogues, of Eaabe's formula. This formula may be 
written {'' Theory of the Gamma Function," § 8) 



log r^ (z + ct I (jt)) dz -^^ a log a — a + - log - . 



We will evaluate 



log V.2 {z '\' a\ o}p 0)2) rf^ and log F^ (;^ + ^* | <^i ^w^) <^^ 



The method which will be employed is the same as that by which Raabe's formula 
itself was originally obtained ; it wtis, in fact, first invented for the proof of the 
present theorem. 

Let 



Wi 



/(a) =: log l\(z + a 0)^, CO.)) dz, 
J 



rn 



Then 



a 



fia 



Ug {z + a) 

log r^ (a 1 6)) + log p^ ((i)o) + 2m7rt S/ (a j o)^) , 



loi 
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and, therefore, on integration^ 



log r2 {^ + a\wi, (u^dz 



a 

log fj {a I cog) da + a log p^ (co^) 





rw, 



+ 2in7n Sx {a \ m^) + log F^ (2; | co^, 0)2) c?.:^. 

^0 



73. We proceed to evaluate the constant term 



Coii 



log To^ (z I a)j^, C^g) C^: 



z_ 



by an application of the multiplication theory for the case m 

We have 



■ft>i 







log F^ {2z) dz 



1 

2 



-toj 



log F^ (z) dz, 



and therefore, by § 65, 



^^ dz jlog F^ (^) + log Fg ( 2; + ^ j + log Fg (^ + 1) + log T^l^z + 



ftjj + a)3 



9, 



3 log pg (wj, wg) — gS/ (22^) log 2 + 3 (m + m^) 27rt ^S/ (0) 



(i2:{ 2 log Tc)^{z) — log ri(2; | o)^) + log p^ (co^) -f 'Zmirt S/{z | (i)^)} . 




Put a = |a>j, |-%5 and I (ajj + ^2) successively in the formula which we have just 

log Fg (2; + a) dZ; and substitute in the formula just written. We 



obtained for 
obtain 



log T2{z)dz 




Jo 



fOil + tOg 







1 

2 



r^i 







[log T^{z\a)^)dz} 



+ 2m7rt 



i^i (c^i I C02) -^ B^[^ 



Mo 



Si(f 



ft)i 






COo 



+ «3 l0gpi(w2) 



s, 









O) 



2 



+ 3c<)^ log ps(wi, cun) + i-3Si(2cuj^) log 2 — 3o)^{m + ^b^) 27^38/(0). 



By means of the formula 



'a 







log Ti {a j co) c^a =::= a log F^ (0^ ] w) — S^ (a 



) + ft) log ^ - - , 
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we see that 



(It's 

^1 



r 2 i 2 r 
Jo J J 



fa)2 W] + WrJ 

2" r^^"2 







I 

'2 



«i 



0^ 



log Ti (0 1 %) dz 



ft) 



ft)i 



2 log Ti ( ^ I w, ) + '^^ log r^ 



ft) 



O 1 --3 



/ 



ft)i "f 6)3 



^ ' % ) + "'-4 "'^ log r/ "-^ ^ "-^ 



ft)j 4" 6)2 



a>o 



/ 






9" log Ti (a>i I (^g) + -|Si (co^_ I a).^) -- S 



ft)i 



<w? 



Si h 



Cl)a 



'^ 2 



W.) 



S- 



I. t^ 



ft) 



3 



+ <^o loff 



ft)i 



'3 -^'-'B 



1 



'X 



m, 



9 



(f) -j- (f)^ 



T^((o, + -^ 1 0)3) r^ (6)3 + ^ h3)r3(^3 + _-j-_-j I ft),) 



2 



AK^i + ^^l^-dpfM 



and this expression in tnrnj by utilising the multiplication formula? when m = 2 for 
the simple and double gamma functions and the simple Bernoullian function, is 
equal to 



2% + J ) log pi (o)^) + 



(Oi 



4 






ofr + 7>' 2S/(ftii + ft)2 I ft>2, ft>^) 



2coa ' 2 



loes: 2 



+ |S/{o|w2) - S, r;riw 



" 1 2' I "'■^ 



+ ^logr,(H^^ 






y 



r 

1 



r„3(-^|wJ/5^3(ft,3) 



1 



+ -rf log ^ - , 



1> 



L 



______ 



«3JA'^- 



o/ft)g 

2 



, |-+Vl0gp2(w3,W3). 



ft)c 



J 



Now, on making oi^ = oi^, we obtain from the multiplication formula 



loff 



(Ur, 



^3^ "9 I ^3) Pi(«a) 



3 log ^.3(03, 0)3) + [38/(0! wj)- 1] log 2. 



Therefore the expression which we have just found reduces to 



2«3 + t) 108/^1(^3) + 



ft)j ftjg ft)i 



9 



4 



'*^1 l_ '^'3 Q ' / 1 \ 



+ I S/ (0 I ^^ 



Q / ^3 



1 



2 



+ 2 -A'C^l^* %) 

%) + 3%l0gP2(ft>2, ft)a). 



► loff 2 



And therefore 



3 



Wl 



JO 



log V^{z\ ft>2, ft)cj)c?2; 

^ 3ft>i log p2(ft)j, ft>^) + Sft^^ log p^^, foa) — 3ca3 log p^io^^ 



3ft>j(m + m!)2wh 281(0) + (1 + 2m?rt) 



Si{?l'-2) + |B/(o|a,3) 



Q 
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since the coefBcient of log 2 is equal to 

*■.■» r.% #'.•» (On / 



«i: I ^1 _ :?!3 I ^3 

iiiCOc) ~c Zi ^J 



We have, then^ finally^ 



'to. 







log r2(2;)(i2; = W|[log p2(0J|, OJ3) 



(m + ?n^) 27^283 '(0)] + ojglog ^^^--"-- ^' 



We thus see that substantially the double Stirling function of m^ and o)^- is 



expressed by - 



a?i J 



log Fg (z) dz^ or, by symmetry, by - log To (0) dz. We have, in 







fact, the relation 



Psl^b %) -- 0^^ + m')27rt 38/(0) 



1 p' 

^1 Jo 



ft). 



log r^ (z) dz ~ -- log A 



W^ 



+ Z' [log -2 



2m7n 



where A is the Glaisher-Kinkelin constant (" Theory of the G Function/' § 3) ; and 
therefore 



/>2(^n %) - (^^^ + ^/-^^T^^^S/ (0) 



Ct)i 



a>| — 0)3'" Jo 



■^1 
log Tg {2;) (is: 



COg 



O).) 



ct)]_'" — a>3^ J 



g Tg (z) & + |;j-^^^ 



^1^2 



1 



OP^ 



Wo 



ft>2 ) "^ a>i 



for^ by § 22, log o).^ — log oj^ —2{m — m^) tti = log --^ 

§ 74. We may readily prove these results by the relation which exists between A 
and p, (co). [The latter is a convenient way of writing p,{a>, c.).] 

For, when each of the parameters is equal to co, we have (^^ Theory of the 
G Function/' § 29) 



V2 -^ (2: I Cu) = G ( - ■ ) (2it)' 2co M 2w» - 3 



Z (z — iof' 



\0) J 



and therefore 



G ii) = Ti 



ft) 



^^ 



(i) 



{27r) 



1 «„^ 5 



Now from the multiplication theorem for the double gamma functions^ when m 
we have 



2 / ft) 

J-2 \ 2 



xj.exiot> 



ftj ) f — 1 =^ p^ (co) 2/ 



12 



G(i) = ^|y^V3-'(a,) 



-* 
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Now it has been seen that {'^ Theory of the G Function," § 15") 



A=: 



2 



""^e^. 



and therefore 



7r"G=*(-|-) 



Pi (^) 



which is the relation which was used at the end of the preceding paragraph. 

§ 75. It is interesting as a verification of the algebra to notice that Alexeiewsky's 
analogue of Raabe's formula {'^ Theory of the G Function/' § 1^) yields the result 
of § 73 in the case where the parameters are equal to one another. 

This theorem is expressed by the formula 



logG(. + l)dz = |logG(i) + J^log^ + J«log2 



1 2 5 



and therefore, utilising the relation between the G and double gamma functions, 



'03 

log T^ (a + oj 1 0)) da 



"" log G (1) + ^ log 27r + -^ log 2 + ^ - ^"^ log 






3 



CO 



If now we express log G (|) in terms of log p^ (w) by the formula of § 74, 
we find 

log Tg (a + w I 0)) da = ^ log p^ (w) + w log ^~j- + :^ , 







12 



a formula which is equivalent to the result of § 73. 

§ 7&. By combining the results of §§ 72 and 73, we may now write down the 
value of 

log V^{z -{- a\ o)^, ojg) dz . 

• 

For we have seen in § 72 that this integral is equal to 

ra rwj 

— log TjL (a I wg) + a log p^ [co^) + 2m7n S^ (a | cog) + log Fg [z \ w^, w^) dz, 

which expression in turn is equal to 

— a log r;^ (a I 0)^) + Sj, (a I w^) [1 + 2m7rt] — w^ log T^ (a + w^ | w^) + (a + w^) log p^ (co^) 

+ ^1 [log p, (^i, ^2) - (m + mO 2.r6 ,S/(o)] + o;, log "'^^^" + 19(1 + ^^^O- 



12 



Thus 





log Fg (2; + a 



^1? <^2 



a log 



Pi (^2) 



+ ^1 [log/>2 (^1. ^2) - {^ + ^0 27rt^S/(o)] 



+ 1 + 2m7rt) S, (a I co,) + j^l ^ co, log -~~7~~T^ 
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Since, by § 64, 

log p^ (co^, o)^) - (m + in) 2iT6 ^8/(0) = i log yiy^y., + 1 [l ^ Sii^^Tl % 2 
we see that this formula may equally be written 



^ ■ 



log r^ (2: + a 



C0|, OJ2) Ct: 



^ 



a 






§7172X3+ Tr[l -- aS/(o)]log2 



-f" (1 + 2m7n) i Sj (a j cog) 4- 



03.) 

12 



o los: -"^ — --^1-^^ . 

P2 (^2) 






This and the corresponding formula, obtained by the interchange of co^ and coo, 
m and t?/, are the analogues for double gamma functions of Raabe's formula for 
simple gamma functions. 

I 77, In the particular case when a is positive with respect to the m8, it is possible 
to obtain more simply the value of 



'Wl 



Jo 



log r^ (a o>|3 (i)^) d 



a 



by means of the contour integrals investigated in Part III. 

We give this method of proof as it leads incidentally to an expression as a contour 



integral for log 



Pi {^?y 



We find, on integrating the expression for log ^ "^AiL__.: gi 



pc,v«i, ft)3 



given in § 45 



Wi 



- 



log Tg {a) da — oy^ log p^ (oj|, cog) 



(M + m + m') 2,rt .So {«) c^« + «i ,S/(o) 2M7rt --- -f f (=-f );^ ^^^^^X " ^^H ^^^^^ dz. 

AiT J L (1 — (3 ^■^^) 







and the right-hand side, by an application of the formulae of §§ 6 and 44, becomes 



(M + ^ + '^i"^') 2?rt 



<^iSi{o) + 






27r 



■+ a)| 28/(0) 2M'?ri 
(•-^ ^^)~" {log (- z) -^ (M + mO 2t. t + 7} 



^1 



(1 — ^-^2^^) 



C/2; 



NoWj reducing the contour to a small circle round the origin, we see that 



27T 



(M 4- m') 27n 



(—- Z)"~''^d: 



S./{o\m,) 



. I 1 



e 



— (anZ 



(M + m') 2m 



M2 



and therefore 
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log Fg (a) da — o)^ log p^ (oj^, w^) + <^i (^'^ + '^^^0 27rt 281^(0) — ^2776 



S/(o|ft)2) 



27r J 1 

<03 



( - g) ^{lo g(~^) + 7 } 
1 — e~^3^ 



dz . 






Cf)f 






we see that, to establish the formula of § 73, it is necessary to 



show that 



L f (- ^)""Ml0g(~ Z) + 7} 



'TT J 1 



1 



e 



MZ 



d: 



Z 



0) 



^ p;j^ (O)) 12 



OJ 



This may be readily done as follows : 
We have 



T Vh(a I m) 
lo2""-^-'^^^ ^ 

Po(ft)) 



Stt . 1 



^-«^(— zy^ { log (— .^) +■ 7} . 



(1 — e'~'"'^f 



dz 



10 



Integrate with respect to a between co and 2a), and we find 



r2w 

log Tc^{a \(t})da — co log p^{co) 



e;-"'^( — ^)-3 { log ( — ^) 4. ry} 



w 



27r J 1 1 — (3""'^^ 

to 

, r (_-,^)~2{iog(-^)-i-^} 



'ia; 



9. 



^ir J 



1 



1 — e 



-OiZ 



dz , 



for it may be readily seen that 



'TT . i 



log(— iS) + 7 



cfe .= 0. 



z-" 



(x> 



[This vanishing contonr integral occurs when Eaabe's formula is proved by means 
of the expression of the simple gamma function as a contour integral 
But as we have deduced in § 75, from Alexeibwsky's theorem, 

logT^{a a))cta— wlogp2,(a)) ='colog~~y~ + -. 

The contour integral has, therefore, the required value. 

We here conclude our investigations of the algebra of the double gamma functions. 
It is evident that the formulae admit of still further development ; they lead, for 
instance, to many curious relations between the integrals obtained in Part III. Such 
considerations are, however, foreign to our immediate purpose. The development of 
the integral formulae and the theories of multiplication and transformation in the case 
of the double Eiemann C function is interesting in that we thus combine many of the 
formula which have been obtained separately for Bernoullian and gamma functions, 
and the algebra by which such developments are obtained by the extension of 
M Ellin's definition of the simple I function is in many ways attractive. Owing, 
however, to the length of this paper, we do not propose to consider it in this place. 
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Part V. 

The Asymptotic Expansion and Transcendentally-transcendental Nature of the 

Double Gamma Function. 



78. There remains now the consideration of two more general characteristics of 



the double gamma function : — 

1. It admits over part of the region near infinity an asymptotic expansion in powers 

of the variable. 

2, It cannot be obtained as the solution of a differential equation whose coefficients 

involve exclusively more simple functions. 

It will afterwards be seen that these characteristics are common to all gamma 
functions. 

Let us consider first the behaviour of r2 {z) near infinity^ We know that its poles 
are given by 



Z + '^'h'^i + ^'^'^2^2 ™ ^ 



^n^ 



m^ 



0. L « . . CO 



I 



3 -'- 3 




Therefore near z :==: co the poles of V^ (z) are massed together between and on the 
negative axes of w^ and co^, so as to form a lacunary space on the equivalent portion 
of the Neumann sphere. Between these axes, therefore, an asymptotic expansion 
cannot represent the function. We have to consider whether such an expansion can 
exist outside this lacunary area, within, that is to say, the non-shaded portion of the 
figure. 

We shall in the first place proceed entirely algebraically. It will be proved that 
within this non-shaded area a quasi-Laurent asymptotic expansion of the form 
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exists, and then it will be shown that, the possibility of such an expansion being 
established, its actual form is readily obtained by a process of difference-integration. 
Subsequently we shall verify the results by an alternative proof by means of contour 
integration, this proof being the natural extension to double gamma functions of the 
one employed for functions of a single parameter ('' Theory of the Gamma Function." 
Part IV.). 

It should be noticed that the expansion under consideration differs materially from 
the expansions obtained in Part III. In those expansions the limits of the number of 
terms of the series and products, the quantities pn and g?i, formed the infinite basis 
terms ; but in the present case that basis is the variable itself 

§ 79. We first write down the asymptotic expansion for 

log r^ (z + a 1 CO, 0)) . 
We have obtained (''Theory of the G Function,'' § 15) the asymptotic expansion 

Z -4- Cl X /(z -4" ft ■"** 1^^ \ 

log G(2; + a) = i\ - log A + — ^^ - log 277 + ( ^-— ^ _ -A- j log 



z 






+1 



n 



=.1 nz^ ^^\ I ' / « Zi'in {2n + 2)z^^ 






this expansion being valid for all values of a and z such that | ^ | is large, and the 
principal value of log z being taken. There is, in the language previously employed, 
a barrier-line along the axis of — 1- 
Now in § 70 it has been seen that 

log g(^^ = ^ log r,(^ I o)) + ^ log 27r -- i^^^ + ij log 0), 
and log A — ^^^ == log "^/^ + \ log co. 



Fi(«) 



Therefore 



T VJz-\- a\(o) (z-\-a'-'(of , \ , ?)Z^ / a 

— log -^ r— ^ — ^r~~^ — -A- 10£C 2: — 7"-:; — 2; -T, 



1^ 
ft) 



+ 1 t:)::^Ls4aio,)+^^ 



- 



where now there is a barrier-line along the axis of •— co. 
so that we have finally 
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loi 



r/ 1 "\ 









00 



f7 , 






,n = i m (m + 1) 



-V 



with a barrier-line along the axis of — co. 

When a :=: w, we see that we have tlie asymptotic expansion 



lo 



VJz-\- ft) I CD ) 



Pzi^) 



1^ 

l"2 






(-)'--^' B,.,2(«) 



mojz^^'' 



§ 80. We will now prove that, provided z be positive with respect to the oi's, there 
exists an asymptotic expansion for log r^ (z) of the form 



■ (1, zY log z -f (I5 zY + 

where (1, s:)'^ denotes a quadratic function of ^. 
In the first place it is evident that we have 



OD 



A = l 



Xr 



/'jV 



z 



CI -y- njG)^ — p lljoOt).-}' 



where a is some finite quantity, and w.j and t2^ are singly or together large positive 
integers. 

Now, from the fundamental diflPerence equations of the double gamma function, it is 
at once seen that 



log r^ (a) — log r2 {a + "^h^i + ^^2) 



n^ no 



log n n (a + ^'^i^i + ^2^2) H" ^^^g ^ 



'!j. Fi (ct 4- ^iO>] I Wg) 



«% = 



Fi (^^2) 



+ log n — -^^ — — _-- -.ii--L _ 2mirrt >v bi (a-f- ^^'^i^^i (^h) 



m2 = 



7l9 



2m Vt 2 Si'' (a + '^^^ 

m2=0 



Cx) 



S 



) 



A term has been neglected which is an integral multiple of 27^, and which is there- 
fore absorbed by a suitable specification of the logarithms involved. The above 
formula may be rewritten 



loo' 



Tj (a + n^«j + flgOg) 



log 1\ (a) — log 



n, 



Tj (a + %ojj I mo) 



+ 



•51, 



log r^ (a) — log IT 



Pi («i) ! 



Pi {^^2) 



n-t Uo 



log r^ (a) 4" l^g n n (a + ^^'^i^i + "^"^h^h) 



riii~0 rii2—0 



+ a (1, r/)^ -|~ (1? ^^)^ 



In the first place, if v/e put ^^l in place of p7i and n^ in place of gn in the procedure 
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of the paragraphs leading up to § 49, we evidently obtain for 



«1 '«2 



log T, (a) + log n n (a + m,co, + m,a>,) , 

11 ... 

an expansion in powers of - and - , every term of which involves a algebraically, and 

of which the non-ultimately -vanishing terms are typified by (1, ny^logn + (1? ^0^ • 
In the second place, consider 



log r, (a) - log n lA ± ;h^>i). 



r/ii—O 



Pi i^d 



We have seen, in § 30, that 



To{a) — 6-^-1 ~a{y,, + ±iogo.,^ ^^]T^{a\a)^) 



00 



•\ xX 



m^ = 1 






the product being absolutely convergent. 
The typical term may be written 



E 



xp. 



,4. 



* i/; W (m-^o)! ! 0)3) + j7 iij'^-%mj<Oi \oj^) + . . . 



O i 



4! 



Therefore {'' Genesis of the Double Gamma Function," §§ 4 and 5) 



Wi 



log Tg (a) —■ log n Tj (a 



■}/ij = 



mjWi 



0) 



.) 



admits an asymptotic expansion of the form. 



00 "^ 

(1, n^^) log n^ + (1, n^f + 2 ^^r. 



r-l % 



each term of which involves a algebraically. 
Combining these results we see that 

log r^ (a + 7/iCOj + 71^0)^) 

11 
admits an asymptotic expansion in powers of - and — , each term of which involves 

a algebraically, and of which the terms which do not ultimately vanish are typified 
by (1, n)^log n + (1, '^0^- 

But log Vo^(a 4- ^^i<^i + %<^s) is ^ function of a + ri^o)^ + ^^2^2- ^^ must then be 

1 
capable of an asymptotic expansion in powers of — , each term of which 

involves a algebraically, and of which the terms which do not ultimately vanish are 
typified by 



376 



MR, E. W. BAENES ON THE THEOEY OF THE 



And now by mere re-arrangement we may include a with the term niO)^ -}- n^ct)^ each 
time that the latter occurs, and we obtain for log Fg (z), where z = a -j- n^w^ + "^-2^2^ 
an asymptotic expansion of the form 

(1, 2:)'^ log Z + (1, zY + t -;7 . 

§ 81. We can readily extend the previous proof to the case where z lies between 
the axes of — co^ and co^, so that it is given by 



z 



a — ^^1^1 + ^^2^3 ° 



By writing the fundamental difference equation in the form 



T,{z ^ CO,) = T,{z) 






^— 2i/i7r(,S_^' (2 — Ml \ co^) 



we readily see that 



log Tg {z - n^o)^ + n^wg) 



logr,«-S"log'>±5f^>' 

Pi \^l) 






+ 



log r,( 



wa - 



2: — W- 



71, 7r9 - 1 



<^1? ^5) +2 S log (0 — ?72^(i)| + ^^^2) 

7«.i = 1 J% ^0 



Wi 



log To (2; — - (ij J -— C0|, oi^) + S lo 

■nil =:= 1 



T^(z — myco-,\(0o) 



Pl((^2) 



n. 



n^ - 1 



mi = 1 7)12 — 

But by the theorems just quoted in the previous paragraph the three expressions 

1 
in the square brackets severally admit of asymptotic expansions in powers of — and 

— , whose terms which do not ultimately vanish are typified by 

(1, n)^ log 71 + (1, nY 

and whose coefficients all involve a algebraically. 

Thus, by a repetition of the previous argument, log Tc^(a — /ij^oj, + ^^2^2) admits 
when z = a — "i^io^i + ^2^2? ^^ asymptotic expansion of the form 

{l,zflogz + {l,zY + t ^-l' . 

A " 1 Z" 



In an exactly similar manner we may show that log 1^2 (^) ^^^^^ admit of an asj^rnp- 



totic expansion of the same form when 
and 0)^, 



z 



is large, z lying between the axes of — w 



2 



DOUBLE GAMMA FUNCTION. 



377 



S 82. But when we come to the case of the negative quasi-qnadrant given by 



% 



a 



n^mi 



%^2' 



it is interesting to notice that the above proof breaks down. 

As before, by the use of the fundamental difference equation, we obtain the 
relation 



log I\(a — n^o}^ — n^^ -» log l\^{a) — log l\{ 



a 



0) 



0) 



1? %/ 



log r2(^* - ^2 1 ^1. -" ^2.) - log V^{ 



a 



(I) 



CO. 



OJp 



OJo) 



log l\{ 



a 



0) 



a)|, 0) 



^) + i log IV« - «v»^%)- 



'I'll, = 1 



f 



log r2(a — ojg 



no 



O) 



2? 



<j) 



l) + 2 log 



m.-, — I 






loQC Vcy{a — CD, — • w 



2 



^2 «2 



(jt)|, — ^ O) 



2) + S S log (a -- TOiMi - iMgai,) 

-//l^ = 1 rill =: 1 



■5?! 



n.-, 






O) 



,) 



TlU — 1 



The several expressions bracketed on the right-hand side of this identity admit of 

of which the terms involve a 



1 1 

asymptotic expansions in powers of ~— and — 



n^ 



71, 



algebraically ; and therefore the whole of the right-hand side admits of an expansion 



of this form. But there remain the non-algebraic terms 



log V^{ci) + log Vo{a -- oj^ 



W|, 0J2) + log f'2 {^^' "^ ^2 1 ^1? 

+ log I\ {a - 



0). 



:) 



0) 



O) 



a) 



1? 



0) 



2)? 



and when we seek to group — n^o)^ — ??2^2 with a, we are forced back on the original 
function r2(a — - n^ct)^ — n^o)ci). Thus as regards the possibility of an expansion, when 
z is negative with regard to the w's, our- results are, as we should expect from § 78, 
entirely negative. The region between the axes of — - coj and — o)^ is a harrie^'-region 
for the asymptotic expansion of the double gamma function. When W| = m^ this 
region closes up into the barrier-line which occurs for the G and simple gamma 
functions. 

§ 83. We can now find the asymptotic expansion of 

I r^ O' " '^ 

1 o \'<^ 6t)|5 ftJoj 



for large values of \z\ which are such that % does not lie in tlie barrier-region, a being 
any complex quantity of finite modulus. 

For such values of z and a we have the expansion 
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^S^'y^ = t^(^)^ +fM1H ^ + <^, («)^ + M<^) +,iyTTj 



where f (o) and <^(a) are algebraical polynomials of degree indicated 
suffixes, and Xr(cO is? ^^ long as r is finite, likewise an algebraical function. 
Now by the fundamental difference equation 



by the 



leu' 



loo" 



Vc,(z 4- a + cdi) 



log 









Again ('^ Theory of the Gamma Function," § 41), we have the asymptotic 
expansion 



Pi K^2) 






where log^. z is that natural logarithm of z which has its principal value with respect 



to the axis of — oj^. It is thus equal to 



log z — 2mVt, 

the latter logarithm having its principal value with respect to the axis of ~-{(i)^"\~ w^). 
We have then, if log z have its principal value with respect to the axis of 



[f,(a + o)^)z +/;(a + CO,) -fi{a)z ^f^{a)']logz 

CO 



, %.(t^ + 6 )j ) - y,(ft) 

r (r + l)z'' 
r= "^ 8/(2: + a\a)^)[logz - 2{m+m')7rL] + z8,^'^{z + a\oj^) + S --";:r^-7-^y^^ - 

If we equate corresponding powers of z on both sides of this result, we find 

Xr{a + 0)1) — Xri^l) — SV + i(6^ I Wg), 



and similar relations among the /'s and cjy's. 

We shall get, in like manner, another set of relations in which w^ and w^ are 
interchanged. Remembering that the/'s, i//'s, and x^ are all algebraical polynomials 
which vanish with a, we thus prove that 



fi{a)z+f^{a) : 

Xr{a) = 

By this process, which may appropriately be called a process of finite integration, 
we obtain the asymptotic expansion 



S^{z) ^ gSo(2: + a) 
^{?So'(^ + cO -- gSoX-)} 
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+ zlSoV^ + a) - SoW + S ^=^^ , 



log ^ having its principal value with respect to the axis of — (ct>i + ah). 

§ 84. We may now obtain the asymptotic expansion for log r^(2; 4- <^)j nnder the 
limitations assigned at the commencement of the preceding paragraph. 

For this purpose integrate the relation just obtained with respect to a between 
the limits o and Wp 

Then, by the formulae of §§ 12 and 76^ we find 



Pl(^2) 



log/>2(o>i, %) — (m + ni)27n.2^{{o)[ 



(1 + 2m7n)\ Bi{z\a)^)-{' 



ft); 



12 



m^log 






^1 % TsC^^) 



1 .,r a /v I \ C! /.A a ^/.A1 i ^ (—") L2^mn(ft >i) — C«)| gb m +](^)J 

ubstitute now the asymptotic expansions lor loff -^-j-^-^- and log -^ z — . . 

of which the former has been quoted in the preceding paragraph, and the latter 
obtained in § 79. 
Then we find 



cO| log 






(1 -f 2m7n) 



Bi{z ft)o) + 



12 



S^(2:|%) + 



S/(() j 0J3) 



5> 



Ct) 



aSiX^) 









w loff..^ 2; 



{log,,+<o,^ — 2(m + m')7rt} 










ft>c 












00 



'^^ 



M 



1 



%=1 



•,i=i (w + 1) « 



B« + 2(«3) ~ «3 log»„ 






1„ 

i 2 



J 



o^^ 

o;^" 



00 



M 



lB,^+3(ft>3) — ft>i jjS to + l(c>) 



/>.?* 

fi/ 



Remember that log;,^^ = log,,^+co,''2 — 2mTn; then we obtain by an easy reduction, the 
asymptotic expansion 



log ^ 



P2 («1> "s) 
28/(2) {l0g„^ + .,2 



2(m + 77^') TTt} + 2 2S/^> (0) + f , 3Si(^> (0) (1 + 1) 



OD 



+ 2 



(-)»3S',„ + i(0) 



m = l WiCwi +1) 2" 
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iV; 



is iarp^e, and z does 



tt> 



v/hicli is the complete asymptotic expansion for log F^ [z) when j z 

not lie within the barrier region negative with respect to the axes of — o)^ and 

If we combine this result with that obtained in sS 8?) we find the more general 

'■J o 

formula 



log 



'-Si (2^ + «) (log.,.-.. 



z 



2(«i+m') ,n}+ z J^f> (a) + J .S,<«' (a) {{ + l) 



2' 



7« ~ 



1 1 m (m + 1) z"^ ' 



valid under the assigned limitations. 

The expansion is written in the precise form ado])ted, in order that the analogy 
with the corresponding formula in the theory of multiple gamma functions may be 
more clearly displaj^ed. 

§ 85. We might now conclude this investigation^ Since, however, this would appear 
to be the first time in analysis in which an asymptotic expansion with a barrier 
region has been obtained, it seems better to give an alternative proof which shall not 
need the difficult argument of §§ 80-82. This proof is the direct extension of that 
previously given for the case of the simple gamma function. "^^ We therefore proceed 
as briefly as possible. 

In the investigation of § 57 it was shown that when \s\ is finite and |i (^) > — ^5 
where k is a positive integer, the series for 4 {^, a\o}^^, wg) ^^^ absolutely convergent. 

Suppose now that s :== cr + ir, where cr > — ^:, and suppose further that z does not 
lie within the region bounded by axes to — w^ and — oj^, and that a is positive wuth 
respect to the mb. 



hen, smce 



1 



CO- 



■ CO:-, 



1> k -^ 1 

s^^^., jc V^^) — ^^ _ ^^^ ^2 ■-- s) w^ 0)3 a' -' ''■ 2(1 -- s) w^ a>2 a' ■ 1 ^ c(^ ^ 



.. \ 



'^ ^ 3 1 1 1 



it is evident that, if ^9 be any positive integer, the absolute value of each term of the 

expression 



sPz' 



sin ITS 



2^-5 



,, {a 



0J|, a)g 



tends to zero as r j tends to infinity. For, by the restrictions on z and a. 



and therefore 



(f 



(r e/'') 



+ IT 



whei'e o < ijj < :i:. 



IT. 



■^' " Theory of the Gamma Function," Part IV, I regret to say that the Lemma of § 40 is faulty ; the 
theorem is evidently only true when a/w is real. A slight modification will, however, establish the truth 
of the main proposition under the conditions enunciated. 
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:V 



a^+^ sin TTS 



^-v^r-Trlrl ^ ^.p J^J^ 



where M is finite, however large 
Hence the absolute value of 



m 



ay be. 



.V ^s 



sin ITS 



I S-~s,k{^^\^^l'> ^^2) 



tends to zero as | r ] tends to infinity ; and this theorem is true if a is rejilaced by 
a + 'in^oy^ + ^'^i<^-2? where mj and in<^ are positive integers. 
Hence the absolute value of 



s^r 



sni ITS 



'Xi^^' "^" '^^^^^1 "~t" '^^'^2%1^'^> ^^) 



(where x is the function introduced for brevity in § 57) tends to zero as 
to infinity. 
But 



4(,s', a I ct)i, coc^ = gS^,,/,(a) 



and therefore 



CO 



00 






s, k), 



t\ tends 



00 CO e/';^"''' 



Now the double series on the right-hand side is absolutely convergent for all 



finite values of 
to infinity. 

Therefore 



, and the absolute value of each term tends to zero as \t\ tends 



S'P7ft^{s,a\(jd^,(D^ 



sin ITS 



remains finite as | s j tends to infinity, ^i(.9) being finite and not greater than 2. 
When %(s) is greater than 2, we have 



^z^^^^ (5, a I Wp Wg) 



sm ITS 



00 CO 



c 



'sP 



^ '^ 



and therefore the expansion on the left-hand side is finite however large \s\ may be 
provided \z\ < 1. 

§ 86. Consider now the integral 

27r6 J 



S sill TTS 



The subject of integration is a uniform function of 6', wherein z^ is to have its 
principal value with respect to the axis of — (co^ + co^)^ ^ is to lie within the region 
bounded by axes to — oj^ and cu2, and a is to be positive with respect to the oj's. 
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3 



"(HH) 




111 the first place let the contour be taken to lie along the real axis, passing from 
+ CO to + ^, a,nd cutting the axis between the points o' =^ 2 and cr ^==^ 3, as the 
contour 1 of the figure. This is equivalent to taking the integral round a contour 
enclosing the points 3 and + o^- 

When 1^1 < 1 the integral is, by the theorem of the preceding paragraph, finite ; 
and by Cauchy's theorem it will be equal in value to the sum of the residues inside 
the contour. 

Now by § 53, when 5^2 + /c, where k is an hiteger, 



4 (^. ^0 



i-y 



(/H"i) 



^jj^^'^ia) 



Hence the value of the integral along the contour 1 is 



CO r/2 + k 

rV 



s 



o(/i;+2)! da~ + ^- 



fP + i 
¥Ii lc>g r, (a) 



and by Taylor's theorem this expression is, under the assigned Hmitations, equal to 

\ (z '+ a) 



1 



Oi 



rs (a) 



^,<» (a) 



o 

^ 



Let us now make the contour expand until it becomes a straight line perpendicular 
to the axis of cr, cutting the axis between the points 2 and 3, and a half circle at 
infinity. The value of the integral will be unaltered, since the contour in expanding 
passes over no poles of the subject of integration. And by the theorem of the previous 
paragraph the part of the integral which is taken along the semicircle at infinity 
vanishes. Hence the integral along the perpendicular line (the contour numbered 
in the figure) is equal to 

log- — -^- — r-" -^ z^jjJ^Ua), when Iz < 1, 

1 3(a) / ^ \ / ! 

But the integral and this expression both remain continuously finite when \ z | becomes 
greater than unity. They are therefore equal to one another for all values of | z 
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Let now the perpendicular contour be distorted into a contour which encloses the 
points 2, 1, 0, . . . — n, and which after the point o- = — ?i again goes off to 
infinity perpendicularly to the real axis. 

So far as the value of the integral is concerned this contour will differ from the 
second contour only by two strips at infinity of length less than (n + 3) parallel to 
the real axis : and by the previous paragraph the integral along these strips will 
vanish. 

By Cauchy's theorem the value of the integral along this third contour will be 
equal to minus the sum of its residues at the points 2, 1, 0, . , . —n, together with 
the integral along the perpendicular line cutting the axis of or between the points 
— n and — (n + 1). 

Now, when 5 — 1 or 2, the residue of the integral is equal to the coefficient 
of 1/e in 



e \ s /It (s 



(_-^).5+l 



(-y+^ A 



1) I'' '^ \ / ' ^^^ __ 2) 



"Y" ... "T" 



1 \ 



s 



^jsr'\{<^) 






1 'yS 



and is therefore equal to 



^iyi i^'^\<^) + (-)■'"' 2(m + m') 776 ,S/'+'){a) } 'f {1 + elog z -\- . . . } 

by I 53, 






s\ 



'X 1 

11 • • • I 

6- 



log 2; 



% p- 



.S^(-^^){a) + ^^[t/.,^^) (a) - 2(m + mO ttl ,S,^-^) {a)\ 



where the logarithm has its principal value with respect to the axes of 
When e =: 0, the residue is the coefficient of 1/e in 



(^1 + %). 



'^---—-■- ---- - \ gS'j (a) + € log -~:^^l-^er^^^''^'^^-^ 



po{o)i,0)c^) 



and is therefore equal to 



V, (a) 



oB\ (a) hcf z + kw --^y^^'^'^^^ ^-A(«)2(m + .H0..^ 



When 5 



w, the residue is 



We therefore have 






. VJz -f a) 
lop- ^'^~~^-:~-^ 



^ i\ (a) 1 1 ^-^ ^ ^ 



loo- i ^^^^^- g.S',(o)2(w-fm')m 



by I 60. 



i 



2^'lict) [logZ — 2(m + 'i^l)TTL 






s 



V (~)"'3'-'' «H-1 W , 



« « 



,=1 m(m 4- ]);>;^'* ^ 



Oi,, Wo) 



s 
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where the quantity log z has its principal value with respect to the axis of — -(w^ + o)^) 
and J^, (2:, ajo)^, co^) is equal to the fundamental integral taken along a perpendicular 
contour cutting the axis of a between — n and — (n + 1). It is evident that the 



integral when 1 2; | is large is of an order of magnitude less than 
We therefore have the asymptotic expansion, when j z is large 



1 



n 



log 



Pi (Wl, Wo) 



^ 






1 
1 



1 

2. 



z 



Zi 1 



1 
1 



,S/ (a) [log 2 - 2(m + to') ^t] + i ^ ^ '"~^''" + i ^''^ 



m 



I m(m 4- 1)?^' 



and the residue after fi terms of the final series have been taken is of the same order 
of magnitude as the final term taken. 

This expansion is evidently the same as that previously obtained. The limitation 
that a must be positive with respect to the o)'s may evidently be removed by 
employing the fundamental difference relations for the double gamma function and 
the asymptotic expansion for logr|(^ -{- a). We are finally left with the essential 
limitation that z shall not lie within the barrier region bounded by the axis to —o)^ 
and •— coq. 

The IVanscenderitciHy -transcendental Nature oj T^iz), 

§ 87. We finally prove the theorem that the double gamma function cannot arise 
as the solution of a differential equation whose coefficients are not generated from the 
function itself Modifying slightly the nomenclature introduced by Moore/^' we 
may say that Tc^iz) is a transcendentally-transcendental function. The proof is a 
slight modification of that given for the G function (§ 30), which in turn was similar 
to the investigation of Part V. of the "Theory of the Gamma Function." 

In the first place it may be proved exactly as before that if the theorem is true for 



d 

cl 

tion of the function 



-log T^{z), it is true for l\^{z). We shall therefore confine ourselves to the considera- 



4>(z) 



cP 



clz 



:a l^g ^'2 (^) » 



By the fundamental difference equations of § 20, we have 

(j) {Z + 0)^) - rjy (z) — Xjj{z\ a)^) 
<f) {Z + Wg) — ^ {Z) == .1// {Z j U)y) 



cP 
where, for convenience, we put i//(^) -=■ — log \\{z). 

* MoORE, 'Math. Ann.,' vol. 48, pp. 49 ei seq. Moore uses tlie term only to describe functions which 
cannot be generated by a differential equation with algebraic coefficients. 
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Suppose that 

y =z cj) (^x) satisfies the, differential equation 

so transformed that it is rational and integral in y and its derivatives. 
Let the terms of class s be symbolically 

Ro(^) Q^^ (y)^ ^j {^) Q'^iy)^ • — ^^/c (^) Q/ (y) ; 

in terms of class (5 —- 1) being 

So (x) Q\.M ,8i{x) QLi {y)> 

and the functions R(x), S{x) being holomorphic. 

If (f){oc) satisfies the differential equation, (jy{x) + xIj[x\ oj^) will satisfy the equation 
in which (x + a)i) is written for x and 4^{x) + \jj{x\o)i) the equation in which [x + oj^) 
is written for x. 

Make the first substitution, divide the equations by Ro(''^) ^^^^ ^^{x + oj^) respec- 
tively, and subtract one from the other. We find 

|||q.'»W] + ■ ■ + |;||q,'[*w] }+q,°ww +*(=t|<»,)] -q.»[.mx)] 

+ terms of lower class == . 
But Q,o [cl>{x) + x}j{x I CO,)] - Q,o [(/>(^)] 

consists solely of terms of lower class than s. 

Hence either the equation, which has been obtained vanishes identically, or we can 
reduce the equation for y to one in which there are fewer terms of class s. 

The equation cannot vanish identically unless the coefficients of the various terms 
of class s all vanish, which necessitates that the ratios 

E-l(^) Bk (x) 

are doubly periodic functions of x of periods w^ and cog. 

The equation for y can thus be always reduced to one of the form 

n{x) [p^{x) q,%y) + . . . +pi{x)Q/'{y)'] 

+ SoHQ,"_i(:y)+ . . . +Mx)Ql-M 

+ terms of lower class = 0. 

where all the coefficients are holomorphic functions, and, in addition, the functions 
p [x) are doubly periodic of periods coj and o),. 
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Divide the equation by Ji{x) and subtract it from the equation which results from 
changing x into x -{- Wi^. We obtain 



QJ'iH^^m 



+ ^^--'^ QUiH^) +^(^h,)}^ lt^QLi{</>(.^)} 



Ji(x + 



CO, 



E(a;) 



'Y' .... 



i j ii tt|»>"i 



E (^ + 0)^) 



;qU{<^(^)+^(^|co,)} ^^^QU{</>H[ 



+ terms of lower class = 0. 



This equation will not vanish identically unless the functions 

^, 7V , . . . :rr7~r a^U satisfy relations of the form 
E(a^) ' E(^) ^ 



n 



f{x + o)y) —f{x) - t pt{x) Q/ixlj{x I a).^[ 



7c ^0 



and therefore by symmetry relations of the form 



n 



f{x+ wg) —f{x) — % 'Pi (x) Qi^ 



i-=ti 



l>(x|6Jj) 



The quantity {n -\- 1) which gives the number of terms on the right-hand side of 
these two relations will not vanish unless the original equation can be reduced to the 
form 

■ + P'2,0 {«^) Qti (y) + ■ ■ ■ + p&, Hx) QLi (y)] 

+ To (x) QU {y)+ . . . + T,, (x) QS (y) 

+ terms of lower class = . , , . . . . • (1)? 

where the coefficients are holomorj)hic functions of x and the p's are doubly periodic 
functions of periods a)^ and co^. 

Either then the original equation can be reduced to this form, or at least one of the 



ratios 



So(^0 



Si(x) . 







n(x)' ■ ■ ■ ' n(x) 

equations of the form 



is composed of an additive number of solutions of difference 



/ (.« + Wj) ™ / (x) = p (x) x\i {x I «2) 
/ {x + Wg) -/ («) = p {x) xfj {x I «i) 



But the most general solution of such a pair of equations i 






q [x) (j) (x) + r (x). 
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where q{x) and r{x) are doubly periodic functions of x of periods co^ and co^. And 

therefore one of the ratios —- , , , . , if—-^ must be a function generated from the 

R{x) ii{x) ^ 

function r2{x\a)i^ w^)* 

The original equation therefore either contains the double gamma function 
implicitly among its coefficients, or it is reducible to the form (1). 

Continue our former procedure, and we see that either at least one of the ratios 

-~~4 , . . . -~7- is composed of an additive number of equations of the type 
B(x) ' E(a)) -^ -^ -^^ 

/(a^ + a)^) -f{x) =p{x) Q^[i//(a?|a),)} 

and is therefore generated from the double gamma function^ or the original equation 
is reducible to one in which the ratios of terms of the three highest classes are doubly 
periodic functions of x of periods Wj^ and w^. 

The successive repetitions of the argument are now evident. Ultimately we 
reduce the equation to one in which either all the coefficients are doubly periodic 
functions (which is absurd), or to one in which the last term is generated from the 
double gamma function. 

Thus the proposition is established. The double gamma function cannot satisfy a 
differential equation in which the coefficients are finite combinations of, e.g., 

( 1 ) Rational or irrational algebraic functions of x, 

(2) Simply or doubly periodic functions, 

(3) Simple gamma functions, 

(4) G functions, 

(5) Theta functions, 

or, in fact, of any functions which are not substantially reducible to or compounded 
of the double gamma function itself. 
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